生理周期（中国剩余定理）

【问题描述】

人生来就有三个生理周期，分别为体力、感情和智力周期，它们的周期长度为23 天、 28 天和33 天。每一个周期中有一天是高峰。在高峰这天，人会在相应的方面表现出色。例如，智力周期的高峰，人会思维敏捷，精力容易高度集中。因为三个周期的周长不同，所以通常三个周期的高峰不会落在同一天。对于每个人，我们想知道何时三个高峰落在同一天。对于每个周期，我们会给出从当前年份的第一天开始，到出现高峰的天数（不一定是第一次高峰出现的时间）。你的任务是给定一个从当年第一天开始数的天数，输出从给定时间开始（不包括给定时间）下一次三个高峰落在同一天的时间（距给定时间的天数）。例如：给定时间为10，下次出现三个高峰同天的时间是12，则输出2（注意这里不是3）。

【输入形式】

输入四个整数：p, e, i 和d。 p, e, i 分别表示体力、情感和智力高峰出现的时间（时间从当年的第一天开始计算）。d 是给定的时间，可能小于p, e, 或 i。 所有给定时间是非负的并且小于365, 所求的时间小于等于21252。

【输出形式】

从给定时间起，下一次三个高峰同天的时间（距离给定时间的天数）。

【样例输入】

0 0 0 0
0 0 0 100
5 20 34 325
4 5 6 7
283 102 23 320
203 301 203 40
-1 -1 -1 -1

【样例输出】

Case 1: the next triple peak occurs in 21252 days.
Case 2: the next triple peak occurs in 21152 days.
Case 3: the next triple peak occurs in 19575 days.
Case 4: the next triple peak occurs in 16994 days.
Case 5: the next triple peak occurs in 8910 days.
Case 6: the next triple peak occurs in 10789 days.

---

转化为数学问题：

读入p,e,i,d 4个整数，已知(n+d)%23=p，(n+d)%28=e， (n+d)%33=i ，求n。

由中国剩余定理，结合求乘法逆元，可以在O(1)内解出本题。

 

①M=23*28*33，M1=M/23，M2=M/28，M3=M/33，设Ni为在(mod Mi)下的乘法逆元

 

②求3个乘法逆元（此处可以手算），以求28*33(mod 23)的乘法逆元为例，即28*33 * x ≡ 1 (mod 23)

Ⅰ.因为28*33与23互质，可以使用拓展欧几里得算法，gcd(28*33, 23) = 1 = 28*33x + 23y，其中的x就是乘法逆元。

Ⅱ.可以使用欧拉定理：a^φ(m) ≡ 1 (mod m)，得出a的逆元为a^(φ(m)-1)，其中φ(m)为m的欧拉函数

则(28*33)^(φ(23)-1)即为乘法逆元，此处φ(23)=22。

 

③则n + d = (M1*N1*p + M2*N2*e + M3*N3*i + M) % M

---

#include<bits/stdc++.h>
using namespace std;
//------------------------------
//三个模两两互质，23,28,33
//由中国剩余定理
// |~  x = a1 (mod n1)
//-|   x = a2 (mod n2)
// |_  x = a3 (mod n3)
//记 N = n1*n2*……*nk
//  Ni = N/ni
//  xi 为(mod ni)下Ni的乘法逆元
//解:
//n1=23 N1=924 x1=6
//n2=28 N2=759 x2=19
//n3=33 N3=644 x3=2
//通解x=∑(ai*xi*Ni) (mod N)
//以下是代码实现
//------------------------------
long long n1=23,N1=924,x1=6;
long long n2=28,N2=759,x2=19;
long long n3=33,N3=644,x3=2;
long long N=n1*n2*n3;
int main()
{
    long long a1,a2,a3,a0,ans;
    int t=1;
    while(true)
    {
        cin>>a1>>a2>>a3>>a0;
        if(a1==-1)break;
        ans=(a1*x1*N1+a2*x2*N2+a3*x3*N3)%N;
        ans=(ans-a0)%N;
        while(ans<=0)ans+=N;
        //cout<<ans-a0<<'\n';
        printf("Case %d: the next triple peak occurs in %lld days.\n",t,ans);
        t++;
    }
    return 0;
}