[POJ1006]Biorhythms（中国剩余定理）

题目描述

题解

代码

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;

int Case,N,M,d,a,b,x,y,n[5],m[5],c[5];

void exgcd(int a,int b,int &x,int &y)
{
if (!b) x=1,y=0;
else exgcd(b,a%b,y,x),y-=a/b*x;
}
int main()
{
m[1]=23;m[2]=28;m[3]=33;M=m[1]*m[2]*m[3];
while (~scanf("%d%d%d%d",&c[1],&c[2],&c[3],&d))
{
if (c[1]==-1&&c[2]==-1&&c[3]==-1&&d==-1) return 0;
for (int i=1;i<=3;++i)
{
a=M/m[i];b=m[i];
exgcd(a,b,x,y);
x=((x%b)+b)%b;
if (!x) x+=b;
n[i]=a*x;
}
N=0;
for (int i=1;i<=3;++i)
N+=n[i]*c[i];
N-=d;
N=(N%M+M)%M;
if (!N) N+=M;
printf("Case %d: the next triple peak occurs in %d days.\n",++Case,N);
}
}