1. 程式人生 > >Dividing POJ - 1014 多重揹包二進位制優化

Dividing POJ - 1014 多重揹包二進位制優化

多重揹包模型  寫的時候漏了一個等號找了半天 i<<=1 !!!!!!

 

 

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 using namespace std;
 5 int a[500];
 6 int dp[120006];
 7 int m[120006];
 8 int main(){
 9 int n=6;
10 int kase=1;
11 while(scanf("%d",&a[0])==1){
12     int
sum=a[0]; 13 int flag=0; 14 for(int j=1;j<=a[0];j<<1){ 15 m[flag++]=j; 16 a[0]-=j; 17 } 18 if(a[0])m[flag++]=a[0]; 19 for(int i=1;i<n;i++){ 20 scanf("%d",&a[i]); 21 sum+=a[i]; 22 for(int j=1;j<=a[i];j<<=1){ 23
m[flag++]=j*(i+1); 24 a[i]-=j; 25 } 26 if(a[i])m[flag++]=a[i]*(i+1); 27 } 28 29 memset(dp,0,sizeof(dp)); 30 if(sum==0)break; 31 sum=0; 32 for(int i=0;i<flag;i++)sum+=m[i]; 33 if(sum&1){ 34 printf("Collection #%d:\nCan't be divided.\n\n",kase++);
35 continue; 36 37 } 38 for(int i=0;i<flag;i++){ 39 for(int j=sum/2;j>=m[i];j--){ 40 dp[j]=max(dp[j],dp[j-m[i]]+m[i]); 41 } 42 } 43 if(dp[sum/2]==sum/2)printf("Collection #%d:\nCan be divided.",kase++); 44 else printf("Collection #%d:\nCan't be divided.",kase++); 45 puts(""); 46 puts(""); 47 48 } 49 return 0; 50 }