蓝桥杯第八届B组包子凑数

标题：包子凑数

小明几乎每天早晨都会在一家包子铺吃早餐。他发现这家包子铺有N种蒸笼，其中第i种蒸笼恰好能放Ai个包子。每种蒸笼都有非常多笼，可以认为是无限笼。

每当有顾客想买X个包子，卖包子的大叔就会迅速选出若干笼包子来，使得这若干笼中恰好一共有X个包子。比如一共有3种蒸笼，分别能放3、4和5个包子。
当顾客想买11个包子时，大叔就会选2笼3个的再加1笼5个的（也可能选出1笼3个的再加2笼4个的）。

当然有时包子大叔无论如何也凑不出顾客想买的数量。比如一共有3种蒸笼，分别能放4、5和6个包子。而顾客想买7个包子时，大叔就凑不出来了。

小明想知道一共有多少种数目是包子大叔凑不出来的。

输入

第一行包含一个整数N。(1 <= N <= 100)
以下N行每行包含一个整数Ai。(1 <= Ai <= 100)

输出

一个整数代表答案。如果凑不出的数目有无限多个，输出INF。

例如，
输入：
2
4
5

程序应该输出：
6

再例如，
输入：
2
4
6

程序应该输出：
INF

样例解释：
对于样例1，凑不出的数目包括：1, 2, 3, 6, 7, 11。
对于样例2，所有奇数都凑不出来，所以有无限多个。

资源约定：
峰值内存消耗（含虚拟机） < 256M
CPU消耗 < 1000ms

请严格按要求输出，不要画蛇添足地打印类似：“请您输入…” 的多余内容。

注意：
main函数需要返回0;
只使用ANSI C/ANSI C++ 标准;
不要调用依赖于编译环境或操作系统的特殊函数。
所有依赖的函数必须明确地在源文件中 #include
不能通过工程设置而省略常用头文件。

提交程序时，注意选择所期望的语言类型和编译器类型。

思路

这一题和数学密切相关， 用到了欧几里得定理（最大公约数）

先看如果是两个数可以凑出几种

1.两个数a和b，如果互质的话，那么存在无限个数可以凑出，最大的不能凑出的数是a*b-a-b（记住就行）

2.如果a和b不互质的话，那么就存在无限个数凑不出。

多个数也是一样，判断是否互质，先判断两个，求出最大公约数，然后再用最大公约数和后面的数求最大公约数， 如果最后最大公约为1，那么说明互质，否则存在无限个数凑不出。

代码

#include <stdio.h>
int a[101];
int f[10001];
int gcd(int a, int b)//判断是否互质的函数，根据欧几里得定理
{
	if (b == 0) return a;
	return gcd(b, a%b);
}
int main()
{
	int n;
	scanf("%d", &n);
	int i, j;
	int g;
	f[0] = 1;
	for (i = 0; i < n; i++) {
		scanf("%d", &a[i]);
		if (i == 0) g = a[i];
		else g = gcd(a[i], g);
		for (j = 0; j <= 10000; j++) if (f[j]) f[j + a[i]] = 1;//找那些凑不出的数， 10000是因为a， b最大为100，所以100*100-200
	}
	if (g != 1) {//不互质
		printf("INF\n");
		return 0;
	}
	int ans = 0;
	for (j = 0; j <= 10000; j++) {
		if (!f[j]) {
			printf("%d ", j);
			ans++;
		}
	}
	printf("\n%d", ans);
	return 0;
}

思路2

用动态规划，和完全背包类似

dp【i】【j】，i 表示前 i 个数，j 表示如果dp【i】【j】有值，则 j 这个数能凑出

所以dp[i][j] = dp[i-1][j]+dp[i-1][j-a[i]]+dp[i-1][j-2a[i]]+dp[i-1][j-3a[i]]+…dp[i-1][j-(n*a[i])>=0]

代码

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