POJ1006 Biorhythms

最新推荐文章于 2019-08-07 08:00:00 发布

forezxl

最新推荐文章于 2019-08-07 08:00:00 发布

阅读量583

点赞数 1

分类专栏： POJ 数论---(EX)CRT 蒟蒻zxl的Blog专栏文章标签： poj 中国剩余定理

本文链接：https://blog.csdn.net/a1799342217/article/details/78367984

版权

蒟蒻zxl的Blog专栏同时被 3 个专栏收录

460 篇文章 5 订阅

订阅专栏

POJ

22 篇文章 0 订阅

订阅专栏

数论---(EX)CRT

4 篇文章 0 订阅

订阅专栏

中国剩余定理

题目传送门

题意可以转化为求解方程组：

$a_1x\equiv 1(mod\ m_1)$
$a_2x\equiv 1(mod\ m_2)$
$a_3x\equiv 1(mod\ m_3)$
其中 $m_1=23$ ， $m_2=28$ ， $m_3=33$ 。
因为 $m1,m2,m3$ 两两互质，所以 $a_1=m_2*m_3$ ， $a_2=m_1*m_3$ ， $a_3=m_1*m_2$ 。
先用EXGCD算出解x，代入三个式子得到 $n1,n2,n3$ ，然后用中国剩余定理求解即可。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const LL m1=23,m2=28,m3=33,m=21252;
LL t=0,a1,a2,a3,n1,n2,n3,x,d;
void exgcd(LL a,LL &x,LL b,LL &y){
    if (b==0){
        x=1; y=0; return;
    }
    exgcd(b,y,a%b,x);
    y-=a/b*x;
}
int main(){
    while ((~scanf("%lld%lld%lld%lld",&a1,&a2,&a3,&d)&&(~a1||~a2||~a3||~d))){
        n1=m2*m3; n2=m1*m3; n3=m1*m2;
        LL x,y; 
        exgcd(n1,x,m1,y); n1=((n1*x)%m+m)%m;
        exgcd(n2,x,m2,y); n2=((n2*x)%m+m)%m;
        exgcd(n3,x,m3,y); n3=((n3*x)%m+m)%m;
        x=a1*n1+a2*n2+a3*n3-d;//中国剩余定理
        x%=m;if (x<=0) x+=m; 
        printf("Case %lld: the next triple peak occurs in %lld days.\n",++t,x);
    }
    return 0;
}