P1450 [HAOI2008]硬币购物

题意:

共有 4 种硬币。面值分别为c1,c2,c3,c4c_1,c_2,c_3,c_4c1​,c2​,c3​,c4​。
某人去商店买东西,去了 n 次,对于每次购买,他带了 did_idi​枚 i 种硬币,想购买 s 的价值的东西。请问每次有多少种付款方法。

题解:

第一感觉是母函数,但写不出来,后来看题解时确实可以做,NTT+多项式求逆+母函数,麻烦而且跑的很慢,不过洛谷可以过,这里不做详细讲解,可以看看代码
母函数做法的代码
我们这里讲下完全背包+容斥的做法
如果题目没有数量的限制,可以完全背包做法
现在有了数量限制,我们可以用差分来求,举个例子:
[2,+∞)−(3,+∞)=[2,+∞)−[4,+∞)=[2,3][2,+∞)−(3,+∞)=[2,+∞)−[4,+∞)=[2,3][2,+∞)−(3,+∞)=[2,+∞)−[4,+∞)=[2,3],显然只要我们求出[2,+∞)[2,+∞)[2,+∞)和(3,+∞)(3,+∞)(3,+∞),就可以算出[2,3]。可以认为[x,+∞)[x,+∞)[x,+∞)表示价值大于等于x的方案数,可以用完全背包去求。
现在我们用完全背包先预处理出没有钱数限制的所有情况,然后减去不合法情况(即某种硬币超出了所给数量),得到的不就是答案:
ans=f[s]−∑i=14f[s−ci∗(di+1)]ans=f[s]-\sum_{i=1}^{4}f[s-c_i*(d_i+1)]ans=f[s]−∑i=14​f[s−ci​∗(di​+1)]
现在的任务就是求∑i=14f[s−ci∗(di+1)]\sum_{i=1}^{4}f[s-c_i*(d_i+1)]∑i=14​f[s−ci​∗(di​+1)],可能会认为直接累加不就可以了,并不是,因为i=1,2,3,4,这四部分并不是完全独立的,换句话说,有可能第一种物品超出限制的同时,第二种也超出了,如果直接累加必然会造成重复计算。这咋整?这个模型很想中学学过的容斥原理:
这个式子应该很熟悉:
三元容斥:

A∪B∪C=(A∪B)∪C=A+B+C−A∩B−A∩C−B∩C+A∩B∩C
扩展一下:
奇加偶减

我们设四种硬币分别是A,B,C,D,card(X)表示X集合元素个数
答案=f[s]=card(A⋃B⋃C⋃D)f[s]=card(A\bigcup B\bigcup C\bigcup D)f[s]=card(A⋃B⋃C⋃D)
=f[s]−ans=f[s]-ans=f[s]−ans
ans=card(A)+card(B)+card(C)+card(D)−card(A⋂B)−card(A⋂C)−card(A⋂D)−card(B⋂C)−card(B⋂D)−card(C⋂D)+card(A⋂B⋂C)+card(A⋂B⋂D)+card(B⋂C⋂D)−card(A⋂B⋂C⋂D)ans=card(A)+card(B)+card(C)+card(D)-card(A\bigcap B)-card(A\bigcap C)-card(A\bigcap D)-card(B\bigcap C)-card(B\bigcap D)-card(C\bigcap D)+card(A\bigcap B\bigcap C)+card(A\bigcap B\bigcap D)+card(B\bigcap C\bigcap D)-card(A\bigcap B\bigcap C\bigcap D)ans=card(A)+card(B)+card(C)+card(D)−card(A⋂B)−card(A⋂C)−card(A⋂D)−card(B⋂C)−card(B⋂D)−card(C⋂D)+card(A⋂B⋂C)+card(A⋂B⋂D)+card(B⋂C⋂D)−card(A⋂B⋂C⋂D)
我们刚才说了奇加偶减,所以二进制枚举(第i位为0表示这个card()种并没有因为i,否则有),直接从0(0000)枚举到15(1111)就可以。
本题难点还是在于将完全背包与容斥结合,两个分开看都不算难

代码:

#include <bits/stdc++.h>
#include <unordered_map>
#define debug(a, b) printf("%s = %d\n", a, b);
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;
clock_t startTime, endTime;
//Fe~Jozky
const ll INF_ll= 1e18;
const int INF_int= 0x3f3f3f3f;
void read(){};
template <typename _Tp, typename... _Tps> void read(_Tp& x, _Tps&... Ar)
{x= 0;char c= getchar();bool flag= 0;while (c < '0' || c > '9')flag|= (c == '-'), c= getchar();while (c >= '0' && c <= '9')x= (x << 3) + (x << 1) + (c ^ 48), c= getchar();if (flag)x= -x;read(Ar...);
}
template <typename T> inline void write(T x)
{if (x < 0) {x= ~(x - 1);putchar('-');}if (x > 9)write(x / 10);putchar(x % 10 + '0');
}
void rd_test()
{#ifdef ONLINE_JUDGE
#elsestartTime = clock ();freopen("data.in", "r", stdin);
#endif
}
void Time_test()
{#ifdef ONLINE_JUDGE
#elseendTime= clock();printf("\nRun Time:%lfs\n", (double)(endTime - startTime) / CLOCKS_PER_SEC);
#endif
}
const int maxn=1e5+9;
ll ans,f[maxn];
int c[5];//面值
int d[5];//数量
void pre(){f[0]=1;for(int i=1;i<=4;i++){for(int j=c[i];j<=100001;j++){f[j]+=f[j-c[i]];}}
}
int s;
void cal(){ans=0;for(int i=1;i<=(1<<4)-1;i++){ll k=0,num=0;for(int j=0;j<4;j++){if(i&(1<<j))//第j位是1{num++;k+=c[j+1]*(d[j+1]+1); } }ll flag=1;if(num%2==0)flag=-1*flag;if(s>=k)ans+=flag*f[s-k];}ans=f[s]-ans;
}
int main()
{//rd_test();for(int i=1;i<=4;i++)cin>>c[i];pre();int t;read(t);while(t--){ans=0;for(int i=1;i<=4;i++)cin>>d[i];cin>>s;cal();cout<<ans<<endl;}//Time_test();
}

P1450 [HAOI2008]硬币购物相关推荐

  1. P1450 [HAOI2008] 硬币购物题解

    题面 题目描述 共有 4 种硬币.面值分别为 ,,,​. 某人去商店买东西,去了 n 次,对于每次购买,他带了 ​ 枚 i 种硬币,想购买 s 的价值的东西.请问每次有多少种付款方法. 输入格式 输入 ...

  2. BZOJ 1042 [HAOI2008]硬币购物

    1042: [HAOI2008]硬币购物 Time Limit: 10 Sec  Memory Limit: 162 MB Submit: 1446  Solved: 845 [Submit][Sta ...

  3. BZOJ 1042: [HAOI2008]硬币购物 [容斥原理]

    1042: [HAOI2008]硬币购物 题意:4种硬币.面值分别为c1,c2,c3,c4.1000次询问每种硬币di个,凑出\(s\le 10^5\)的方案数 完全背包方案数? 询问太多了 看了题解 ...

  4. bzoj 1042: [HAOI2008]硬币购物(dp+容斥)

    1042: [HAOI2008]硬币购物 Time Limit: 10 Sec  Memory Limit: 162 MB Submit: 2555  Solved: 1537 [Submit][St ...

  5. BZOJ1042 [HAOI2008]硬币购物

    Description 硬币购物一共有4种硬币.面值分别为c1,c2,c3,c4.某人去商店买东西,去了tot次.每次带di枚ci硬币,买s i的价值的东西.请问每次有多少种付款方法. Input 第 ...

  6. 1042: [HAOI2008]硬币购物 - BZOJ

    Description 硬币购物一共有4种硬币.面值分别为c1,c2,c3,c4.某人去商店买东西,去了tot次.每次带di枚ci硬币,买si的价值的东西.请问每次有多少种付款方法. Input 第一 ...

  7. [bzoj1042][HAOI2008]硬币购物

    有三种硬币,每种有自己的币值. 然后有n次询问,每次都给出每种硬币的数量和要付的钱s,求有多少种付法.n<=1000 s<=100000 ------ 不考虑限制,就是个简单dp.... ...

  8. 洛谷P1450:硬币购物(背包、容斥)

    解析 呜呜呜不废啊 我只会跑n遍多重背包 感觉非常神仙的一道题 之所以只是蓝的可能是因为代码实现难度太低了吧 但感觉思想真的很难想到 也可能是我太菜了 容斥相关还是需要加强啊qwq 考虑如果没有硬币个 ...

  9. HAOI2008 硬币购物

    题目传送门 据说\(NOIp\)前写题解会\(\mathcal{RP}\)++ 看数据范围,肯定不能写多重背包,会\(T\)飞~ 如果每种硬币没有个数限制,就可以用完全背包了. 正难则反,我们可以先用 ...

最新文章

  1. 图像滤镜艺术---(Lightleaks Filter)漏光滤镜
  2. J2EE详细入门教程--人员登入
  3. 写在NOIP2018后
  4. NGUI 3.5教程(二)Label 标签 (Hello world)、多行文本
  5. python中为什么训练数据_python训练数据时打乱训练数据与标签
  6. json 对象 数组
  7. mysql指定某行或者某列的排序
  8. AspNetCore中使用Ocelot之 IdentityServer4(1)
  9. 【Express】—post传递参数
  10. Ubuntu18.04安装Android Studio
  11. MAC OS Finder 中快速定位指定路径
  12. 递归算法经典实例python-python-动态规划的递归、非递归实现
  13. 计算机设备维修与及日常保养,电脑主机日常的维护保养计划
  14. python词频统计三国演义_python实例:三国演义TXT文本词频分析
  15. 武汉纺织大学计算机科学校区在哪,武汉纺织大学阳光校区在哪 环境好不好
  16. 国内linux内核代码贡献率,[图]AMD为Linux内核贡献27.5万行代码 确认Van Gogh APU支持DDR5和VCN3...
  17. 如何达到高效的网络信息传播
  18. 用AI算法起中文名字 ---- AI起名的评价
  19. zabbix监控深信服_zabbix snmp类型 无需安装agent也能监控(51)
  20. C++Directx11开发笔记二:Direct3D基础设备的初始化

热门文章

  1. 令人惋惜的天才新秀:16岁上剑桥大学,27岁就出名,数学事业一路畅通无阻,但自从结婚后,人生从此翻天覆地······
  2. 与毒”共舞30年!清华美女研究生为何放弃高薪,选择特招入伍?背后的原因令人泪崩......
  3. 计算机试题dddd,数据库系统概论试题及答案dddd_Image_Marked.pdf
  4. 超线程_超线程加核显 i310100+梅捷H410超值爆款组合
  5. java基本类型与引用数据类型_java基本数据类型与引用数据类型总结
  6. c++ 隐藏进程_Linux 查看进程的动态信息
  7. 楚乔传手游 服务器维护,楚乔传手游网络异常进不去游戏怎么办?楚乔传网络异常解决方法...
  8. 经纬度 c代码中定义_如何将TXT文本格式的批量经纬度值导入到奥维成为标签
  9. java locale.us_JAVA实现国际化
  10. php curl 使用方法,php curl使用方法与步骤