bzoj 1042: [HAOI2008]硬币购物

2016-06-20

完全背包DP(做的时候都忘记自己当初怎么学得了  囧。。),题目非同寻常的地方是n很小,价值固定,但有多组测试数据,我们首先想到将之按完全背包预处理出来,每一组时再减去不合法的,

减去不合法的可以用容斥原理,第一个不合法的方案数就是让第一个先取d1+1个其余随便取   f[s-(d1+1)*c1],那就减去1,2,3,4加上12,13,14,23,24,34减去123,124,134,234加上1234。

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<iostream>
 4 #include<cstdlib>
 5 #include<algorithm>
 6 #include<queue>
 7 #define M 100008
 8 #define ll long long
 9 #define Mo 12345678
10 using namespace std;
11 ll read()
12 {
13     char ch=getchar();
14     ll x=0,f=1;
15     for(;ch<'0'||ch>'9';ch=getchar())
16       if(ch=='-')
17         f=-1;
18     for(;ch>='0'&&ch<='9';ch=getchar())
19       x=x*10+ch-'0';
20     return x*f;
21 }
22 int c[5],n,m,a[5],b[5];
23 ll f[M];
24 int main()
25 {
26     f[0]=1;
27     for(int i=1;i<=4;i++)
28       c[i]=read();
29     n=read();
30     for(int i=1;i<=4;i++)
31       for(int j=c[i];j<=100000;j++)
32         f[j]+=f[j-c[i]];
33     for(int i=1;i<=n;i++)
34       {
35           for(int j=1;j<=4;j++)
36             a[j]=read();
37           m=read();
38           ll an=f[m];
39         for(int j=1;j<=4;j++)
40             if(m-(a[j]+1)*c[j]>=0)
41               an-=f[m-(a[j]+1)*c[j]];
42           for(int j=1;j<=4;j++)
43             for(int k=j+1;k<=4;k++)
44               if(m-(a[j]+1)*c[j]-(a[k]+1)*c[k]>=0)
45                 an+=f[m-(a[j]+1)*c[j]-(a[k]+1)*c[k]];
46         for(int j=1;j<=4;j++)
47             for(int k=j+1;k<=4;k++)
48               for(int l=k+1;l<=4;l++)
49               if(m-(a[l]+1)*c[l]-(a[j]+1)*c[j]-(a[k]+1)*c[k]>=0)
50                 an-=f[m-(a[l]+1)*c[l]-(a[j]+1)*c[j]-(a[k]+1)*c[k]];
51         if(m-(a[1]+1)*c[1]-(a[2]+1)*c[2]-(a[3]+1)*c[3]-(a[4]+1)*c[4]>=0)
52           an+=f[m-(a[1]+1)*c[1]-(a[2]+1)*c[2]-(a[3]+1)*c[3]-(a[4]+1)*c[4]];
53         printf("%lld\n",an);
54       }
55     return 0;
56 }

 

posted @ 2016-06-20 15:23  xiw5  阅读(...)  评论(... 编辑 收藏