包子凑数
题目分两步, (1)判断结果是否为INF;(2)如果不是INF,统计数量。考点是“数论gcd+简单DP”。
-
- 什么时候答案不是INF?gcd(a,b)==1
什么时候答案不是INF?也就是说,除了少数一些整数无法组合得到,其他所有的整数都能得到。
首先看2个数a、b的情况,结论是:若gcd(a,b)=d,则答案不是INF。
以两个数4、5为例,任取x个4和y个5(x,y≥0),它们能组合得到的数是:
4x+5y=c
若把c看成常数,这是一个二元一次方程。
罗勇军老师在其“线性丢番图方程 --算法竞赛专题解析(21):数论”博文中给出了二元一次方程ax+by=c有整数解的 定理:设a,b是整数且gcd(a,b)=d,如果d不能整除c,那么方程ax+by=c没有整数解,如果d能整除c,那么存在无穷多个整数解。
显然,如果gcd(a,b)=1,由于1能整除所有整数,此时ax+by能得到所有整数。
在例子4x+5y=c中,因为gcd(4,5)=1,那么不管c是什么整数,都存在整数解x、y ,也就是说答案不是INF。
不过,x、y的解可能是负整数,而本题要求x、y是非负整数,所以仍然有不能表达的数,且最大的不能表达的数不会超过 x*y
- 求多个的数的公因子
int res=0;
for(int i=1;i<=n;i++){
cin>>a[i];
res=gcd(a[i],res);
}
因为观察gcd函数可知,两个数求公因子,一个数为0,直接返回另一个数作为公因子
或者
in>>a[1];
res=a[1];
for(int i=2;i<=n;i++){
cin>>a[i];
res=gcd(a[i],res);
}
- 求凑不出的数——正难则反,求能求出俩的数
- 对于任意两个满足gcd(x,y)=1的x,y,其不能表示出来的最大数为(x-1)*(y-1)-1,而100以内 互质的最大的两个数为99和98,因此不能表示的最大数不超过10000,额,不知道咋证明,总之对于
gcd(x,y)=1
,大于x*y
的数都一定可以被表示
dp[0][0]=true;
for(int i=1;i<=100;i++){
for(int j=0;j<=10000;j++){
dp[i][j]=dp[i-1][j];
if(j>=a[i])
dp[i][j]|=dp[i][j-a[i]];
}
}
dp[0]=true; //勿忘初始化
for(int i=1;i<=n;i++){
for(int j=0;j<=10000;j++){
// dp[i][j]=dp[i-1][j];滚动数组直接继承i-1时对应的j情况
if(j>=a[i])
dp[j]|=dp[j-a[i]];
}
}
dp[0]=true; //勿忘初始化
for(int i=1;i<=n;i++){
for(int j=0;j+a[i]<=10000;j++){
// dp[i][j]=dp[i-1][j];
if(dp[j])
dp[j+a[i]]=true;
}
}
#include <iostream>
using namespace std;
int n;
const int N=105;
int a[N];
bool dp[10010];//不要下意识就写[N][N]
int gcd(int a,int b){
if(!b)return a;
else return gcd(b,a%b);
}
int main()
{
cin>>n;
int res=0;
cin>>a[1];
res=a[1];
for(int i=2;i<=n;i++){
cin>>a[i];
res=gcd(a[i],res);
}
if(res!=1){
cout<<"INF";
}
else{
dp[0]=true; //勿忘初始化
for(int i=1;i<=n;i++){
for(int j=0;j<=10000;j++){
// dp[i][j]=dp[i-1][j];
if(j>=a[i])
dp[j]|=dp[j-a[i]];
}
}
int num=0;
for(int i=0;i<=10000;i++){
if(!dp[i])num++;
}
cout<<num;
}
return 0;
}