poj1006 生理周期

题目：
Description

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

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

当p = e = i = d = -1时，输入数据结束。
Output

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

采用以下格式：
Case 1: the next triple peak occurs in 1234 days.

注意：即使结果是1天，也使用复数形式“days”。
Sample Input

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
Sample Output

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.

代码如下：

#include<iostream>
#include<cstdio>
using namespace std;
#define MOD 21252
int ans;
int exgcd(int aa,int bb,int &x,int &y)
{
    if(!bb){
        x = 1;
        y = 0;
        return aa;
    }
    ans = exgcd(bb,aa % bb,x,y);
    int temp = x;
    x = y;
    y = temp - (aa / bb) * y;
    return ans;
}
int CRT(int a[],int m[])
{
    int res = 0,x,y;
    for(int i = 0;i < 3;i++){
        int aa = MOD / m[i];
        exgcd(aa,m[i],x,y);
        res += aa * x * a[i];
    }
    return res;
}
int main()
{
    int num = 1;
    int p,e,i,d,a[4],m[4] = {23,28,33};
    while(~scanf("%d%d%d%d",&p,&e,&i,&d)){
        if(p == -1 && e == -1 && i == -1 && d == -1) break;
        a[0] = p;
        a[1] = e;
        a[2] = i;
        int sum = CRT(a,m) - d;
        sum %= MOD;
        if(sum <= 0) sum += MOD;
        printf ("Case %d: the next triple peak occurs in %d days.\n", num++,sum);
    }
    return 0;
}

题意：
中文题，就不描述了。
思路：
刚刚学了中国剩余式定理，先来讲一下什么是中国式剩余定理：

证明我过几天在别的博客会补上的(提前给自己挖坑…）
回到这道题目上，这里把周期长度23，28，33当成mi，输入的三个数当作ai。因为要对结果取m1 * m2 * m3，所以k*(m1 * m2 * m3)这一部分就为0，只需要计算后面一部分就可以了。算出来的结果减去起始时间然后对232833取模，如果为负数还要加上23 * 28 * 33。因为题目对输出有范围要求。

最后把中国剩余定理的模板给大家看一下吧！(这是b站里一位大佬写的ppt)

ll Sunzi(ll *m,ll *a,int len){
    ll lcm = 1;
    ll ans = 0;
    for (int i=0;i<len;i++){
        ll k0,ki;
        ll d = ext_gcd(lcm,m[i],k0,ki);
        if ((a[i]-ans)%d!=0) return -1;
        else {
            ll t = m[i]/d;
            k0 = ( k0*(a[i]-ans)/d%t + t)%t;
            ans = k0*lcm + ans;
            lcm = lcm/d*m[i];
        }
    }
    return ans;
}