中国剩余定理(孙子定理)

中国国剩余定理，又名孙子定理，是中国古代求解一次同余式组的方法。当前可用于解决密码学中的秘密共享问题，也可以用来解决布隆过滤器的false positive问题。

1 中国剩余定理

如果，有如下一元线性同余方程组$(p_i为素数)$
$$\{\begin{array}{c}x\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{p}_1)\\ x\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{p}_2)\\ .\\ .\\ .\\ x\equiv a_n(\mathrm{m}\mathrm{o}\mathrm{d}{p}_n)\end{array}\text{\hspace{1em}(1)}$$
方程的解为
$$x=(\sum _{i=1}^n{a}_i{t}_i{M}_i)+kM\text{\hspace{1em}(2)}$$
其中$M={\prod }_{i=1}^n{p}_i,M_i=M/p_i,t_i=M_i^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}{p}_i)$
显然解$x$满足上面的一次同余方程组，所以$x$是方程组的解。接下来只需要证明同余方程组的解的周期为$M$即可。
假设有两个解$x_1,x_2$，那么$x_1-x_2\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}{p}_i)$，可得$x_1-x_2$整除$M$，所以解得周期为$M$，得证.

2 扩展欧几里得算法

辗转相除法可以很方便的求解两个数的最大公约数，如果把求解过程中记录下来就可以用于求解方程$ax+by=gcd(a,b)$，这个方法可以用于求解上面涉及到的$t_i$，$t_i=M_i^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}{p}_i)$即
$$t_i\times M_i\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}{p}_i)\text{\hspace{1em}(3)}$$
所以只需要求解
$$x{M}_i+y{p}_i=1=gcd(M_i,p_i)\text{\hspace{1em}(4)}$$
$$t_i=x\text{\hspace{1em}(5)}$$

3 参考代码

#include <stdio.h>
#include <stdlib.h>

typedef int int32;
// ax + by = gcd(a, b)
int32 exgcd(int32 a, int32 b, int32 *x, int32 *y) 
{
	if(b == 0) {
		*x = 1;
		*y = 0;
		return a;
	} else {
		int32 r = exgcd(b, a % b, x, y);
		int32 t = *x;
		*x = *y;
		*y = t - a / b * (*y);
		return r;
	}
}

int32 func(int32 *a, int32 *p, int32 len)
{
	int32 i, *t;
	int32 M = 1, *m, temp, result = 0;

	m = (int32 *)malloc(len * sizeof(int32));
	t = (int32 *)malloc(len * sizeof(int32));
	for(i = 0; i < len; i++) {
		m[i] = 1;
		M *= p[i];
	}
	printf("\n\nM = %d\nMi ", M);
	for(i = 0; i < len; i++) {
		m[i] = M / p[i];
		printf("%5d  ", m[i]);
		exgcd(m[i], p[i], &t[i], &temp);
	}
	printf("\nTi ");
	for(i = 0; i < len; i++) {
		result += a[i] * t[i] * m[i];
		printf("%5d  ", t[i]);
	}
	printf("\n");
	while (result <= 0)
		result += M;
	free(m);
	free(t);
	return result % M;
}

int main()
{
	int32 i, len = 4, result;
	int32 a[10] = {1, 2, 4, 1}, p[10] = {5, 3, 7, 11};
	printf("primer : ");
	for(i = 0; i < len; i++) {
		printf("%3d ", p[i]);
	}
	printf("\nremaind: ");
	for(i = 0; i < len; i++) {
		printf("%3d ", a[i]);
	}
	result = func(a, p, len);
	printf("\nresult = %d\n\n", result);
	system("pause");
}