mi不互质时解同余方程组
存一发模板吧QAQ
讲解看看:这里
那个为什么要mod (m2/d)还是想了很久。。
大概就是最后要在mod(m2*m1/d)意义下,然后带入到
x=m1∗t1+a1
中乘了
m1
,所以前面
mod(m2∗m1/d)
就可以满足条件。
【代码】
#include <cstdio>
#include <iostream>
#include <algorithm>
#define N 2335
#define INF 0x7fffffff
using namespace std;
typedef long long ll;
ll read()
{
ll x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
int n;
ll M,R;
ll m[N],r[N];
void Exgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
if(b) {
Exgcd(b,a%b,d,y,x);
y-=(a/b)*x;
}
else x=1,y=0,d=a;
}
ll Ex_Crt()
{
M=m[1],R=r[1];ll x,y,d;
for(int i=2;i<=n;i++)
{
Exgcd(M,m[i],d,x,y);
if((r[i]-R)%d) return -1;
R=(r[i]-R)/d*x%(m[i]/d)*M+R;
M*=m[i]/d;
R%=M;
}
return (R+M)%M;
}
int main()
{
while(scanf("%d",&n)!=EOF)
{
for(int i=1;i<=n;i++) m[i]=read(),r[i]=read();
printf("%lld\n",Ex_Crt());
}
ll x,y,d;
Exgcd(24,6,d,x,y);
printf("%lld %lld\n",x,y);
return 0;
}