HDU 5779 Tower Defence

求有多少个n个点的、边权为1的无向图满足从1号点到其他任意一个点的最短路都不等于k，不连通距离为inf。

n,k≤60.

从前有个蒟蒻，连个O(n^4)的dp都写不好，循环分不清南北==

#include<set>
#include<ctime>
#include<queue>
#include<cstdio>
#include<bitset>
#include<cctype>
#include<bitset>
#include<cstdlib>
#include<cassert>
#include<cstring>
#include<iostream>
#include<algorithm>
#define inf (1<<30)
#define INF (1ll<<62)
#define fi first
#define se second
#define rep(x,s,t) for(register int x=s,t_=t;x<t_;++x)
#define per(x,s,t) for(register int x=t-1,s_=s;x>=s_;--x)
#define prt(x) cout<<#x<<":"<<x<<" "
#define prtn(x) cout<<#x<<":"<<x<<endl
#define pb(x) push_back(x)
#define hash asfmaljkg
#define rank asfjhgskjf
#define y1 asggnja
#define y2 slfvm
using namespace std;
typedef long long ll;
typedef pair<int,int> ii;
template<class T>void sc(T &x){
int f=1;char c;x=0;
while(c=getchar(),c<48)if(c=='-')f=-1;
do x=x*10+(c^48);
while(c=getchar(),c>47);
x*=f;
}
template<class T>void nt(T x){
if(!x)return;
nt(x/10);
putchar(x%10+'0');
}
template<class T>void pt(T x){
if(x<0)putchar('-'),x=-x;
if(!x)putchar('0');
else nt(x);
}
template<class T>void ptn(T x){
pt(x);putchar('\n');
}
template<class T>void pts(T x){
pt(x);putchar(' ');
}
template<class T>inline void Max(T &x,T y){if(x<y)x=y;}
template<class T>inline void Min(T &x,T y){if(x>y)x=y;}

const int maxn=65;
ll dp[maxn][maxn][maxn];//dp[i][j][k] 共j个点 第i层有k个点
ll dp2[maxn][maxn];
ll dp3[maxn][maxn];
const int mod=1e9+7;
ll cc[maxn][maxn];
ll tt[maxn*maxn];
ll qpow(ll a,ll b){
ll c=1;
for(;b;b>>=1,a=(a*a)%mod)
if(b&1)c=(a*c)%mod;
return c;
}
int main(){
rep(i,0,maxn){
rep(j,1,i)cc[i][j]=(cc[i-1][j-1]+cc[i-1][j])%mod;
cc[i][0]=cc[i][i]=1;
}
tt[0]=1;
rep(i,1,maxn*maxn)
tt[i]=(tt[i-1]<<1)%mod;
//预处理2的幂

dp[1][1][1]=1;//1层的方案仅一种
for(int i=1;i+1<maxn;i++){//正推
for(int j=i;j<maxn;j++){//枚举前i层有几个
for(int k=1;k<=j-i+1;k++){//枚举....=，=写反了呢
if(!dp[i][j][k])continue;
for(int s=1;j+s<maxn;s++){
dp[i+1][j+s][s]+=dp[i][j][k]*cc[j+s-1][s]%mod*qpow(tt[k]-1,s)%mod*tt[s*(s-1)/2]%mod;
dp[i+1][j+s][s]%=mod;
}
}
}
}
for(int i=0;i<maxn;i++)
for(int j=0;j<maxn;j++)
for(int k=0;k<maxn;k++){
dp2[i][j]+=dp[i][j][k];
dp2[i][j]%=mod;
}

for(int i=1;i<maxn;i++){//层数上限
for(int j=1;j<maxn;j++){//个数上限
for(int k=1;k<=i;k++){//
for(int s=1;s<=j;s++){
dp3[i][j]+=dp2[k][s]*cc[j-1][j-s]%mod*tt[(j-s)*(j-s-1)/2]%mod;
dp3[i][j]%=mod;
}
}
}
}

int cas,a,b;sc(cas);
while(cas--){
sc(a);sc(b);
ptn(dp3[b][a]);
}
return 0;
}

蒟蒻觉得自己很可怜，于是又打了个表。

puts("A[60][60]={");
for(int a=0;a<60;a++){
printf("{");
for(int b=0;b<60;b++){
cout<<dp3[b+1][a+1];
if(b!=59)printf(",");
else printf("}");
}
puts(a!=59?",":"");
}
puts("};");

但因为这是蒟蒻打的第一个二维的表，所以写了博客。