POJ 1006 Biorhythms .

裸的中国剩余定理

求三个数的周期就是求他们的最小倍数

而且这三个数正好互素，所以肯定是中国剩余定理

#include<iostream>
using namespace std;
int gcdEx(int a,int b,int& x,int& y)
{ //求ax+by=gcd(a,b) 的整数解 返回gcd(a,b);
if(b==0) {
x=1; y=0; return a;
}
int x1,y1;
int gcd=gcdEx(b,a%b,x1,y1);
x=y1;
y=x1-a/b*y1;
return gcd;
}
int n0;
int CRT(int* a,int* n,int k)
{ //a[i]表示第i个余数，n[i]表示第i个除数，一共有k个数，返回x值
n0=1;
for(int i=0;i<k;i++) n0*=n[i]; //n0=n1*n2....*nk;
int ans=0,x,y,m;
for(int i=0;i<k;i++)
{
m=n0/n[i];
gcdEx(m,n[i],x,y);
ans=(ans+a[i]*x*m)%n0;
}
if(ans<0) ans+=n0;
return ans;
}
int main()
{
int a[3],day,kase=0;
int n[3]={23,28,33};
while(cin>>a[0]>>a[1]>>a[2]>>day&&a[0]!=-1)
{
int ans=CRT(a,n,3)-day;
if(ans<=0) ans+=n0;
cout<<"Case "<<++kase<<": the next triple peak occurs in "<<ans<<" days."<<endl;
}
return 0;
}