感觉我只会做这题了。。。。。。。
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;
}