包子凑数(完全背包,凑数)

包子凑数

题目链接
小明几乎每天早晨都会在一家包子铺吃早餐。
他发现这家包子铺有 N 种蒸笼，其中第 i 种蒸笼恰好能放 Ai 个包子。
每种蒸笼都有非常多笼，可以认为是无限笼。
每当有顾客想买 X 个包子，卖包子的大叔就会迅速选出若干笼包子来，使得这若干笼中恰好一共有 X 个包子。
比如一共有 3 种蒸笼，分别能放 3、4 和 5 个包子。
当顾客想买 11 个包子时，大叔就会选 2 笼 3 个的再加 1 笼 5 个的（也可能选出 1 笼 3 个的再加 2 笼 4 个的）。
当然有时包子大叔无论如何也凑不出顾客想买的数量。
比如一共有 3 种蒸笼，分别能放 4、5 和 6 个包子。
而顾客想买 7 个包子时，大叔就凑不出来了。
小明想知道一共有多少种数目是包子大叔凑不出来的。

输入格式

第一行包含一个整数 N。
接下来 N 行，每行包含一个整数 Ai。

输出格式

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

数据范围

1≤N≤100,
1≤Ai≤100

输入样例1：

2
4
5

输出样例1：

6

输入样例2：

2
4
6

输出样例2：

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

算法分析

首先两个数凑不出来的最大数问题
引理：给定a，b，若d=gcd(a,b)>1,则一定不能凑出最大数
结论：如果 a,b 均是正整数且互质，那么由 ax+by,x≥0,y≥0不能凑出的最大数是 (a−1)(b−1)−1==a*b-a-b。
这个结论也是可以用到n个数之间的
之后就是可以无限选看能不能凑整=>完全背包的问题
我们直接用bool 来记录当前数是否能被凑出来,true就是能被凑出来
状态转移方程跟完全背包一样的,只不过不是数量而已
状态转移方程,dp[i][j]用到前i个组能凑出来j的数.

	dp[i][j]=dp[i-1][j];
    if(j>=a[i])
    {
         dp[i][j]|=dp[i][j-a[i]];
    }

由于dp只有两种状态,i层和i-1层而且,i层是继承i-1层的状态的
所以可以直接用一维来进行优化

代码实现

二维

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
using namespace std;
int n;
const int maxn=1e3+5;
bool dp[maxn][maxn];
int a[maxn];
int gcd(int a,int b)
{
    return b?gcd(b,a%b):a;
}
int main()
{
    cin>>n;
    int d;
    for(int i=1;i<=n;i++)
    {
        cin>>a[i];
        if(i==1)
            d=a[i];
        if(i>1)
            d=gcd(d,a[i]);
    }
    if(d!=1)    cout<<"INF"<<endl;
    else
    {
        dp[0][0]=true;
        for(int i=1;i<=n;i++)
        {
            for(int j=0;j<=maxn;j++)
            {
                dp[i][j]=dp[i-1][j];
                if(j>=a[i])
                {
                    dp[i][j]|=dp[i][j-a[i]];
                }
            }
        }
        int res=0;
        for(int i=1;i<=maxn;i++)
            if(!dp[n][i])
                res++;
        cout<<res<<endl;
    }
    return 0;
}

一维优化

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
using namespace std;
int n;
const int maxn=1e4+5;
bool dp[maxn];
int a[maxn];
int gcd(int a,int b)
{
    return b?gcd(b,a%b):a;
}
int main()
{
    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
    {
        dp[0]=true;
        for(int i=1;i<=n;i++)
        {
            for(int j=a[i];j<=maxn;j++)
            {
                    dp[j]|=dp[j-a[i]];
            }
        }
        int res=0;
        for(int i=1;i<=maxn;i++)
            if(!dp[i])
                res++;
        cout<<res<<endl;
    }
    return 0;
}