openjudge 百练 4148 生理周期 中国剩余定理

中国剩余定理水体，我个智障竟然调了一个多小时，原因竟然是…（我并不是个标题党）…求逆元求错了。。

【题面】

–> 题目链接 <–

4148:生理周期

总时间限制:1000ms 内存限制:65536kB

---

描述

人生来就有三个生理周期，分别为体力周期、感情周期和智力周期，它们的周期长度分别为23天、28天和33天。每一个周期中有一天是高峰。在高峰这天，人会在相应的方面表现出色。例如，在智力周期的高峰，人会思维敏捷，注意力容易高度集中。因为三个周期的长度不同，所以通常三个周期的高峰不会落在同一天。

对于每个人，想知道何时三个高峰落在同一天。对于每个周期，会给出从当前年份的第一天开始，到出现高峰的天数（不一定是第一次高峰出现的时间）。给定一个从当年第一天开始的天数，你的任务是输出从给定时间开始（不包括给定时间），下一次三个高峰落在同一天的时间（距给定时间的天数）。

例如：给定时间为10，下次出现三个高峰同一天的时间是12，则输出2（注意这里不是3）。

输入

输入包含多组数据，每一组数据由四个整数组成，数据以-1 -1 -1 -1 结束。

对于四个整数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.

【思路】

利用中国剩余定理，可以求出从模23余p，模28余e，模33余i的第一个非负整数。然后平移到大于d即可。

【代码】

#include<cstdio>
#include<cmath>
using namespace std;

#define LLint long long
void exgcd(LLint a, LLint b, LLint& x, LLint& y){
    if(b == 0){
        x = 1; y = 0;
    }else{
        LLint x0, y0;
        exgcd(b, a%b, x0, y0);
        x = y0;
        y = x0 - a/b * y0;
    }
}

LLint Inv(LLint a, LLint p){
    LLint x, y;
    exgcd(a, p, x, y);
    return (x % p + p) % p;
}

int main(){
    const int M = 23 * 28 * 33;
    const int M1 = M/23, M2 = M/28, M3 = M/33;
    int iM1 = Inv(M1, 23), iM2 = Inv(M2, 28), iM3 = Inv(M3, 33);
    int K1 = iM1*M1, K2 = iM2*M2, K3 = iM3*M3;
    int p, e, i, d, qcnt = 0;
    while(scanf("%d%d%d%d", &p, &e, &i, &d) && (p != -1)){
        p %= 23; e %= 28; i %= 33;
        int x = (p*K1 + e*K2 + i*K3) % M;
        int t = d/M;
        int n = t*M + x; while(n <= d) n += M;
        printf("Case %d: the next triple peak occurs in %d days.\n", ++qcnt, n-d);
    }
    return 0;
}

我之所以调试了这么长时间就是因为求逆元写错了。我最开始想试试欧拉定理的快速幂法求逆元，但是智障的我把28和33当成了素数。导致我把他们的phi函数当成了他们自身减一。。吸取经验教训，以后直接用拓展欧几里得求逆元。。