hihoCoder 1147 时空阵

感觉我只会做这题了。。。。。。。

Clj给的题解只有一句话QAQ，Holy High

只能自己YY了。

由于K<=100，所以构造分层图。

于是乎边只能和上一层以及本层的节点相连。

设状态(i,j,k)为处理到第i层，共处理了j个结点，本层有k个结点。

于是就可以从(i-1,j-k,x)转移过来，转移的时候k个结点有在本层中可以有k*(k-1)/2条边，上一层点集的非空子集有2^x-1个，本层每个结点与上一层各子集中所有的点都连边，这里的方案就是(2^x-1)^k。

于是有f[i][j][k]=sigma(f[i-1][j-k][x]*(2^x-1)^k*(2^(k*(k-1))/2)*C[n-j+k-1][k]。

然后最后一层特判一下，因为最后一层必须有结点n。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int N=100+5;
const int mod=1000000007;
ll f[N][N][N],C[N][N],bin[N*N];
ll mul(ll a,int n){
	if(!n)return 1;
	ll tmp=mul(a,n/2);
	tmp=(tmp*tmp)%mod;
	if(n&1)tmp=tmp*a%mod;
	return tmp;
}
void build(){
	C[0][0]=1;
	for(int i=1;i<=100;i++){
		C[i][0]=1;
		for(int j=1;j<=i;j++)
		C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
	}
	bin[0]=1;
	for(int i=1;i<=10000;i++){
		bin[i]=(bin[i-1]<<1)%mod;
	}
}
int main(){
	build();
	int n,K;scanf("%d%d",&n,&K);
	f[0][1][1]=1;
	for(int i=1;i<=K;i++)
	for(int j=i;j<=n-K+i;j++)
	for(int k=1;k<=j-i;k++){
		ll t=i==K?C[n-j+k-1][k-1]:C[n-j+k-1][k];
		for(int x=1;x<=j-k-i+1;x++){
			f[i][j][k]+=f[i-1][j-k][x]*t%mod*mul((bin[x]-1+mod)%mod,k)%mod*bin[(k-1)*k/2]%mod;
			f[i][j][k]%=mod;
		}
	}
	ll ans=0;
	for(int i=K+1;i<=n;i++)
	for(int j=1;j<=i-K;j++){
		ans+=f[K][i][j]*bin[j*(n-i)+(n-i)*(n-i-1)/2]%mod;
		ans%=mod;
	}
	printf("%lld",ans);
	return 0;
}