中国剩余定理（CRT）（同余）

首先给出模板：

	int CRT(int a[],int m[],int n)
	{
	    int M = 1;
	    int 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];
	        extend_Euclid(Mi, m[i], x, y);
	        ans = (ans + Mi * x * a[i]) % M;
	    }
	    if(ans < 0) ans += M;
	    return ans;
	}

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

设正整数m1,m2,mk两两互素，则同余方程组

有整数解。并且在模M=m1m2…*mk下的解是唯一的，用上面的代码求。

例题：
有些人认为一个人生命中有三个周期从他或她出生的那一天开始。这三个周期是身体，情绪和智力周期，它们的周期分别为23,28和33天。每个周期中有一个峰值。在一个周期的高峰期，一个人在相应的领域（身体，情感或精神）中表现最佳。由于三个循环具有不同的周期，因此三个循环的峰值通常在不同的时间发生。我们想确定任何人何时出现三峰（所有三个周期的峰值出现在同一天）。对于每个周期，您将获得从当前年度开始的一个峰值（不一定是第一个）发生的天数。您还将获得一个日期，表示为从当年年初开始的天数。您的任务是确定从给定日期到下一个三峰的天数。给定日期不计算在内。例如，如果给定日期为10且下一个三峰值出现在第12天，则答案为2，而不是3.如果在给定日期出现三峰值，则应给出下一次出现的天数三峰。

输入：
您将获得许多案例。每种情况的输入由一行四个整数p，e，i和d组成。值p，e和i分别是从当前年份开始的物理，情绪和智力周期达到峰值的天数。值d是给定日期，并且可以小于p，e或i中的任何一个。所有值均为非负数且最多为365，您可以假设在给定日期的21252天内将出现三峰值。输入的结束由一条线表示，其中p = e = i = d = -1。

由于28,23，33是互素的，可以用上面的中国剩余定理来求出，多少天后将会达到峰值。

	#include <iostream>
	#include <string.h>
	#include <stdio.h>
	 
	using namespace std;
	 
	int a[4], m[4];
	 
	void extend_Euclid(int a, int b, int &x, int &y)
	{
	    if(b == 0)
	    {
	        x = 1;
	        y = 0;
	        return;
	    }
	    extend_Euclid(b, a % b, x, y);
	    int tmp = x;
	    x = y;
	    y = tmp - (a / b) * y;
	}
	 
	int CRT(int a[],int m[],int n)
	{
	    int M = 1;
	    int 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];
	        extend_Euclid(Mi, m[i], x, y);
	        ans = (ans + Mi * x * a[i]) % M;
	    }
	    if(ans < 0) ans += M;
	    return ans;
	}
	 
	int main()
	{
	    int p, e, i, d, t = 1;
	    while(cin>>p>>e>>i>>d)
	    {
	        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);
	        if(ans <= d)
	            ans += 21252;
	        cout<<"Case "<<t++<<": the next triple peak occurs in "<<ans - d<<" days."<<endl;
	    }
	    return 0;
	}

上面的要求是两两互素，如果不是两个互素呢，下面的代码解决问题：

	#include <iostream>
	#include <string.h>
	#include <stdio.h>
	 
	using namespace std;
	typedef long long LL;
	const int N = 1005;
	 
	LL a[N], m[N];
	 
	LL gcd(LL a,LL b)
	{
	    return b? gcd(b, a % b) : a;
	}
	 
	void extend_Euclid(LL a, LL b, LL &x, LL &y)
	{
	    if(b == 0)
	    {
	        x = 1;
	        y = 0;
	        return;
	    }
	    extend_Euclid(b, a % b, x, y);
	    LL tmp = x;
	    x = y;
	    y = tmp - (a / b) * y;
	}
	 
	LL Inv(LL a, LL b)
	{
	    LL d = gcd(a, b);
	    if(d != 1) return -1;
	    LL x, y;
	    extend_Euclid(a, b, x, y);
	    return (x % b + b) % b;
	}
	 
	bool merge(LL a1, LL m1, LL a2, LL m2, LL &a3, LL &m3)
	{
	    LL d = gcd(m1, m2);
	    LL c = a2 - a1;
	    if(c % d) return false;
	    c = (c % m2 + m2) % m2;
	    m1 /= d;
	    m2 /= d;
	    c /= d;
	    c *= Inv(m1, m2);
	    c %= m2;
	    c *= m1 * d;
	    c += a1;
	    m3 = m1 * m2 * d;
	    a3 = (c % m3 + m3) % m3;
	    return true;
	}
	 
	LL CRT(LL a[], LL m[], int n)
	{
	    LL a1 = a[1];
	    LL m1 = m[1];
	    for(int i=2; i<=n; i++)
	    {
	        LL a2 = a[i];
	        LL m2 = m[i];
	        LL m3, a3;
	        if(!merge(a1, m1, a2, m2, a3, m3))
	            return -1;
	        a1 = a3;
	        m1 = m3;
	    }
	    return (a1 % m1 + m1) % m1;
	}
	 
	int main()
	{
	    int n;
	    while(scanf("%d",&n)!=EOF)
	    {
	        for(int i=1; i<=n; i++)
	            scanf("%I64d%I64d",&m[i], &a[i]);//这里的输入输出格式忽略。
	        LL ans = CRT(a, m, n);
	        printf("%I64d\n",ans);
	    }
	    return 0;
	}