Poj 1006 中国剩余定理

Description

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

#include"cstdio"
#include"iostream"
#include"cmath"
using namespace std;
typedef int LL;
int e_gcd(int a,int b,int &x,int &y)    //求a*x+m*y=1的特解x
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    else
    {
        int res=e_gcd(b,a%b,x,y);
        LL temp;
        temp=x;
        x=y;
        y=temp-a/b*y;
        return res;
    }


}
LL cal(LL a,LL b,LL c)  //求乘法逆元
{
    LL s,k;
    LL ans=e_gcd(a,b,s,k);
    if(c%ans!=0)
        return -1;
    s*=c/ans;
    b/=ans;
    if(b<0)
        b=abs(b);
    return (s+b)%b;    //有特解得出最小解
}
int main()
{
    int beg,s,n1,n2,n3,t1=23,t2=28,t3=33;
    while(scanf("%d %d %d %d",&beg,&n1,&n2,&n3)!=EOF)
    {
        int y1=t2*t3,y2=t1*t3,y3=t1*t2;  //y1,y2,y3为衍数（局部公倍数）
        //printf("%d %d %d\n",y1,y2,y3);
        int c1=cal(y1,t1,1);   //ci为乘率（乘法逆元）
        int c2=cal(y2,t2,1);
        int c3=cal(y3,t3,1);
        s=(n1*y1*c1+n2*y2*c2+n3*y3*c3)%(t1*t2*t3);
        printf("%d\n",(s-beg)%21252);
    }
}