[多重背包+二进制优化]HDU1059 Dividing

 1 #include <cstdio>
 2 #include <algorithm>
 3 #include <cstring>
 4 int dp[200005],n[7],V[105],W[100];
 5 
 6 int main(){
 7     int k=0;
 8     while(scanf("%d%d%d%d%d%d",&n[1],&n[2],&n[3],&n[4],&n[5],&n[6])!=EOF && (n[1]||n[2]||n[3]||n[4]||n[5]||n[6])){
 9         printf("Collection #%d:\n",++k);
10         int i,j,sum=0,mid;
11         for(i=1;i<=6;i++) sum+=i*n[i];
12         if(sum%2){
13             printf("Can‘t be divided.\n\n");
14             continue;
15         }
16         else mid = sum/2;
17         memset(dp,0,sizeof(dp));  //切记不要忘了初始化 
18         int count=0,temp,flag=1;
19         //这块是二进制优化 
20         for(i=1;i<=6;i++){
21             temp=1;
22             while(n[i]>=temp){
23                 V[++count]=i*temp;
24                 W[count]=i*temp;
25                 n[i]-=temp;
26                 temp*=2;
27             }
28             if(n[i]>0){
29                 V[++count]=i*n[i];
30                 W[count]=i*n[i];
31             }
32         }
33         //上面是二进制优化 
34         for(int i=1;i<=count && flag;i++)
35             for(int j=mid;j>=V[i] && flag;j--){
36                 dp[j] = std::max(dp[j],dp[j-V[i]]+W[i]);
37                 if(dp[j]==mid){   
38                     printf("Can be divided.\n\n");
39                     flag=0;
40                 }
41             }
42         if(flag){
43             printf("Can‘t be divided.\n\n");
44         }
45     }
46     return 0;
47 }