中国剩余定理
题意可以转化为求解方程组:
a1x≡1(mod m1)
a2x≡1(mod m2)
a3x≡1(mod m3)
其中
m1=23
,
m2=28
,
m3=33
。
因为
m1,m2,m3
两两互质,所以
a1=m2∗m3
,
a2=m1∗m3
,
a3=m1∗m2
。
先用EXGCD算出解x,代入三个式子得到
n1,n2,n3
,然后用中国剩余定理求解即可。
代码:
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const LL m1=23,m2=28,m3=33,m=21252;
LL t=0,a1,a2,a3,n1,n2,n3,x,d;
void exgcd(LL a,LL &x,LL b,LL &y){
if (b==0){
x=1; y=0; return;
}
exgcd(b,y,a%b,x);
y-=a/b*x;
}
int main(){
while ((~scanf("%lld%lld%lld%lld",&a1,&a2,&a3,&d)&&(~a1||~a2||~a3||~d))){
n1=m2*m3; n2=m1*m3; n3=m1*m2;
LL x,y;
exgcd(n1,x,m1,y); n1=((n1*x)%m+m)%m;
exgcd(n2,x,m2,y); n2=((n2*x)%m+m)%m;
exgcd(n3,x,m3,y); n3=((n3*x)%m+m)%m;
x=a1*n1+a2*n2+a3*n3-d;//中国剩余定理
x%=m;if (x<=0) x+=m;
printf("Case %lld: the next triple peak occurs in %lld days.\n",++t,x);
}
return 0;
}