POJ-1006 Biorhythms 中国剩余定理

POJ-1006 Biorhythms

原题链接🔗，点击进入

题目大意：

由于题目过于冗长，在此只描述大致题意：人出生就伴随：体力，情感和智力三个生理周期，他们的周期分别为23天，28天，33天，每个周期内有一天为该阶段的峰值，这一天人在对应的方面的表现的最好。一般情况下，这三天不会在同一天，现在分别给出三个日期，分别对应体力§，情感(e)，智力(i)出现峰值的日期，以及一个起始日期(d)，问：从起始日期开始，最少经过多少天三个方面同时到达峰值。

思路：

关于中国剩余定理的相关内容，可以参考知乎相关内容

设三重峰值的出现为K，则
(K+d) % p = 0,
(K+d) % e = 0,
(K+d) % i = 0,
因此我们只需要求出K+d，即可求出最后的答案。对于这种同余求值的问题，就需要考虑使用中国剩余定理(CRT)来进行解决。

推导细节

令m1 = p, m2 = e, m3 = i

对于本题，先构造同余方程组：

① K + d ≡ p (mod 23)
② K + d ≡ e (mod 28)
③ K + d ≡ i (mod 33)

令: m=m1m2m3=23 * 28 * 33 = 21252 (最小公倍数)
令 M1 = m / m1, M2 = m / m2, M3 =m / m3,
得：

M1 = 28 * 33 = 924 M2 = 23 * 33 = 759 M3 = 23 * 28 = 644

再求出M1,M2,M3的逆元t1,t2,t3： 即：

① M1 * t1 ≡ 1 (mod m1)
② M2 * t2≡ 1 (mod m2)
③ M3 * t3 ≡ 1 (mod m3)

计算大致如下：

M1对23求余为1: M1 * a % 23 = 1 得：t1 = 6, 6 * 924 = 5544;
M2对28求余为1: M2 * b % 28 = 1 得：t2 = 19, 759 * 28 = 14421
M3对33求余为1: M3 * b % 33 = 1 得：t3 = 2, 644 * 33 = 1288

所以， K + d = (M1 * t1 * m1 + M2 * t2 * m2 + M3 * t3 * m3) % m

代入得： K + d = (5544 * p *+14421 * e + 1288 * i) % 21252
变式得： K = (5544 * p *+14421 * e + 1288 * i - d) % 21252

因为题目需要的是最小解，因此:
K = (5544 * p *+14421 * e + 1288 * i - d + 21252) % 21252

AC代码

#include <iostream>
using namespace std;
int main() {
	int p,e,i,d;
	int cnt = 1;
	while(cin >> p >> e >> i >> d){
		if(p != -1 && e != -1 && i != -1 && d != -1){
			int ans = (5544 * p +14421 * e + 1288 * i - d + 21252) % 21252;
			if(ans == 0){
				cout << "Case " << cnt++<< ": the next triple peak occurs in 21252 days." << endl;
			}
			else{
				cout << "Case " << cnt++ << ": the next triple peak occurs in "<< ans << " days." << endl;
			}
		}
		else{
			break;
		}
	}
	return 0;
}