首先给出模板:
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;
}