总理同学的编程尝试

为大家见证传说中的初学者暴力编程(C/C++)

openjudge 百练 4148 生理周期 中国剩余定理

中国剩余定理水体,我个智障竟然调了一个多小时,原因竟然是…(我并不是个标题党)…求逆元求错了。。

【题面】

–> 题目链接 <–

4148:生理周期

总时间限制:1000ms 内存限制:65536kB


描述

人生来就有三个生理周期,分别为体力周期、感情周期和智力周期,它们的周期长度分别为23天、28天和33天。每一个周期中有一天是高峰。在高峰这天,人会在相应的方面表现出色。例如,在智力周期的高峰,人会思维敏捷,注意力容易高度集中。因为三个周期的长度不同,所以通常三个周期的高峰不会落在同一天。

对于每个人,想知道何时三个高峰落在同一天。对于每个周期,会给出从当前年份的第一天开始,到出现高峰的天数(不一定是第一次高峰出现的时间)。给定一个从当年第一天开始的天数,你的任务是输出从给定时间开始(不包括给定时间),下一次三个高峰落在同一天的时间(距给定时间的天数)。

例如:给定时间为10,下次出现三个高峰同一天的时间是12,则输出2(注意这里不是3)。

输入

输入包含多组数据,每一组数据由四个整数组成,数据以-1 -1 -1 -1 结束。

对于四个整数p, e, i和d,p, e, i分别表示体力、情感和智力高峰出现的时间(时间从当年的第一天开始计算)。d是给定的时间,可能小于p, e或i。所有给定时间是非负的并且小于或等于365,所求的时间小于或等于21252。

输出

从给定时间起,下一次三个高峰同一天的时间(距离给定时间的天数)。

样例输入

0 0 0 0
0 0 0 100
5 20 34 325
4 5 6 7
283 102 23 320
203 301 203 40
-1 -1 -1 -1

样例输出

Case 1: the next triple peak occurs in 21252 days.
Case 2: the next triple peak occurs in 21152 days.
Case 3: the next triple peak occurs in 19575 days.
Case 4: the next triple peak occurs in 16994 days.
Case 5: the next triple peak occurs in 8910 days.
Case 6: the next triple peak occurs in 10789 days.

【思路】

利用中国剩余定理,可以求出从模23余p,模28余e,模33余i的第一个非负整数。然后平移到大于d即可。

【代码】

#include<cstdio>
#include<cmath>
using namespace std;

#define LLint long long
void exgcd(LLint a, LLint b, LLint& x, LLint& y){
    if(b == 0){
        x = 1; y = 0;
    }else{
        LLint x0, y0;
        exgcd(b, a%b, x0, y0);
        x = y0;
        y = x0 - a/b * y0;
    }
}

LLint Inv(LLint a, LLint p){
    LLint x, y;
    exgcd(a, p, x, y);
    return (x % p + p) % p;
}

int main(){
    const int M = 23 * 28 * 33;
    const int M1 = M/23, M2 = M/28, M3 = M/33;
    int iM1 = Inv(M1, 23), iM2 = Inv(M2, 28), iM3 = Inv(M3, 33);
    int K1 = iM1*M1, K2 = iM2*M2, K3 = iM3*M3;
    int p, e, i, d, qcnt = 0;
    while(scanf("%d%d%d%d", &p, &e, &i, &d) && (p != -1)){
        p %= 23; e %= 28; i %= 33;
        int x = (p*K1 + e*K2 + i*K3) % M;
        int t = d/M;
        int n = t*M + x; while(n <= d) n += M;
        printf("Case %d: the next triple peak occurs in %d days.\n", ++qcnt, n-d);
    }
    return 0;
}

我之所以调试了这么长时间就是因为求逆元写错了。我最开始想试试欧拉定理的快速幂法求逆元,但是智障的我把28和33当成了素数。导致我把他们的phi函数当成了他们自身减一。。吸取经验教训,以后直接用拓展欧几里得求逆元。。

阅读更多
版权声明:文章纯属版主手敲,请同学们尊重版主的知识产权。 https://blog.csdn.net/GGN_2015/article/details/79144640
文章标签: 中国剩余定理
个人分类: 算法导论
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

不良信息举报

openjudge 百练 4148 生理周期 中国剩余定理

最多只允许输入30个字

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭