多重背包二进制优化的思考

其实说白了我们最朴素的多重背包做法是将有数量限制的相同物品看成多个不同的0-1背包。这样的时间复杂度为O(V*Σn(i)) V为空间容量，n(i)为每种背包的数量限制。如果这样会超时，我们就得考虑更优的拆分方法，由于拆成1太多了，我们考虑拆成二进制数，对于13的数量，我们拆成1，2，4，6（有个6是为了凑数）。 19 我们拆成1，2，4，8，4（最后的4也是为了凑和为19）。经过这个样的拆分我们可以组合出任意的小于等于n(i)的数目（二进制啊，必然可以）。j极大程度缩减了等效为0-1背包时候的数量。 大概可以使时间复杂度缩减为O(V*lg(ΣN(i))；

给个经典题目　POJ 1014 

也是进行6个物品的数量二进制拆分，等效为0-1背包。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int dp[110056];
int v[500];
int w[500];
int main()
{
    int a[7],TT=0;
    while(cin>>a[1]>>a[2]>>a[3]>>a[4]>>a[5]>>a[6])
    {   if(a[1]==0&&a[2]==0&&a[3]==0&&a[4]==0&&a[5]==0&&a[6]==0)break;
        TT++;
         memset(dp,0,sizeof(dp));
        int flag=0;
        int sum=a[1]+a[2]*2+a[3]*3+a[4]*4+a[5]*5+a[6]*6;
        if(sum%2==1) flag=0;
        else
        {
            sum=sum/2;
            int num=1;
            for(int i=1; i<=6; i++)
            {
                int k=1,c=a[i];
                for(k=1; k<=c; k=k*2)
                {
                    v[num]=i*k;
                    w[num]=v[num];
                    num++;
                    c-=k;
                }
                if(c>0)
                {
                    v[num]=i*c;
                    w[num]=v[num];
                    num++;
                }
            }
            num--;
            for(int i=1; i<=num; i++)
            {
                if(flag)break;
                for(int j=sum; j>=v[i]; j--)
                {
                    dp[j]=max(dp[j],dp[j-v[i]]+w[i]);
                    if(dp[sum]==sum)
                    {
                        flag=1;
                        break;
                    }
                }
            }

        }
        cout<<"Collection #"<<TT<<":"<<endl;
        if(flag)cout<<"Can be divided."<<endl;
        else cout<<"Can't be divided."<<endl;
       cout<<endl;

    }
}