hdu 5713(状态压缩DP)

要进行两次dp，

第一个，dp[i],1<=i<=(1<<n)

其中用i的二进制形式表示已选择的点。

dp[i] 用来保存i中的点构成一个连通块，边集多少种可能。

转移方程：

        save[] = ;//这里用save[i]表示dp[i]
        for(int i=;i<(<<n);i++)
        {
            int ti = i-lowbit(i); //一定选择最后一个点，使之有序
            int j = ti;
            long long tmp=;//问题转化为对立问题，记录有多少不连通边集的情况，all[i]表示i中的点，边集的所有可能情况。
            for(;j>;j=(j-)&ti)
            {
                tmp += save[i-j]*all[j];
                tmp %= MOD;
            }
            save[i] = (all[i]-tmp+MOD)%MOD;
        }    

第二次dp，套路基本相同。

dp[i][j] 表示恰有（i+1）个连通块，且含j中的点的所有可能。

转移方程：

　　　　　for(int i=;i<kk;i++)
        {
            for(int j=;j<(<<n);j++)
            {
                int tj= j-lowbit(j);
                for(int k=tj; k>;k = (k-)&tj )
                {
                    dp[i][j] += dp[i-][k]*save[j-k];
                    dp[i][j] %= MOD;
                }
            }
        }

//
//  main.cpp
//  Astar160529
//
//  Created by 陈加寿 on 16/5/29.
//  Copyright © 2016年 chenhuan001. All rights reserved.
//

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <algorithm>
using namespace std;

#define MOD 1000000009

long long dp[][<<];
int mat[][];
long long save[<<];
long long s2[];
long long all[<<];

int lowbit(int x)
{
    return x&(-x);
}

void test()
{
    int cnt=;
    for(int i=;i<(<<);i++)
    {
        for(int j=i;j>;j=(j-)&i)
        {
            cnt++;
        }
    }
    cout<<cnt<<endl;
    cout<<(<<)<<endl;
}

int main() {

    s2[] = ;
    for(int i=;i <= *;i++)
        s2[i] = (*s2[i-])%MOD;

    //test();
    int T;
    cin>>T;

    int tt=;
    while(T--)
    {
        int kk,n,m;
        scanf("%d%d%d",&n,&m,&kk);
        int cnt0=;
        memset(mat,,sizeof(mat));
        for(int i=;i<m;i++)
        {
            int a,b;
            scanf("%d%d",&a,&b);
            a--; b--;
            if(a==b) cnt0++;
            else
            {
                mat[a][b] = mat[b][a] = ;
            }
        }

        memset(save,,sizeof(save));
        memset(all,,sizeof(all));
        memset(dp,,sizeof(dp));

        for(int i=;i<(<<n);i++)
        {
            int tcnt=;
            for(int j=;j<n;j++)
            {
                if( ((<<j)&i) !=)
                {
                    for(int k=j+;k<n;k++)
                    {
                        if( ((<<k)&i) != )
                        {
                            tcnt += mat[j][k];
                        }
                    }
                }
            }
            all[i] = s2[tcnt];
        }

        save[] = ;
        for(int i=;i<(<<n);i++)
        {
            int ti = i-lowbit(i);
            int j = ti;
            long long tmp=;
            for(;j>;j=(j-)&ti)
            {
                tmp += save[i-j]*all[j];
                tmp %= MOD;
            }
            save[i] = (all[i]-tmp+MOD)%MOD;
            dp[][i] = save[i];
        }

        for(int i=;i<kk;i++)
        {
            for(int j=;j<(<<n);j++)
            {

                int tj= j-lowbit(j);
                for(int k=tj; k>;k = (k-)&tj )
                {
                    dp[i][j] += dp[i-][k]*save[j-k];
                    dp[i][j] %= MOD;
                }
            }
        }
        printf("Case #%d:\n",tt++);

        cout<<dp[kk-][(<<n)-]*s2[cnt0]%MOD<<endl;

    }

    return ;
}

CODE