网络安全-数论入门-孙子定理(CRT)

1. 简介

孙子定理是中国古代求解一次同余式组（见同余）的方法。是数论中一个重要定理。又称中国余数定理。一元线性同余方程组问题最早可见于中国南北朝时期（公元5世纪）的数学著作《孙子算经》卷下第二十六题，叫做“物不知数”问题，原文如下：
有物不知其数，三三数之剩二，五五数之剩三，七七数之剩二。问物几何？即，一个整数除以三余二，除以五余三，除以七余二，求这个整数。《孙子算经》中首次提到了同余方程组问题，以及以上具体问题的解法，因此在中文数学文献中也会将中国剩余定理称为孙子定理。

2.公式

3. 例子

• 由以下方程组求x
  $$\{\begin{array}{l}x\equiv 1mod2\\ x\equiv 2mod3\\ x\equiv 3mod5\\ x\equiv 5mod7\end{array}$$
  解:
  $$M=\prod _{i=1}^k{m}_i=2\ast 3\ast 5\ast 7=210$$
  $$\begin{gathered} M_1=\frac{M}{m_1}=105，M_1=\frac{M}{m_2}=70 \\ {M}_1=\frac{M}{m_3}=42，M_1=\frac{M}{m_4}=30 \end{gathered}$$
  $$\begin{gathered} e_1\equiv M_1^{-1}mod2\equiv 1，e_2\equiv M_2^{-1}mod3\equiv 1 \\ {e}_3\equiv M_3^{-1}mod5\equiv 3，e_4\equiv M_4^{-1}mod7\equiv 4 \end{gathered}$$
  $$xmod210\equiv (105\ast 1\ast 1+70\ast 1\ast 2+42\ast 3\ast 3+30\ast 4\ast 5)mod210\equiv 173$$

4. 代码实现 (Java)

public class crt {
	
	public static void main(String[] args) {
		// mod 左边数的集合数组		
		int[] a = {1,2,3,5};
		// mod 右边数的集合数组	
		int[] m = {2,3,5,7};
		int x = crt(a,m);
		System.out.println("方程组的解为" + x);
	}
	
	/**
	 * 中国剩余定理(CRT)
	 * @param a
	 * @param m
	 * @return
	 */
	public static int crt(int[] a,int[] m) {
		int M = arrayProduct(m);
		int[] e = new int[m.length];
		int sum = 0;
		for (int i = 0; i < m.length; i++) {
			e[i] = gcd(M/m[i],m[i]);
			sum += M/m[i]*e[i]*a[i];;
		}
		return sum % M;
	}
	
	/**
	 * 整数数组的元素之积
	 * @param array
	 * @return
	 */
	public static int arrayProduct(int[] array) {
		int product = 1;
		for (int i = 0; i < array.length; i++) {
			product *= array[i]; 
		}	
		return product;
	}

	/**
	 * 乘法逆元（ f^-1 mod d）
	 * @param f
	 * @param d
	 * @return
	 */
	public static int gcd(int f, int d) {
		int x1 = 1, x2 = 0, x3 = d;
		int y1 = 0, y2 = 1, y3 = f;
		while (true) {
			if (y3 == 0) {
				return 0;
			}
			if (y3 == 1) {
				return (y2 + d) % d;
			}
			int q = x3 / y3;
			int t1 = x1 - q * y1;
			int t2 = x2 - q * y2;
			int t3 = x3 - q * y3;

			x1 = y1;
			x2 = y2;
			x3 = y3;
			y1 = t1;
			y2 = t2;
			y3 = t3;
		}
	}
}