POJ 1006 Biorhythms （中国剩余定理）

#include<stdio.h>

int exgcd(int a,int b,int &x,int &y)
{
	if(b==0)
	{
		x=1;
		y=0;
		return a;
	}
	int ans=exgcd(b,a%b,x,y);
	int temp=x;
	x=y;
	y=temp-a/b*y;
	return ans;
}
// ans==ai mod mi ....
int crt(int a[],int m[],int n)
{
	int i,M=1;
	int ans=0;
	for(i=1;i<=n;i++)
		M*=m[i];
	for(i=1;i<=n;i++)
	{
		int x,y;
		int Mi=M/m[i];
		int gcd=exgcd(Mi,m[i],x,y);
		ans=(ans+Mi*x*a[i])%M;
	}
	if(ans<0) ans+=M;
	return ans;
}
int main()
{
	int cas=1;
	int p,e,i,d;
	int a[10],m[10];
	while(scanf("%d %d %d %d",&p,&e,&i,&d)!=EOF)
	{
		if(p==-1 && e==-1 && i==-1 && d==-1)
			break;
		m[1]=23;
		m[2]=28;
		m[3]=33;
		a[1]=p;
		a[2]=e;
		a[3]=i;
		int ans=crt(a,m,3);
		ans=ans-d;
		while(ans<=0)
			ans+=21252;
		printf("Case %d: the next triple peak occurs in %d days.\n",cas++,ans);

	}
return 0;
}