POJ 1006 Bioryhthms

中国剩余定理的应用

 

23，28，33两两互素

x=p(mod 23)

x=e(mod 28)

x=i(mod 33)  (两条线的等号代替三条线的等号。。。）

 

最后x-d再对21252取模，注意x-d可能为负数或0

 

由于23 28 33已经固定，可以先做预处理，算出Mi和Mi^-1 ，M*M^-1=1(mod ni)

代码：

#include<iostream> using namespace std; struct node { int d; int x; int y; }; int m[4],mm[4],a[4];//a[i]x=m[i](mod n[i]) int n[4]={23,28,33,0}; struct node extended_eculid(int a,int b) { struct node t1,t2; if(b==0) { t1.d=a; t1.x=1; t1.y=0; return t1; } t1=extended_eculid(b,a%b); t2.d=t1.d; t2.x=t1.y; t2.y=t1.x-(a/b)*t1.y; return t2; } int mles(int a,int b,int n) { struct node t; t=extended_eculid(a,n); return t.x*(b/t.d); } int main() { int i,ans,mul,ca=0; mul=23*28*33; while(scanf("%d%d%d%d",&a[0],&a[1],&a[2],&a[3])) { if(a[0]==-1&&a[1]==-1&&a[2]==-1&&a[3]==-1) break; for(i=0;i<3;i++) { a[i]%=n[i]; m[i]=mul/n[i]; } ans=0; for(i=0;i<3;i++) { int te; te=mles(m[i],1,n[i]); te=te+n[i]; ans+=te*a[i]*m[i]; } ans-=a[3]; ans=(ans+21252)%21252; //ans%=21252; if(ans==0) ans=21252; ca++; printf("Case %d: the next triple peak occurs in %d days./n",ca,ans); } return 0; }