中国剩余定理

（中国剩余定理CRT）

设m1,m2,...,mk是两两互素的正整数，即gcd(mi, mj) =1, i≠j, i,j = 1,2,...,k

则同余方程组：

　　x≡b1 mod m1

　　x≡b2 mod m

　　...

　　x≡bk mod mk

模[m1,m2,...,mk]有唯一解，即在[m1,m2,...,mk]的意义下，存在唯一的x,满足：x≡bi mod [m1,m2,...,mk], i = 1,2,...,k

 

---

 

实例：

　　有一个年级的同学，每9人一排多5人，每7人一排多1人，每5人一排多2人，问这个年级至少有多少人 ？（幸福123老师问的题目）题中9、7、5三个数两两互质。则〔7，5〕=35；〔9，5〕=45；〔9，7〕=63；〔9，7，5〕=315。为了使 35被9除余1，用35×8=280；使45被7除余1，用45×5=225；使63被5除余1，用63×2=126。然 后，280×5＋225×1＋126×2=1877，因为，1877>315，所以，1877－315×5=302，就是所求的数

 

算法实例：

题目见：ZOJ Problem Set - 1160  Biorhythms

 

#include <iostream> using namespace std; int main() { int i = 2, c1 = 0, c2 = 0, c3 = 0; int p1 = 23 * 28; int p2 = 23 * 33; int p3 = 28 * 33; int t = 23 * 28 * 33; while( i++ ){ if( (p1 * i) % 33 == 1 && c1 == 0 ) c1 = p1 * i; if( (p2 * i) % 28 == 1 && c2 == 0 ) c2 = p2 * i; if( (p3 * i) % 23 == 1 && c3 == 0 ) c3 = p3 * i; if( c1 && c2 && c3 ) break; } int N; cin >> N; while( N-- ){ int p, e, i, d, c = 1, s; while( cin >> p >> e >> i >> d, d != -1 ){ s = ( c1 * i + c2 * e + c3 * p ) % t - d; if( s <= 0 ) s += t; cout << "Case " << c++ << ": the next triple peak occurs in " << s << " days." << endl; } if( N ) cout << endl; } return 0; }