时间限制:C/C++ 1秒,其他语言2秒
空间限制:C/C++ 32768K,其他语言65536K
64bit IO Format: %lld
空间限制:C/C++ 32768K,其他语言65536K
64bit IO Format: %lld
题目描述
uu遇到了一个小问题,可是他不想答。你能替他解决这个问题吗?
问题:给你k对a和r是否存在一个正整数x使每队a和r都满足:x mod a=r,求最小正解x或无解。
问题:给你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;
}