poj 1006中国剩余定理模板

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

设正整数两两互素，则同余方程组

有整数解。并且在模下的解是唯一的，解为

其中，而为模的逆元。

模板：

int crt(int a[],int m[],int n)
{
    int e=;
    int mod=;
    for(int i=; i<=n; i++) mod*=m[i];
    for(int i=; i<=n; i++)
    {
        int temp=mod/m[i];
        e=(e+a[i]*temp*inv(temp,m[i]))%mod; // 这里求逆元 是对mi求，，
    }
    return e;
}

poj 1006ac代码：

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
int mod;
int exgcd(int a,int b,int &x,int &y)
{
    if(b==)
    {
       x=;
       y=;
       return a;
    }
    int temp=exgcd(b,a%b,y,x);
    y-=(a/b)*x;
    return temp;
}
int inv(int a,int b)
{
    int x,y;
    int temp=exgcd(a,b,x,y);
    if(x<) x+=b/temp;
    return x;
}
int crt(int a[],int m[],int n)
{
    int e=;
    int mod=;
    for(int i=; i<=n; i++) mod*=m[i];
    for(int i=; i<=n; i++)
    {
        int temp=mod/m[i];
        e=(e+a[i]*temp*inv(temp,m[i]))%mod; // 这里求逆元 是对mi求，，
    }
    return e;
}
int main()
{
    int a[],m[];
    m[]=;
    m[]=;
    m[]=;
    int s;
    int Case=;
    while(cin>>a[]>>a[]>>a[]>>s)
    {
        if(a[]==- && a[]==- && a[]==- && s==-) break;
        int e=;
        Case++;
        e=crt(a,m,);
        if(e <= s) e+=;
        printf("Case %d: the next triple peak occurs in %d days.\n",Case,e-s);
    }
    return ;
}