东华大学 22上机题进阶-89 包子凑数(动态规划）

晚乔最美

于 2022-03-07 19:40:34 发布

阅读量163

点赞数

分类专栏：考研复试文章标签：算法 c++ 动态规划

本文链接：https://blog.csdn.net/s1050712899/article/details/123337700

版权

考研复试专栏收录该内容

6 篇文章 1 订阅

订阅专栏

本文介绍了一个有趣的编程题目，涉及到数学中的最大公约数和完全背包策略。小明发现包子铺的包子数量有特定规律，通过计算所有包子种类的最大公约数，判断能否凑出所有可能的数量。如果最大公约数不为1，则存在无限个无法凑出的数字；否则，通过动态规划求解哪些数量无法凑出，输出这些数量的个数。

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

                   89 包子凑数

小明几乎每天早晨都会在一家包子铺吃早餐。他发现这家包子铺有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。

思路：如果所有包子的最大公因子不为0，那么肯定有无穷个数是凑不出来的，如果最大公因子为0，那么问题就变成一个完全背包，只需要求容量1-1e4时候的最大价值dp[i]，如果dp[i]=i,说明能凑出来，否则就不能凑出来

ac代码：

#include<bits/stdc++.h>
#include<bitset>
#include<unordered_map>
#define pb push_back
#define bp __builtin_popcount
#define TIME cout << "RuningTime: " << clock() << "ms\n", 0
#define ls x<<1
#define rs x<<1|1
using namespace std;
typedef  long long ll;
const int inf=0x3f3f3f3f;
const int maxn=1e4+10;
const int MOD=9901;
const int mod = 1e9+7;
const double PI=3.14;
int lowbit(int x){return x&-x;}
ll gcd(ll x, ll y){ return y == 0 ? x : gcd(y, x%y); }
ll lcm(ll x, ll y){ return x / gcd(x, y)*y; }
inline ll dpow(ll a, ll b){ ll r = 1, t = a; while (b){ if (b & 1)r = (r*t) % MOD; b >>= 1; t = (t*t) % MOD; }return r; }
inline ll fpow(ll a, ll b,ll p){ ll r = 1, t = a; while (b){ if (b & 1)r = (r*t) % p; b >>= 1; t = (t*t) % p; }return r; }
priority_queue<int, vector<int>, greater<int> > q;
int dp[maxn],a[maxn];
int  main()
{
    int n;
    cin>>n;


    for(int  i=1;i<=n;i++)
    {
        cin>>a[i];
    }
    int x=__gcd(a[1],a[2]);

    for(int i=3;i<=n;i++)
        x=__gcd(x,a[i]);


    if(x!=1)
    {
        cout<<"INF"<<endl;

        return 0;
    }

    for(int i=1;i<=n;i++)
    {
        for(int j=a[i];j<=maxn-1;j++)
        {
            dp[j]=max(dp[j],dp[j-a[i]]+a[i]);
        }
    }

    int ans=0;
    for(int i=1;i<maxn;i++)
    {
        if(dp[i]!=i)ans++;
    }

    cout<<ans<<endl;






    return 0;
}