poj 1006 Biorhythms [osily]

这个题以前做过，但用的是枚举法， 想起当时做时有人说可用中国剩余定理，重看此题，一起研究后终于明白了。

#include <iostream> #include <cstring> #include <cmath> #include <cstdlib> #include <algorithm> #include <string> #include <cstdio> #include <climits> using namespace std; int b[3], w[3] = {23, 28, 33}; int ext_gcd(int a, int b, int &x, int &y) { int d; if (b == 0) { x = 1; y = 0; return a; } else { d = ext_gcd(b, a % b, y, x); y -= (a / b) * x; } return d; } int ch_re(int len) { int d, x, y, m, n, ret; ret = 0; n = 1; for (int i=0; i < len; i++) n *= w[i]; for(int i=0; i < len; i++) { m = n / w[i]; d = ext_gcd(w[i], m, x, y); ret = (ret + y*m*b[i]) % n; } return (n + ret%n) % n; } int main() { int d1, d2, d3, d; b[0] = 1, b[1] = 0, b[2] = 0; int t1 = ch_re(3); b[0] = 0, b[1] = 1, b[2] = 0; int t2 = ch_re(3); b[0] = 0, b[1] = 0, b[2] = 1; int t3 = ch_re(3); b[0]=b[1]=b[2]=0; int t = w[0]*w[1]*w[2]; int n = 1; while (scanf("%d%d%d%d", &d1, &d2, &d3, &d), d1!=-1) { int result = (t1*d1+t2*d2+t3*d3 + t - d) % t; if (result == 0) result = t; printf("Case %d: the next triple peak occurs in %d days./n", n++, result); } return 0; }

这是我写的完整代码，虽然可以把那几个数在本地算出来直接用，但那样的代码很让人费解。

最重要的一句是：

int result = (t1*d1+t2*d2+t3*d3 + t - d) % t;

 

这是让 t1*d1+t2*d2+t3*d3 和 w[i] 取模分别是d1, d2, d3。

看来数论的题目重在灵活运用定理，而不是证明定理。