POJ.1006 Biorhythms (拓展欧几里得+中国剩余定理)

最新推荐文章于 2020-05-04 19:43:53 发布

原创最新推荐文章于 2020-05-04 19:43:53 发布 · 368 阅读

2 ·

CC 4.0 BY-SA版权

数论---拓展欧几里得同时被 2 个专栏收录

5 篇文章

订阅专栏

数论---中国剩余定理

2 篇文章

订阅专栏

本文介绍了解决POJ.1006Biorhythms问题的方法，通过使用拓展欧几里得算法与中国剩余定理来找到特定日期。代码实现了输入周期并计算下一个三重峰值出现的日期。

部署运行你感兴趣的模型镜像

POJ.1006 Biorhythms (拓展欧几里得+中国剩余定理)

题意分析

(1)

不妨设日期为x，根据题意可以列出日期上的方程：

(2)

化简可得：

(3)

根据中国剩余定理求解即可。

代码总览

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
typedef int ll;
ll p,e,i,d;
void exgcd(ll a, ll b, ll& d, ll& x, ll &y)
{
    if(!b){
        d = a,x = 1,y = 0;
    }else{
        exgcd(b, a % b, d, y, x);
        y -= x * (a / b);
    }
}
int CRT(int a[], int m[])
{
    int ans = 0;
    int M = 1;
    int x,y,gcd;
    for(int i = 1; i<=3;++i)M*=m[i];
    for(int i =1 ;i<=3;++i){
        int Mi = M/m[i];
        exgcd(Mi,m[i],gcd,x,y);
        ans = (ans + Mi * a[i] * x) % M;
    }
    if(ans<=0) ans+=M;
    if(ans<d) ans+=M;
    return ans;
}
int main()
{
    //freopen("in.txt","r",stdin);
    int kase = 0;
    while(scanf("%d %d %d %d",&p,&e,&i,&d)){
        if(p == -1 && e == -1 && i == -1 && d == -1) break;
        int a[] = {0,p,e,i};
        int m[] = {0,23,28,33};
        int ans = CRT(a,m);
        printf("Case %d: the next triple peak occurs in %d days.\n",++kase,ans-d);
    }
    return 0;
}