poj 2411 状态压缩dp

最新推荐文章于 2018-08-21 21:27:13 发布

转载最新推荐文章于 2018-08-21 21:27:13 发布 · 288 阅读

状压dp 专栏收录该内容

11 篇文章

订阅专栏

本文提供了一道来自POJ平台的题目2441的解题思路及完整代码实现，该题涉及到状态压缩动态规划，通过预处理状态转移来求解最优解。

转自：http://www.cnblogs.com/chasetheexcellence/archive/2012/04/16/poj2441.html

他还推荐了一个论文。。

反正我是不会敲的。。

他给出了之后我还觉得这个是有点像是DAG。。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<string>
#include<cstring>
#include<iomanip>
#include<iostream>
#include<stack>
#include<cmath>
#include<map>
#include<vector>
#define ll long long
#define inf 0x3f3f3f3f
#define INF 1e9
#define bug1 cout<<"bug1"<<endl;
#define bug2 cout<<"bug2"<<endl;
#define bug3 cout<<"bug3"<<endl;
#include<sstream>
using namespace std;
#define sf scanf
#define pf printf
#define mem(a,b) memset(a,b,sizeof(a));
int n,m;
int bull[25][25];
int node[1<<25];
int dp[1<<25];
int main(){
    while(~sf("%d%d",&n,&m)){
        for(int i=1,tmp;i<=n;++i){
            sf("%d",&tmp);
            bull[i][0]=tmp;
            for(int j=1;j<=bull[i][0];++j){
                sf("%d",&bull[i][j]);
                bull[i][j]--;
            }
        }
        int M=1<<m;
        for(int i=0;i<M;++i){
            int num=0;
            for(int j=0;j<m;++j){
                if(i&(1<<j))num++;
            }
            node[i]=num;

        }
        for(int i=1;i<=bull[1][0];++i){
            int x=bull[1][i];
            dp[1<<x]=1;
        }
        for(int i=2;i<=n;++i){
            for(int j=M-1;j>=0;j--){
                if(node[j]!=i)continue;
                for(int k=1;k<=bull[i][0];++k){
                    int  x=bull[i][k];
                    if(!(j&(1<<x)))continue;
                    dp[j]+=dp[j-(1<<x)];
                }
            }
        }
        int ans=0;
        for(int i=0;i<M;++i){
            if(node[i]==n){
                ans+=dp[i];
            }
        }
        pf("%d\n",ans);

    }
}