2017年第八届蓝桥杯c/c++B组 包子凑数 DP + 数学方法

题目提交点

点我进入官网提交该题

题目

思路

首先我们考虑，什么情况下要输出INF。
我们由裴蜀定理可以知道：$a\ast x+b\ast y=d$, 那么我们可以扩展到n个数之间： $x_1\ast a_1+x_2\ast a_2+x_3\ast a_3\dots \dots +x_n\ast a_n=d$,d为他们的最大公约数，那么我们就可以知道：能凑成数只有d的倍数，当$d=1$的时候能凑成的数是有限个（因为我们的$x_i$得是正数），当$d\ne 1$的时候，他只能凑出d的倍数，不能凑出来的数是无穷多个，此时输出INF。
然后我们考虑有限多个，因为包子是无限多个，所以我们就是完全背包的题目。那我们最多凑不出多少个包子呢？
做过 [蓝桥杯][2013年第四届真题]买不到的数目的朋友应该知道，求两个不互质的数凑不出来的数的数量格式应该是$(p-1)\ast (q-1)-1$。
当有n个数互质的时候凑不出来的数只会越来越来少，所以他最大不会超过$(100-1)\ast (100-1)-1$，索性我们直接设置为10000。

DP分析：

有三重循环，会超时，所以我们得要将我们的背包优化一下：
$f[i,j]=f[i-1,j]|f[i-1,j-a]|f[i-1,j-2\ast a]|\dots \dots |f[i-1,j-k\ast a]$
$f[i,j-a]=f[i-1,j-a]|f[i-1,j-2\ast a]|\dots \dots |f[i-1,j-k\ast a]$
所以可以得出：$f[i,j]=f[i-1,j]|f[i,j-a]$。

初始化：$f[0][0]=1$

AC代码

#include<bits/stdc++.h>
using namespace std;

#define _for(i, a, b) for (int i = (a); i < (b); ++i)
#define _rep(i, a, b) for (int i = (a); i <= (b); ++i)
#define For(i, a, b) for (int i = (a); i >= (b); --i)
#define debug(a) cout << #a << " = " << a << endl;
#define mod(x) (x) % MOD
#define ENDL "\n"
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<int> vi;

const int N = 10000 + 10, M = 100 + 10;
int a[M];
bool f[N];

int gcd(int a, int b) {
    return b ? gcd(b, a % b) : a;
}

int main()
{
#ifdef LOCAL
    freopen("data.in", "r", stdin);
#endif
    ios::sync_with_stdio(false); // 取消cin与stdin 的同步
    cout.tie(0), cin.tie(0);

    int n, d=0;
    cin >> n;
    _rep(i, 1, n) {
        cin >> a[i];
        d = gcd(d, a[i]);
    }

    if (d != 1) cout << "INF" << ENDL;
    else {
        f[0] = true;
        _rep(i, 1, n) _rep(j, a[i], N - 1) f[j] |= f[j - a[i]];

        int ans = 0;
        _rep(i, 1, N - 1) if (!f[i]) ans++;
        cout << ans << ENDL;
    }
    return 0;
}