hdu 6006 Engineer Assignment（状压）

题目： here

主要思路：离散化之后进行状压dp

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;
typedef long long ll;
const int maxn=1e5+5;
int engi[15][5];
int pro[15][5];
int dp[15][1200];
int vis[105];
int proall[15];
int engiall[15];
vector<int>ve[15];
int main()
{
    int l,T,n,m,t=0;
    scanf("%d",&T);
    while(T--)
    {
        t++;
        memset(vis,0,sizeof(vis));
        memset(dp,0,sizeof(dp));
        memset(engiall,0,sizeof(engiall));
        memset(proall,0,sizeof(proall));
        scanf("%d %d",&n,&m);
        for(int i=1;i<=n;i++)
        {
            ve[i].clear();
            scanf("%d",&engi[i][0]);
            for(int j=1;j<=engi[i][0];j++)
            {
                scanf("%d",&engi[i][j]);//无视数组名字，这是工程
            }
        }
         for(int i=1;i<=m;i++)
        {
            scanf("%d",&pro[i][0]);
            for(int j=1;j<=pro[i][0];j++)
            {
                scanf("%d",&pro[i][j]);//无视数组名字，读入工程师
                vis[pro[i][j]]=1;      //标记值，用于离散化
            }
        }
        int flag=0;
        for(int i=1;i<=100;i++)
        {
            if(vis[i])
            {
                flag++;
                vis[i]=flag;         //将值进行编号，用于状态压缩
            }
        }
        for(int i=1;i<=m;i++)
        {
            for(int j=1;j<=pro[i][0];j++)
            {
                proall[i]+=(1<<(vis[pro[i][j]]-1));    //将工程师的expert areas在离散化之后 压缩
            }
        }
        for(int i=1;i<=n;i++)                                 //找出每个工程 能完成的 工程师组合 的种类
        {
            for(int j=1;j<=engi[i][0];j++)
            {
                if(vis[engi[i][j]])
                   engiall[i]+=(1<<(vis[engi[i][j]]-1));   //将工程师的expert areas在离散化之后 压缩
                else
                {
                    engiall[i]=-1;                       //如果有工程师的没有的expert areas，那么这项工程肯定无法完成
                    break;
                }
            }
            if(engiall[i]!=-1)
            {
                for(int j=0;j<(1<<m);j++)               
                {
                    int now=0;
                    int mem3=0;
                    for(int k=0;k<30;k++)
                    {
                       if(mem3>3) break;               //一个工程最多只要三个工程师
                       if(j&(1<<k))
                       {
                           now|=proall[k+1];
                           mem3++;
                       }
                    }
                    if(mem3>3) continue;
                    if((now&engiall[i])==engiall[i])
                    {
                        ve[i].push_back(j);               //将一个工程能对应完成的工程师的状态压缩成数放入
                    }
                }
            }
        }
        for(int i=1;i<=n;i++)                           //dp 第一维是指 前i个工程，第二维是 工程师的使用状态，dp值是完成的工程数
            for(int j=0;j<(1<<m);j++)
            {
                for(int k=0;k<ve[i].size();k++)
                {
                    if(!(j&ve[i][k]))
                    {
                        dp[i][j|ve[i][k]]=max(dp[i][j|ve[i][k]],dp[i-1][j]+1);
                    }
                }
                dp[i][j]=max(dp[i][j],dp[i-1][j]);
            }
        int ans=0;
        for(int j=0;j<(1<<m);j++)
        {
            ans=max(ans,dp[n][j]);
        }
      printf("Case #%d: %d\n",t,ans);
    }
    return 0;
}