中国剩余定理

中国剩余定理（CRT）的表述如下

设正整数为两两互素的数，那么同余方程组存在一组解：并且在模意义下有唯一解。解的计算公式：其中：，为模意义下的逆元。

代码：

void EXGCD(int a,int b,int &x,int &y){
    if(!b){ x=1; y=0; return; }
    EXGCD(b,a%b,x,y);
    int tmp=x; x=y;
    y=tmp-(a/b)*y;
}
int CRT(int a[],int m[],int n){
    int ans=0,M=1;
    for(int i=1;i<=n;i++) M*=m[i];
    for(int i=1;i<=n;i++){
        int x,y;
        int Mi=M/m[i];
        EXGCD(Mi,m[i],x,y);//x 即为Mi在模m[i]意义下的逆元 
        ans=(ans+Mi*x*a[i])%M;
    }
    while(ans<0) ans+=M;
    return ans;
}

题目： POJ1006
题意：人自出生起就有体力，情感和智力三个生理周期，分别为23，28和33天。一个周期内有一天为峰值，在这一
天，人在对应的方面（体力，情感或智力）表现最好。通常这三个周期的峰值不会是同一天。现在给出三个日
期，分别对应于体力，情感，智力出现峰值的日期。然后再给出一个起始日期，要求从这一天开始，算出最少
再过多少天后三个峰值同时出现。

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
int a[4],m[4]; 
void EXGCD(int a,int b,int &x,int &y){
    if(!b){ x=1; y=0; return; }
    EXGCD(b,a%b,x,y);
    int tmp=x; x=y;
    y=tmp-(a/b)*y;
}
int p,e,i,d;
int CRT(int a[],int m[],int n){
    int M=1,ans=0;
    for(int i=1;i<=n;i++) M*=m[i];
    for(int i=1;i<=n;i++){
        int x,y;
        int Mi=M/m[i];
        EXGCD(Mi,m[i],x,y);
        ans=(ans+x*Mi*a[i])%M;
    }
    while(ans<0) ans+=M;
    return ans;
}
int main(){
    int t=1;
    while(scanf("%d%d%d%d",&p,&e,&i,&d)==4){
        if(p == -1 && e == -1 && i == -1 && d == -1)  break;
        a[1]=p;a[2]=e;a[3]=i;
        m[1]=23;m[2]=28;m[3]=33;
        int ans=CRT(a,m,3);
        while(ans<=d) ans+=21252;

        printf("Case %d: the next triple peak occurs in %d days.\n",t++,ans-d);
    }
    return 0;
}