硬币购物
【问题描述】
硬币购物一共有
4
种硬币。
面值分别为
某人去商店买东西,去了
tot
次。
每次带
di
枚
ci
硬币,买
si
的价值的东西。
请问每次有多少种付款方法。
【输入描述】
第一行
c1
,
c2
,
c3
,
c4
,
tot
。
下面
tot
行
d1
,
d2
,
d3
,
d4
,
s
, 其中
【输出描述】
每次的方法数。
【输入样例】
1 2 5 10 2
3 2 3 1 10
1000 2 2 2 900
【输出样例】
4
27
【Solution】
参见Byvoid 。
动态规划+容斥原理:
设 fi 为不考虑每种硬币的数量限制的情况下,得到面值 i 的方案数。
则状态转移方程为:
为避免方案重复,要以
k
为阶段递推,边界条件为
接下来对于每次询问,奇妙的解法如下:
根据容斥原理,答案为 得到面值
S
的超过限制的方案数 - 第
当第
1
种硬币超过限制时,只要要用到
当且仅当 Si−(D1+1)⋅C1>=0 ,否则方案数为 0 。
其余情况类似,每次询问只用问
【Code】
- #include <iostream>
- #include <cstdio>
- #define LL long long
- using namespace std;
- LL c[10],tot,ans,tmp;
- LL f[100010];
- int main(){
- scanf(”%lld%lld%lld%lld%lld”,&c[1],&c[2],&c[3],&c[4],&tot);
- f[0]=1;
- for(LL i=1;i<=4;i++)
- for(LL j=c[i];j<=100000;j++)f[j]+=f[j-c[i]];
- for(LL i=1;i<=tot;i++){
- LL d1,d2,d3,d4,s;
- scanf(”%lld%lld%lld%lld%lld”,&d1,&d2,&d3,&d4,&s);
- ans=f[s];
- // ———————————————————-
- tmp=0;
- if(s-(d1+1)*c[1]>=0)tmp+=f[s-(d1+1)*c[1]];
- if(s-(d2+1)*c[2]>=0)tmp+=f[s-(d2+1)*c[2]];
- if(s-(d3+1)*c[3]>=0)tmp+=f[s-(d3+1)*c[3]];
- if(s-(d4+1)*c[4]>=0)tmp+=f[s-(d4+1)*c[4]];
- ans-=tmp;
- // ———————————————————-
- tmp=0;
- if(s-(d1+1)*c[1]-(d2+1)*c[2]>=0)tmp+=f[s-(d1+1)*c[1]-(d2+1)*c[2]];
- if(s-(d1+1)*c[1]-(d3+1)*c[3]>=0)tmp+=f[s-(d1+1)*c[1]-(d3+1)*c[3]];
- if(s-(d1+1)*c[1]-(d4+1)*c[4]>=0)tmp+=f[s-(d1+1)*c[1]-(d4+1)*c[4]];
- if(s-(d2+1)*c[2]-(d3+1)*c[3]>=0)tmp+=f[s-(d2+1)*c[2]-(d3+1)*c[3]];
- if(s-(d2+1)*c[2]-(d4+1)*c[4]>=0)tmp+=f[s-(d2+1)*c[2]-(d4+1)*c[4]];
- if(s-(d3+1)*c[3]-(d4+1)*c[4]>=0)tmp+=f[s-(d3+1)*c[3]-(d4+1)*c[4]];
- ans+=tmp;
- // ———————————————————-
- tmp=0;
- if(s-(d1+1)*c[1]-(d2+1)*c[2]-(d3+1)*c[3]>=0)tmp+=f[s-(d1+1)*c[1]-(d2+1)*c[2]-(d3+1)*c[3]];
- if(s-(d1+1)*c[1]-(d2+1)*c[2]-(d4+1)*c[4]>=0)tmp+=f[s-(d1+1)*c[1]-(d2+1)*c[2]-(d4+1)*c[4]];
- if(s-(d1+1)*c[1]-(d3+1)*c[3]-(d4+1)*c[4]>=0)tmp+=f[s-(d1+1)*c[1]-(d3+1)*c[3]-(d4+1)*c[4]];
- if(s-(d2+1)*c[2]-(d3+1)*c[3]-(d4+1)*c[4]>=0)tmp+=f[s-(d2+1)*c[2]-(d3+1)*c[3]-(d4+1)*c[4]];
- ans-=tmp;
- // ———————————————————-
- tmp=0;
- if(s-(d1+1)*c[1]-(d2+1)*c[2]-(d3+1)*c[3]-(d4+1)*c[4]>=0)tmp+=f[s-(d1+1)*c[1]-(d2+1)*c[2]-(d3+1)*c[3]-(d4+1)*c[4]];
- ans+=tmp;
- // ———————————————————-
- printf(”%lld\n”,ans);
- }
- return 0;
- }