第八届蓝桥杯省赛C++A/B组 包子凑数（完全背包）

题目链接：包子凑数

---

大致题意:

有几个物品，每个物品无限个，每个物品选任意个，能否凑到某个重量。

问一共有多少数字是凑不出来的，如果有无限个，输出INF

---

解题思路:

如果有无限个凑不出来，说明这些数gcd不是1
裴蜀定理 ，任意两个数的组合必定是他们gcd的倍数，同样可以推广到更多数：如果这些数的gcd是d，那么他们的组合是d的倍数，如果d不是1，那么必然有无限个数无法被组合出来
当gcd为1时，最大不能表示出来的 数是：(a−1)(b−1)−1，而99和98是100内最大的互质的数，所以这个上界选择10000

dp(i,j)表示前选i项物品任意个，重量为j，属性为能否达到重量j
集合分析法：由最后不同的一步来划分集合，在重量为j的时候，第i个物品选了几件
dp(i,j)=dp(i−1,j)||dp(i−1,j−w(i))||…||dp(i−1,j−k∗w(i))(k=j/w(i)) (k=j/w(i))
dp(i,j−w(i))=dp(i−1,j−w(i))||dp(i−1,j−2∗w(i))…||dp(i−1,j−k∗w(i))(k=j/w(i))
dp(i,j−w(i))的状态就是dp(i,j)后面一大部分
最终方程为:dp(i,j)=dp(i−1,j)||dp(i,j−w(i))(w(i)<=j)
初始化：dp(0,0)=true 0件物品，0重量的时候肯定是可以凑出来的。
优化：由于转移方程是由前一个状态转移过来的，所以状态数组可以压成一维:dp(j)=dp(j)||dp(j−w(i))(j>=w(i))dp(j)=dp(j)||dp(j−w(i))(j>=w(i))

---

AC代码:

#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;
const int N = 10010;
int f[110][N];
int a[110];
int gcd(int a, int b) {
	return b == 0 ? a : gcd(b, a % b);
}
int main() {
	int n; cin >> n;
	int d = 0;
	for (int i = 1; i <= n; ++i) {
		cin >> a[i];
		d = gcd(d, a[i]);
	}
	if (d != 1)cout << "INF" << endl;
	else {
		f[0][0] = 1;
		for (int i = 1; i <= n; ++i) {
			for (int j = 0; j <= 10000; ++j) {
				f[i][j] = f[i - 1][j];
				if (j >= a[i])f[i][j] |= f[i][j - a[i]];
			}
		}
		int res = 0;
		for (int i = 0; i <= 10000; ++i)
			if (!f[n][i])res++;
		cout << res << endl;
	}
	return 0;
}

//优化
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long ll;
const int N = 10010;
int f[N];
int a[110];
int gcd(int a, int b) {
	return b == 0 ? a : gcd(b, a % b);
}
int main() {
	int n; cin >> n;
	int d = 0;
	for (int i = 1; i <= n; ++i) {
		cin >> a[i];
		d = gcd(d, a[i]);
	}
	if (d != 1)cout << "INF" << endl;
	else {
		f[0] = 1;
		for (int i = 1; i <= n; ++i) {
			for (int j = a[i]; j <= 10000; ++j) {
				f[j] |= f[j - a[i]];
			}
		}
		int res = 0;
		for (int i = 0; i <= 10000; ++i)
			if (!f[i])res++;
		cout << res << endl;
	}
	return 0;
}

END