题目链接:包子凑数
大致题意:
有几个物品,每个物品无限个,每个物品选任意个,能否凑到某个重量。
问一共有多少数字是凑不出来的,如果有无限个,输出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;
}