在一千多年前的《孙子算经》中,有这样一道算术题:有物不知其数,三三数之剩二,五五数之剩三,七七数之剩二。问物几何?即,一个整数除以三余二,除以五余三,除以七余二,求这个整数。《孙子算经》中首次提到了同余方程组问题,以及以上具体问题的解法,因此在中文数学文献中也会将中国剩余定理称为孙子定理。
宋朝数学家秦九韶于1247年《数书九章》卷一、二《大衍类》对“物不知数”问题做出了完整系统的解答。明朝数学家程大位在《算法统宗》中将解法编成易于上口的《孙子歌诀》:三人同行七十希,五树梅花廿一支,七子团圆正半月,除百零五便得知。
这个歌诀给出了模数为3、5、7时候的同余方程的秦九韶解法。意思是:将除以3得到的余数乘以70,将除以5得到的余数乘以21,将除以7得到的余数乘以15,全部加起来后再减去105或者105的整数倍,得到的数就是答案(除以105得到的余数则为最小答案)。
因此按歌诀求出的结果就是23.
这个题目其实可以用中国剩余定理来解决。用现代数学的语言来说明的话,中国剩余定理给出了以下的一元线性同余方程组:
中国剩余定理说明:假设整数m1, m2, … , mn其中任两数互质,则对任意的整数:a1, a2, … , an,方程组(S)有解,并且通解可以用如下方式构造得到:
不妨设
则方程组(S)的通解为:
我们先来看一道杭电的题目HDU1370
题目大意是:有3个循环周期,周期天数分别为23、28、33。对于某一年,已知某年这3个周期的某一峰值分别是当年的第p、e、i天, 问从第d天开始到最近一个满足3个周期都达到峰值的日期还有多少天?
分析:
我们假设这里d = 0,那么题目就变得和上面孙子算经的题目一样了,即寻找一个数使得该数字除以23,28,33之后的余数相同。题目好像数据好像比较小,可以直接打表,先看看打表的方法,如下:
#include <iostream>
#include <cstdio>
using namespace std;
int a[50][50][50];
int main()
{
int t;
cin >> t;
int k = 1;
int p, e, i, d;
for (int i = 0 ; i < 21252 ; i++)
{
a[i%23][i%28][i%33] = i;
}
while (cin >> p >> e >> i >> d)
{
if (p == -1 && e == -1 && i == -1 && d == -1)
{
break;
}
int x = a[p % 23][e % 28][i % 33] - d;
printf("Case %d: the next triple peak occurs in %d days.\n", k++, x <= 0 ? x + 21252 : x);
}
return 0;
}
.