Poj 1006 中国剩余定理

来源:互联网 发布:赚流量的软件 编辑:程序博客网 时间:2024/05/22 10:49

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);    }}

0 0