POJ 1006

2018/8/29

POJ 1006

#include <iostream>
#include <cstring>
#include <stdio.h>
#include <stdlib.h>
#include <cmath>
#include <map>
#include <algorithm>
using namespace std;

/*
中国剩余定理公式：(a+k*b)%b=c
X=p+x1*23=e+x2*28=i+x3*33
需要找到一个X=x1+x2+x3，如下：
(x1+x2+x3)%23=p x2和x3分别是23的倍数
(x1+x2+x3)%28=e x1和x3分别是28的倍数
(x1+x2+x3)%33=i x1和x2分别是33的倍数
综上：
x1是28和33的公倍数，且x1%23=p
x2是23和33的公倍数，且x2%28=e
x3是23和28的公倍数，且x3%33=i
由公式：a%b=c，那么（a*k）%b=a%b+a%b+…+a%b=c+c+…+c=kc（k>0）得：
x1%23=1,x2%28=1,x3%33=1，再分别乘以p、e、i即可
*/
int main()
{
    int k = 0;
    while(1){
        int p, e, i, d;
        scanf("%d%d%d%d", &p, &e, &i, &d);
        if(p == -1 && e == -1 && i == -1 && d == -1)
            break;
        //x(i)代表给定两数的公倍数
        int j, x1 = 28 * 33, x2 = 23 * 33, x3 = 23 * 28, r = 23 * 28 * 33;
        for(j = 1; ;j++)
            if(x1 * j % 23 == 1)
                break;
        x1 = x1 * j * p;
        for(j = 1; ;j++)
            if(x2 * j % 28 == 1)
                break;
        x2 = x2 * j * e;
        for(j = 1; ;j++)
            if(x3 * j % 33 == 1)
                break;
        x3 = x3 * j * i;
        //x为三个峰值同时出现的那一天
        int x = x1 + x2 + x3;
        //n=0表示给出三个峰值同时出现的时间，那么下一次同时出现的天数间隔为r
        int n = (x - d + r) % r;
        if(n == 0)
            n = r;
        printf("Case %d: the next triple peak occurs in %d days.\n", ++k, n);
    }

    return 0;
}