求解同余方程组（难度：2颗星）

问题描述：

有一个同余方程组，有N个同余方程组成（N由用户输入），另外每个同余方程的a[i]和m[i]也又用户指定，如下所示：

x≡a[1](mod m[1])
x≡a[2](mod m[2])
x≡a[3](mod m[3])
x≡a[4](mod m[4])
…
x≡a[n](mod m[n])

求x的最小正整数解。

问题分析：

参考我的另外一篇文章“中国剩余定理”
http://blog.csdn.net/yi_ming_he/article/details/72848318

参考代码：

#include <stdio.h>

typedef __int64 int64;
int n, a[100], m[100];//m[100]是mod数组

int64 ExtendGcd(int64 a, int64 b, int64* px, int64* py)
{
    int64 nRet, temp;
    if (b == 0)
    {
        if (px && py)
        {
            *px = 1;
            *py = 0;
        }
        return a;
    }

    nRet = ExtendGcd(b, a % b, px, py);
    if (px && py)
    {
        temp = *px;
        *px = *py;
        *py = temp - a / b * (*py);
    }
    return nRet;
}

int64 ChinaRemainder()
{//中国剩余定理只适合解决m数组两两互质的情况
    int64 M = 1, res = 0, x, y, tempM;
    int i;
    for (i = 0; i < n; i++)
        M *= m[i];

    for (i = 0; i < n; i++)
    {
        tempM = M / m[i];
        ExtendGcd(tempM, m[i], &x, &y);
        res = (res + tempM * x * a[i]) % M;
    }
    res = (res + M) % M;
    return res;
}

int main()
{
    scanf_s("%d", &n);//同余方程组的方程个数
    int i;
    for (i = 0; i < n; i++)
        scanf_s("%d%d", &a[i], &m[i]);

    printf("%I64d\n", ChinaRemainder());
    return 0;
}

运行结果：