POJ 2923 Relocation 装车问题 【状态压缩DP】+【01背包】

#include<iostream>  
#include<algorithm>  
#include<string>  
#include<cstring>  
#include<cstdio>  
#include<cmath>  
#define mem(a,x) memset(a,x,sizeof(a))  
#define inf 1<<30  
using namespace std;  
const double PI = acos(-1.0);  
typedef long long LL;  
int state[1030];  
int dp[1030];  
bool vis[1005];  
int n,v1,v2,tot;  
int c[12];  
bool ok(int x)//判断一种状态是否可行（可以一次运走）  
{  
    int sun = 0;  
    mem(vis,0);  
    vis[0] = 1;  
    for (int i = 0;i < n;i++)  
    {  
        if ((x>>i)&1)  
        {  
            sun += c[i];  
            for (int j = v1;j >= c[i];j--) //这个真的非常巧妙 开始看半天都不懂，自己模拟一遍才懂  
            {                              //比如说此状态有c1、c2、c3,3个体积，第一次操作把体积c1标记为1，  
                  if (vis[j-c[i]])         //第二次操作把c2和c1+c2两种体积标记为1，第三次把c3和前面的组合标记为1，  
                    vis[j] = 1;            //最后这些体积能组合成的所有体积就都被标记成了1  
            }  
        }  
    }  
    if (sun > v1+v2)//装不下  
        return 0;  
    //总体积小不代表一定装得下，拆分成2份要2份都装得下  
    for (int i = 0;i <= v1;i++)  
    {  
        if (vis[i]&&sun-i <= v2)//如果存在（i,sun-i）这样的组合  
            return 1;           //满足i可以被v1装下(前面for循环是对于v1的，vis[i]表示体积i可以被v1装下)，sun-i可以被v2装下  
    }  
    return 0;  
}  
void init()//初始化找到满足条件的状态  
{  
    tot = 0;  
    for (int i = 1;i < (1<<n);i++)  
    {  
        dp[i] = inf;  
        if (ok(i))  
            state[tot++] = i;  
    }  
}  
int main()  
{  
    int T;  
    cin>>T;  
    int oo = 0;  
    while (T--)  
    {  
        cin>>n>>v1>>v2;  
        for (int i = 0;i < n;i++)  
            scanf("%d",&c[i]);  
        init();  
        int V = (1<<n) - 1;//V是n个1的二进制数  
        dp[0] = 0;//没有物品当然是0次运走  
        for (int i = 0;i < tot;i++)  
        {  
            for (int j = V;j >= 0;j--)  
            {  
                if (dp[j] == inf)  
                    continue;    //原版的背包是dp[j] = min(dp[j],dp[j-c[i]]+w[i])  
                                  //但是显然二进制不好表示减，但是可以用|抽象加  
                                  //这就相当于背包改版成dp[j+c[i]] = min(dp[j+c[i]],dp[j] + w[i])  
                if ((j&state[i])==0) //当然2种状态不能有交集  
                {  
                    dp[j|state[i]] = min(dp[j|state[i]] ,dp[j] + 1);  
                }  
  
            }  
        }  
        printf("Scenario #%d:\n%d\n",++oo,dp[V]);  
        if (T) puts("");  
    }  
    return 0;  
}