[背包dp+容斥定理] luogu 1450 硬币购物

kosf_

于 2020-10-14 23:28:20 发布

阅读量125

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/kosf_/article/details/109087481

版权

数论专栏收录该内容

19 篇文章 0 订阅

订阅专栏

这篇博客主要探讨了一道编程竞赛题目，涉及到动态规划和完全背包问题的解决。作者通过预处理不同面额硬币的背包问题，并利用二进制枚举处理大量询问，避免超时。在每个询问中，通过容斥原理计算出所有可能的组合，得出最终答案。文章还展示了完整的C++代码实现。

摘要由CSDN通过智能技术生成

题目

在这里插入图片描述
题目链接：https://www.luogu.com.cn/problem/P1450

思路

如果只询问次数很少可以完全背包
不过这个询问次数很多要超时
那么就预处理 $1 — 1 e 5$ 的背包的值
每次带了d[i]枚硬币那么d[i]+1之后的硬币就不行
$s - c [i] * (d [i] + 1)$ 就是第i个硬币能取的最大钱钱
$d p [s - c [i] * (d [i] + 1)]$ 就是第i个硬币的可能性
不过这道题还要考虑一下容斥的东西
在这里插入图片描述
用二进制枚举一下每种情况。

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<cctype>
#include<ctime>
#include<iostream>
#include<string>
#include<map>
#include<queue>
#include<stack>
#include<set>
#include<vector>
#include<iomanip>
#include<list>
#include<bitset>
#include<sstream>
#include<fstream>
#include<complex>
#include<algorithm>
#if __cplusplus >= 201103L
#include <unordered_map>
#include <unordered_set>
#endif
#define int long long
using namespace std;
const int INF = 0x3f3f3f3f;
int c[10],d[10];
int dp[100010];
signed main(){
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
	int n;
	cin>>c[1]>>c[2]>>c[3]>>c[4]>>n;
	dp[0]=1;
	for(int i=1;i<=4;i++){
		for(int j=c[i];j<=100000;j++){
			dp[j]+=dp[j-c[i]];
		}
	}
	while(n--){
		for(int i=1;i<=4;i++){
			cin>>d[i];
		}
		int s;
		cin>>s;
		int ans=dp[s];
		for(int i=1;i<=15;i++){
			int k=0,num=0;
			for(int j=0;j<4;j++){
				if(i&(1<<j)){
					num++;
					k+=c[j+1]*(d[j+1]+1);
				}
			}
				int f=1;
				if(num%2) f=-1;
				if(s>=k) ans+=f*dp[s-k];
		}
		cout<<ans<<endl;
	}
    return 0;
}