包子凑数(如果gcd(a,b)=1,由于1能整除所有整数,此时ax+by能得到所有整数)

包子凑数

题目分两步, (1)判断结果是否为INF;(2)如果不是INF,统计数量。考点是“数论gcd+简单DP”。

    1. 什么时候答案不是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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值