P1450 [HAOI2008]硬币购物
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=14f[s−ci∗(di+1)]
现在的任务就是求∑i=14f[s−ci∗(di+1)]\sum_{i=1}^{4}f[s-c_i*(d_i+1)]∑i=14f[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]硬币购物相关推荐
- P1450 [HAOI2008] 硬币购物题解
题面 题目描述 共有 4 种硬币.面值分别为 ,,,. 某人去商店买东西,去了 n 次,对于每次购买,他带了 枚 i 种硬币,想购买 s 的价值的东西.请问每次有多少种付款方法. 输入格式 输入 ...
- BZOJ 1042 [HAOI2008]硬币购物
1042: [HAOI2008]硬币购物 Time Limit: 10 Sec Memory Limit: 162 MB Submit: 1446 Solved: 845 [Submit][Sta ...
- BZOJ 1042: [HAOI2008]硬币购物 [容斥原理]
1042: [HAOI2008]硬币购物 题意:4种硬币.面值分别为c1,c2,c3,c4.1000次询问每种硬币di个,凑出\(s\le 10^5\)的方案数 完全背包方案数? 询问太多了 看了题解 ...
- bzoj 1042: [HAOI2008]硬币购物(dp+容斥)
1042: [HAOI2008]硬币购物 Time Limit: 10 Sec Memory Limit: 162 MB Submit: 2555 Solved: 1537 [Submit][St ...
- BZOJ1042 [HAOI2008]硬币购物
Description 硬币购物一共有4种硬币.面值分别为c1,c2,c3,c4.某人去商店买东西,去了tot次.每次带di枚ci硬币,买s i的价值的东西.请问每次有多少种付款方法. Input 第 ...
- 1042: [HAOI2008]硬币购物 - BZOJ
Description 硬币购物一共有4种硬币.面值分别为c1,c2,c3,c4.某人去商店买东西,去了tot次.每次带di枚ci硬币,买si的价值的东西.请问每次有多少种付款方法. Input 第一 ...
- [bzoj1042][HAOI2008]硬币购物
有三种硬币,每种有自己的币值. 然后有n次询问,每次都给出每种硬币的数量和要付的钱s,求有多少种付法.n<=1000 s<=100000 ------ 不考虑限制,就是个简单dp.... ...
- 洛谷P1450:硬币购物(背包、容斥)
解析 呜呜呜不废啊 我只会跑n遍多重背包 感觉非常神仙的一道题 之所以只是蓝的可能是因为代码实现难度太低了吧 但感觉思想真的很难想到 也可能是我太菜了 容斥相关还是需要加强啊qwq 考虑如果没有硬币个 ...
- HAOI2008 硬币购物
题目传送门 据说\(NOIp\)前写题解会\(\mathcal{RP}\)++ 看数据范围,肯定不能写多重背包,会\(T\)飞~ 如果每种硬币没有个数限制,就可以用完全背包了. 正难则反,我们可以先用 ...
最新文章
- 图像滤镜艺术---(Lightleaks Filter)漏光滤镜
- J2EE详细入门教程--人员登入
- 写在NOIP2018后
- NGUI 3.5教程(二)Label 标签 (Hello world)、多行文本
- python中为什么训练数据_python训练数据时打乱训练数据与标签
- json 对象 数组
- mysql指定某行或者某列的排序
- AspNetCore中使用Ocelot之 IdentityServer4(1)
- 【Express】—post传递参数
- Ubuntu18.04安装Android Studio
- MAC OS Finder 中快速定位指定路径
- 递归算法经典实例python-python-动态规划的递归、非递归实现
- 计算机设备维修与及日常保养,电脑主机日常的维护保养计划
- python词频统计三国演义_python实例:三国演义TXT文本词频分析
- 武汉纺织大学计算机科学校区在哪,武汉纺织大学阳光校区在哪 环境好不好
- 国内linux内核代码贡献率,[图]AMD为Linux内核贡献27.5万行代码 确认Van Gogh APU支持DDR5和VCN3...
- 如何达到高效的网络信息传播
- 用AI算法起中文名字 ---- AI起名的评价
- zabbix监控深信服_zabbix snmp类型 无需安装agent也能监控(51)
- C++Directx11开发笔记二:Direct3D基础设备的初始化
热门文章
- 令人惋惜的天才新秀:16岁上剑桥大学,27岁就出名,数学事业一路畅通无阻,但自从结婚后,人生从此翻天覆地······
- 与毒”共舞30年!清华美女研究生为何放弃高薪,选择特招入伍?背后的原因令人泪崩......
- 计算机试题dddd,数据库系统概论试题及答案dddd_Image_Marked.pdf
- 超线程_超线程加核显 i310100+梅捷H410超值爆款组合
- java基本类型与引用数据类型_java基本数据类型与引用数据类型总结
- c++ 隐藏进程_Linux 查看进程的动态信息
- 楚乔传手游 服务器维护,楚乔传手游网络异常进不去游戏怎么办?楚乔传网络异常解决方法...
- 经纬度 c代码中定义_如何将TXT文本格式的批量经纬度值导入到奥维成为标签
- java locale.us_JAVA实现国际化
- php curl 使用方法,php curl使用方法与步骤