【问题描述】
人生来就有三个生理周期,分别为体力、感情和智力周期,它们的周期长度为23 天、 28 天和33 天。每一个周期中有一天是高峰。在高峰这天,人会在相应的方面表现出色。例如,智力周期的高峰,人会思维敏捷,精力容易高度集中。因为三个周期的周长不同,所以通常三个周期的高峰不会落在同一天。对于每个人,我们想知道何时三个高峰落在同一天。对于每个周期,我们会给出从当前年份的第一天开始,到出现高峰的天数(不一定是第一次高峰出现的时间)。你的任务是给定一个从当年第一天开始数的天数,输出从给定时间开始(不包括给定时间)下一次三个高峰落在同一天的时间(距给定时间的天数)。例如:给定时间为10,下次出现三个高峰同天的时间是12,则输出2(注意这里不是3)。
【输入形式】
输入四个整数: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.
转化为数学问题:
读入p,e,i,d 4个整数,已知(n+d)%23=p,(n+d)%28=e, (n+d)%33=i ,求n。
由中国剩余定理,结合求乘法逆元,可以在O(1)内解出本题。
①M=23*28*33,M1=M/23,M2=M/28,M3=M/33,设Ni为在(mod Mi)下的乘法逆元
②求3个乘法逆元(此处可以手算),以求28*33(mod 23)的乘法逆元为例,即28*33 * x ≡ 1 (mod 23)
Ⅰ.因为28*33与23互质,可以使用拓展欧几里得算法,gcd(28*33, 23) = 1 = 28*33x + 23y,其中的x就是乘法逆元。
Ⅱ.可以使用欧拉定理:a^φ(m) ≡ 1 (mod m),得出a的逆元为a^(φ(m)-1),其中φ(m)为m的欧拉函数
则(28*33)^(φ(23)-1)即为乘法逆元,此处φ(23)=22。
③则n + d = (M1*N1*p + M2*N2*e + M3*N3*i + M) % M
#include<bits/stdc++.h>
using namespace std;
//------------------------------
//三个模两两互质,23,28,33
//由中国剩余定理
// |~ x = a1 (mod n1)
//-| x = a2 (mod n2)
// |_ x = a3 (mod n3)
//记 N = n1*n2*……*nk
// Ni = N/ni
// xi 为(mod ni)下Ni的乘法逆元
//解:
//n1=23 N1=924 x1=6
//n2=28 N2=759 x2=19
//n3=33 N3=644 x3=2
//通解x=∑(ai*xi*Ni) (mod N)
//以下是代码实现
//------------------------------
long long n1=23,N1=924,x1=6;
long long n2=28,N2=759,x2=19;
long long n3=33,N3=644,x3=2;
long long N=n1*n2*n3;
int main()
{
long long a1,a2,a3,a0,ans;
int t=1;
while(true)
{
cin>>a1>>a2>>a3>>a0;
if(a1==-1)break;
ans=(a1*x1*N1+a2*x2*N2+a3*x3*N3)%N;
ans=(ans-a0)%N;
while(ans<=0)ans+=N;
//cout<<ans-a0<<'\n';
printf("Case %d: the next triple peak occurs in %lld days.\n",t,ans);
t++;
}
return 0;
}