hdoj 5779 Tower Defence ??dp

dp[i][j][k]连通块 大小为i 最大长度为j 有k个点的长度为j的方案数
每加入k个长度为j的点 都要与前面的j-1的点相连
这k个点是从i个点里挑出来的
这k个点互相的直接相连是随机的
不在这连通块的各个点相连是随机的
然后推个公式乱七八糟的统计就行了
注意别溢出

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<iostream>
#define maxn 65
#define LL long long
#define mod 1000000007ll
using namespace std;
LL dp[maxn][maxn][maxn] = {0},c[maxn][maxn],mi[maxn*maxn],cnt[maxn][maxn];
void init()
{
    for(int i=0;i<maxn;i++)
    {
        c[i][0] = c[i][i] = 1;
        for(int j=1;j<i;j++)
            c[i][j] = (c[i-1][j-1]+c[i-1][j])%mod;
    }
    mi[0] = 1;
    for(int i=1;i<=maxn*(maxn-1)/2;i++)mi[i] = (mi[i-1]*2)%mod;
    for(int i=1;i<maxn;i++)
    {
        cnt[i][0] = 1;
        for(int j=1;j<maxn;j++)
            cnt[i][j] = cnt[i][j-1]*(mi[i]-1)%mod;
    }
    dp[1][0][1] = 1;//总数1,距离0,个数为1
    for(int i=2;i<maxn;i++)//总数
        for(int j=1;j<i;j++)//距离
            for(int k=1;k<=i;k++)//个数
                for(int u=0;u+k<=i;u++)
                {
                    dp[i][j][k] += dp[i-k][j-1][u]*c[i-1][k]%mod*cnt[u][k]%mod*mi[(k-1)*k/2]%mod;
                    dp[i][j][k] %= mod;
                }
}
int main()
{
    init();
    int t;
    scanf("%d",&t);
    while(t--)
    {
       LL n,k1,ans = 0;
       scanf("%I64d %I64d",&n,&k1);
       for(int i=1;i<=n;i++)//总个数
           for(int j=0;j<k1;j++)//最大长度
               for(int k=1;k<=i;k++)//最大个数
               {
                   ans+=dp[i][j][k]*c[n-1][i-1]%mod*mi[(n-i-1)*(n-i)/2]%mod;
                   ans%=mod;
               }
        printf("%I64d\n",ans);
    }
    return 0;
}