一个小问题(中国剩余定理非互质版)

时间限制:C/C++ 1秒,其他语言2秒
空间限制:C/C++ 32768K,其他语言65536K
64bit IO Format: %lld

题目描述

uu遇到了一个小问题,可是他不想答。你能替他解决这个问题吗?
问题:给你k对a和r是否存在一个正整数x使每队a和r都满足:x mod a=r,求最小正解x或无解。


输入描述:

第一行是正整数k(k<=100000)
接下来k行,每行有俩个正整数a,r(100000>a>r>=0)

输出描述:

在每个测试用例输出非负整数m,占一行。
如果有多个可能的值,输出最小的值。
如果没有可能的值,则输出-1。

示例1

输入

2
8 7
11 9

输出

31

练习了一下非互质版本的中国剩余定理~~~附证明图片


#include<stdio.h>
#include <algorithm>
#include<iostream>
#include<string.h>
#include<vector>
#include<stdlib.h>
#include<math.h>
#include<queue>
#include<deque>
#include<ctype.h>
#include<map>
#include<set>
#include<stack>
#include<string>
#include<algorithm>
#define ll long long
#define INF 0x3f3f3f3f
using namespace std;

ll gcd(ll a,ll b)
{
    return b?gcd(b,a%b):a;
}

ll exgcd(ll a,ll b,ll &x,ll &y)
{
    if(!b)
    {
        x=1;
        y=0;
        return a;
    }
    ll r=exgcd(b,a%b,y,x);
    y-=a/b*x;
    return r;
}

ll inv(ll a,ll n)
{
    ll x,y;
    ll t=exgcd(a,n,x,y);
    if(t!=1)
        return -1;

    return (x%n+n)%n;
}

ll china(ll *m,ll *a,int n)
{
    ll M=1,ans=0;

    for(int i=0;i<n;i++) M*=m[i];

    for(int i=0;i<n;i++)
    {
        ll x,y,mi=M/m[i];
        exgcd(mi,m[i],x,y);
        ans+=(mi*x*a[i])%M;
    }
    return (ans+M)%M;
}

int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        ll a1,n1,a2,n2;
        scanf("%lld %lld",&n1,&a1);

        int flag=0;
        for(int i=1;i<n;i++)
        {
            scanf("%lld %lld",&n2,&a2);

            ll x,y;
            ll c=a2-a1;
            ll d=exgcd(n1,n2,x,y);

            if(c%d)
                flag=1;

            x=(x*(c/d)%(n2/d)+(n2/d))%(n2/d);

            a1=a1+n1*x;
            n1=n1*n2/d;
        }

        if(flag)
            printf("-1\n");
        else
            printf("%lld\n",a1);
    }
    return 0;
}

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值