hdu 1573 X问题

利用模板求出最小解,如果最小解大于n或者无最小解(此模板返回-1),则无解

否则输出(n-ans)/dg+1(dg为a[i]的最小公倍数,ans为最小解),所有满足dg*i+ans的都符合要求,ans为最小的满足的数

#include<stdio.h>
using namespace std;
#define  ll long long

ll  A[11],B[11];//B[i]为余数 
ll dg,ans;//dg为A[i]的最小公倍数    ans 为最小解 
void exgcd(ll a, ll b, ll &d, ll&x, ll &y)
{
    if (!b) {d=a; x=1; y=0;}
    else
    {
        exgcd(b, a%b, d, y, x);
        y-=x*(a/b);
    }
}
ll gcd(ll a, ll b)
{
    if (!b) return a;
    else gcd(b, a%b);
}
ll china(ll n)
{
    ll a,b,d,x,y,dm;
    ll c,c1,c2;
    a=A[0]; c1=B[0];
    for (int i=1; i<n; i++)
    {
        b=A[i]; c2=B[i];
        exgcd(a, b, d, x, y);
        dm=b/d;
        c=c2-c1;
        if (c%d) return -1;
        x=((x*c/d)%dm+dm)%dm;//x可能为负
        c1=a*x+c1;
        a=a*b/d;
    }
    dg=a;//dg是最大公约数
    if (!c1)//考虑c1为0的情况
    {
        c1=1;
        for (int i=0; i<n; i++)
        {
            c1=c1*A[i]/gcd(c1, A[i]);
        }
        dg=c1;//此时dg为最小公倍数
    }
    return c1;//c1为最小的X
}
int main(){
	int t;
	scanf("%d",&t);
	while(t--){
		int n,m;
		scanf("%d%d",&n,&m);
		
		for(int i=0;i<m;i++)
		scanf("%lld",&A[i]);
		for(int i=0;i<m;i++)
		scanf("%lld",&B[i]);
		ans=china(m);
		if(ans==-1||ans>n)
		printf("0\n");
		else printf("%d\n",(n-ans)/dg+1);
	}
	return 0;
}

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值