hdu 5713 图去多少条边变为指定连通块 dp ...

图去多少条边变为指定连通块

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<cstdlib>

using namespace std;
int low_bit(int x) {return x&(-x) ;}
const int mod = 1000000009 ;
typedef long long ll ;
int n ;
ll bit[120] , edge[20][20] , cnt[20000] , g[200000];
ll dp[20000][15] , f[20000] ;
int cal(int s , int p){
    int res = 0 ;
    p = lower_bound(bit , bit+15 , p) - bit ;
    for(int i = 0 ; i < n ; i ++ ) if(s & bit[i]) res += edge[i][p] ;
    return res ;
}
void work(int K){
    memset(cnt , 0 , sizeof(cnt)) ; memset(dp , 0 , sizeof(dp)) ;
    memset(f , 0 , sizeof(f)) ; memset(g , 0 , sizeof(g)) ;
    for(int i = 1 ; i < bit[n] ; i ++) cnt[i] = cnt[i - low_bit(i)] + cal(i , low_bit(i)) ;
    for(int i = 1 ; i < bit[n] ; i ++){
        int j = i ;
        do{
            j = j - 1&i ;
            if(j & low_bit(i)) g[i] = (g[i] + bit[cnt[i - j]]*f[j]+mod)%mod ;
        }while(i != j) ;
        f[i] = (bit[cnt[i]] - g[i] + mod ) % mod ;
    }
    dp[0][0] = 1 ;
    for(int i = 1 ; i < bit[n] ; i ++ ){
        for(int k = 1 ; k <= K ; k ++ ){
            int j = i ;
            do{
                j = j - 1&i ;
                if(j & low_bit(i)) dp[i][k] = (dp[i][k] + dp[i - j][k - 1]*f[j]%mod+mod)%mod ;
            }while(i != j) ;
        }
    }
    printf("%lld\n" , (dp[bit[n]-1][K] + mod)%mod) ;
}
int main(){
    bit[0] = 1 ;
    for(int i = 1 ; i <= 115 ; i ++) bit[i] = (bit[i-1]<<1)%mod ;
    int kase ; scanf("%d" , &kase) ; int t = 1 ;
    int m , k , a , b ;
    while(kase --){
        memset(edge , 0 , sizeof(edge)) ;
        scanf("%d %d %d" , &n , &m , &k) ;
        for(int i = 0 ; i < m ; i ++){
            scanf("%d %d", &a , &b) ;
            edge[a - 1][b - 1] ++ ;
            if(a != b) edge[b - 1][a - 1] ++ ;
        }
        printf("Case #%d:\n", t ++ ) ;
        work(k) ;
    }
    return 0 ;
}