CRT和EXCRT学习笔记

蒟蒻maomao终于学会\(CRT\)啦!发一篇博客纪念一下(还有防止忘掉)

\(CRT\)要解决的是这样一个问题:

$$x≡a_1​(mod m_1​)$$

$$x≡a_2​(mod m_2​)$$

$$x≡a_3​(mod m_3​)$$

$$...$$

$$x≡a_k​(mod m_k​)​$$

其中,\(m\)之间两两互质。这个问题有一个通解是\(\sum a_i * M * t_i / m_i\),其中\(t_i\)代表方程\(M * t_i / m_i ≡ 1\)的最小正整数解。

为什么它是对的呢?对于任意一个式子\(x≡a_j(mod m_j)\),通解中\(i = j\)的部分会贡献\(a_i\)的余数,而其它部分会贡献\(0\)的余数。

更一般的,我们来考虑如果\(m\)之间不互质的情况,由于打公式很累,所以详细请参考这个博客

发一下\(exCRT\)的板子。

#include <bits/stdc++.h>
using namespace std;

#define int long long


const int N = 100010;

int n, bi[N], ai[N];

int add (int a, int b, int mod) {
	return ((a + b) % mod + mod ) % mod;
}

int mul (int a, int b, int mod) {
    int res = 0;
    while (b > 0) {
        if (b & 1) {
			res = (res + a) % mod;
		}
        a = (a + a) % mod;
        b >>= 1;
    }
    return res;
}

int exgcd (int a, int b, int &x, int &y) {
    if (b == 0) {
		x = 1;
		y = 0;
		return a;
	}
    int gcd = exgcd (b, a % b, x, y);
    int xx = y, yy = x - (a / b) * y;
    x = xx, y = yy;
	return gcd;
}

int excrt () {
    int x, y;
    int M = ai[1], ans = bi[1]; //通解是b[1] + a[1] * t ≡b[2] (mod a[2]); 
	for(int i = 2; i <= n; ++i) {
		//M * x + a[i] * y = b[i] - ans;
		//其中 ans + M * x % lcm (M, b[i]) 就是新的通解 
        //求出来的x是对于gcd (M, a[i])而言,所以要乘上c / gcd (M, a[i]); 
		int a = M, b = ai[i], c = add (bi[i], -ans, b);
        int gcd = exgcd (a, b, x, y), bg = b / gcd;
        x = mul (x, c / gcd, ai[i]);
        ans += x * M;//更新前k个方程组的答案
        M *= bg;//M为前k个m的lcm
        ans = (ans %M + M) % M;
    }
    return ans;
}

signed main () {
    cin >> n;
	for (int i = 1; i <= n; ++i) {
		cin >> ai[i] >> bi[i]; //b是余数,a是模数。 
	}
	cout << excrt () << endl;
    return 0;
}
posted @ 2019-03-06 16:50  maomao9173  阅读(195)  评论(0编辑  收藏  举报