三个线性同余方程组的计算机解决方案（C程序）

三个线性同余方程组的计算机解决方案（C程序）

在我国的《孙子算经》中有着这样一个的题目：有物不知其数，三个一数余二，五个一数余三，七个一数又余二，问该物总数几何？在我国还流传着这样一个故事：相传汉高祖刘邦问大将军韩信统御兵士多少，韩信答说，每3人一列余1人、5人一列余2人、7人一列余4人、13人一列余6人……。刘邦茫然而不知其数。关于这类问题的解决方案，不同教育阶段都有所不同。小学阶段主要的解决方案可能就是寻找最大公倍数，联立等式。而在高等教育阶段，学习了《初等数论》的相关知识后，我们可以把这样的问题抽象成多个线性同余方程组来解决。

首先介绍一下线性同余方程以及其求解方法。

若某未知数ax与常数b对常数n取模后结果相同，可以表示成ax≡b (mod n)的方程。此方程有解当且仅当 b 能够被 a 与 n 的最大公约数整除(记作 gcd(a,n) | b)。这时，如果 x0 是方程的一个解，那么所有的解可以表示为:

{x0+kn/d|(k∈z)}

其中 d 是a 与 n 的最大公约数。在模 n 的完全剩余系 {0,1,…,n-1} 中，恰有 d 个解。

通过这样的方程，回归到文章的最开始，对于《孙子算经》中的问题，可以联立如下方程组：

那么，如何解决这个方程组？

依据中国余数定理，设有一数N，分别被两两互素的几个数a1、a2、……an相除得余数R1、R2、……Rn，即：

N≡Ri(mod ai),(i=1、2、……n),

只需求出一组数K，使满足

R1(mod ai),(i=1、2、……n),

即可求出适合该线性同余方程组的最小整数解（P是整数，M=a1×a2×……×an）。

看到这里相信大家对于线性同余方程组已经有了初步的理解。线性同余方程组本质上就是寻找多个关键数，这些关键数满足的性质一定是：若除数集中元素共有n个，则关键数集中任意一个元素都可以被除数集中n-1个元素整除。对于这个问题，我们可以找到3个关键数，分别是70，21，15。观察发现，70被3除余1而被5，7均可整除；21被5整除余1而被3，7均可整除；15被7除余1而被3，5均可整除。由此可得式：70a+21b+15c。这个式子满足题目，但它并不是最小解，接下来可以减去105N（N=1,2,3……）。

当然这种解决方案属于数理领域，对于信息化的今天，我们完全可以利用计算机来更加高效快捷的计算出这样的线性同余方程组的答案。下文将通过C程序来提供任意三个线性同余方程组的计算工具。

还是观察上文。一个数被3除余2，被5除余3，被7除余2。中国余数定理给出的解决方案计算庞杂且不易理解，丢番图方程的解决方案可能无法求出给定范围内的全部解。而使用计算机，通过遍历的思想，我们可以很轻松的求出给定范围内的全部解。

算法核心思路：设置if条件判断，当且仅当i对3取模为2且i对5取模为3且i对7取模余2，输出一次结果，否则不输出。调用C程序中的if语句以及取模函数，可以得到核心代码段：

if (i % 3 == 2 && i % 5 == 3 && i % 7 == 2)
			printf("该整数为%d\n", i);

核心问题已经解决了，接下来的问题就是：如何将该算法封装成一个一般性的函数接口，每次调用只需要输入上限值，除数与相对应的余数就可以直接得出结果，而不用每次都在函数体内进行修改？

这里我想到了利用数组来实现。定义两个数组D和R，分别用于存放除数集与余数集，通过for循环嵌套scanf语句的方案，依次对除数与相应余数进行输入，从而实现对函数的封装。

整体代码如下：

#include <stdio.h>
#include <stdlib.h>
int ChinaRemainderTheorem(int n) {
	int D[3];	//D数组用于存放除数
	int R[3];	//R数组用于存放余数
	for (int i = 0; i <= 2; i++) {
		D[i] = 0;
		R[i] = 0;
	}			//初始化两个数组
	printf("请依次输入三个除数（结束请按回车）\n");
	for (int i = 0; i <= 2; i++) {
		scanf_s("%d", &D[i]);
	}
	printf("请依次输入与三个除数对应的三个余数（结束请按回车）\n");
	for (int i = 0; i <= 2; i++) {
		scanf_s("%d", &R[i]);
	}
	for (int i = 0; i < n; i++)
		if (i % D[0] == R[0] && i % D[1] == R[1] && i % D[2] == R[2])
			printf("该整数为%d\n", i);
	printf("\n\n");
	return 0;
}
int main()
{
	int a;
	int k;
	k = 1;
	while (k==1) {
		printf("请输入上限。\n");
		scanf_s("%d", &a);
		ChinaRemainderTheorem(a);
	}
	return 0;
}

运行成果展示：