dp题目
首先很容易想到用d[s ][ j ] 代表集合s内的点组成j和联通集合的方式数目,转移为d[ s ][ j ] = (d[ s - s1 ][ j - 1] * f[s1]) 其中s1必须包含s中编号最小的一个元素,这是做限定,不然不好转移,f[s1 ]代表结合s1构成独立联通集合的所有方法数目。
下面求f[ s 1], 计cnt[ s1 ]为s1内点互相之间所有连边的数目, 那么
f[ s1 ] = (1<<cnt[s1] ) - sum((1<<cnt[s1 - s2 ] )* f[s2] ) 其中f[s2]必须包含s1中编号最小的边
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;
typedef long long ll;
#define rep(i,n) for(int i = 0 ; i<(int)n;i++)
#define rep1(i,x,y) for(int i=(int)x;i<=(int)y;i++)
#define lowbit(x) (x&-x)
const int mod = 1000000009;
const int N = 14;
int d[1<<N][N];
int ma[N][N];
int n,m,K,cnt[1<<N];
int wei(int x){
int ans = 0;
while(x > 0){
x>>=1;
ans++;
}
return ans - 1;
}
int mi[N * N + N];
int f[1<<N],g[1<<N];
int main()
{
mi[0] = 1;
rep1(i , 1 , N*N) mi[i] = mi[i - 1] * 2 % mod;
int T;
int kase = 1;
scanf("%d",&T);
while(T--){
scanf("%d %d %d",&n,&m,&K);
memset(ma , 0, sizeof(ma));
rep(i , m){
int x , y;
scanf("%d %d",&x,&y);
--x; --y;
ma[x][y] = ma[y][x] = 1;
}
cnt[0] = 0;
for(int i =1; i < (1<<n); i++){
cnt[i] = cnt[i-lowbit(i)];
int k = wei(lowbit(i));
rep(j , n){
if((i>>j)&1)
cnt[i] += ma[j][k];
}
}
f[0] = 0; g[0] = 1;
for(int s = 1; s < (1<<n) ; s++){
g[s] = 0;
for(int j =(s - 1) & s ; j != 0 ; j = (j - 1)& s) if(j&(lowbit(s))){
g[s] = (g[s] + (ll)g[s - j] * f[j] % mod) % mod;
}
f[s] = (mi[cnt[s]] - g[s] + mod) % mod;
g[s] = (g[s] + f[s]) % mod;
}
d[0][0] = 1;
for(int s = 1; s<(1<<n) ; s++){
for(int k = 0; k<=K;k++){
d[s][k] = 0;
if(k == 0) continue;
for(int j = s ; j ; j=(j-1)&s)if(j&(lowbit(s))){
d[s][k] = (d[s][k] + (ll)d[s - j][k - 1] * f[j] % mod) % mod;
}
}
}
printf("Case #%d:\n",kase++);
printf("%d\n",(d[(1<<n)-1][K] + mod) %mod);
}
return 0;
}