POJ 1006 Biorhythms（中国剩余定理）

题目传送门

http://poj.org/problem?id=1006

---

思路

中国剩余定理果题。

求$x\equiv a_i(\mod m_i)$的最小正整数解。

简要说一下步骤，证明什么的网上有很多，不再赘述。

首先令$M=\prod m_i$，令$M_i=\frac{M}{m_i}$，求出$M_i$的某个倍数使得其满足模$m_i$等于$a_i$，即$kM_i \equiv a_i(\mod m_i)$，发现$k\equiv a_iM_i^{-1}(\mod m_i)$，于是$kM_i\equiv a_iM_iM_i^{-1}(\mod m_i)$，其中$M_i^{-1}$是$M_i$关于$m_i$的逆元。（正确性其实是显然的）

将所有求得的数相加得到一解，尽可能减去所有$m_i$的$lcm$即得最小解。

答案就是一条公式：

$$ans=\sum a_iM_iM_i^{-1} \mod M$$

注意答案要大于0。（话说本题暴力也可以）

---

代码

#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>

using namespace std;

int d, Case, ans, MOD = 21252;
int m[5] = {0, 23, 28, 33}, M[5] = {0, 21252/23, 21252/28, 21252/33};
int A[5], InvM[5];

void Exgcd(int a, int b, int& x, int& y){
    if(!b){
        x = 1;
        y = 0;
        return;
    }
    Exgcd(b, a % b, y, x);
    y -= x * (a / b);
}

int CRT(){
    int x, y;
    for(int i = 1; i <= 3; i++){
        Exgcd(M[i], m[i], x, y);
        InvM[i] = (x % m[i] + m[i]) % m[i];
    }
    int res = 0;
    for(int i = 1; i <= 3; i++)
        res = (res + A[i] * M[i] % MOD * InvM[i] % MOD) % MOD;
    return res;
}

int main(){

    for(;;){
        for(int i = 1; i <= 3; i++)  scanf("%d", &A[i]);
        if(A[1] == -1)  break;
        scanf("%d", &d);
        ans = CRT() - d;
        if(ans <= 0)  ans += MOD;
        printf("Case %d: the next triple peak occurs in %d days.\n", ++Case, ans);
    }

    return 0;
} 

---

接下来要写三模数NTT。