题
小明几乎每天早晨都会在一家包子铺吃早餐。这家包子铺有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
提示
对于样例,凑不出的数目包括:1, 2, 3, 6, 7, 11。
思路:
欧几里得定理 + 完全背包
【欧几里得定理】
对于不完全为 0 的整数 a,b,gcd(a,b)表示 a,b 的最大公约数。那么一定存在整数 x,y 使得 gcd(a,b)=ax+by。
如果包子的gcd(最大公约数)不为一,那么凑不出来的情况就有无数种。
如果为一,则可以利用完全背包去求解有多少种。
首先分析数据不是很大,最大可以构成的数字是100*100=10000。
所以我们可以设一个数组dp[10004],运用背包的思想,如果dp[j]可以,那么dp[j+a[i]]也是可以表示的。
最后再遍历整个dp数组,记录哪些不能表示即可。
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std ;
int gcd( int a , int b )
{
return b == 0 ? a : gcd(b , a % b ) ;
}
int main()
{
int n , a[105] , dp[10004] ;
cin >> n ;
for( int i = 0 ; i < n ; ++i )
{
cin >> a[i] ;
dp[a[i]] = 1 ;
}
int tem = gcd( a[0] , a[1] ) ;
for( int i = 2 ; i < n ; ++i )
{
tem = gcd( tem , a[i] ) ;
}
if( tem != 1 )
{
cout << "INF" << endl ;
return 0 ;
}
for( int j = 0 ; j < 10002 ; ++j )
{
if( dp[j] != 0 )
{
for( int i = 0 ; i < n ; ++i )
{
if( j + a[i] < 10003 )
dp[j+a[i]] = 1 ;
}
}
}
int num = 0 ;
dp[0] = 1 ;
for( int i = 0 ; i <= 10000 ; ++i )
{
if( dp[i] == 0 )
num++ ;
}
cout << num << endl ;
return 0 ;
}