洛谷 P4777 【模板】扩展中国剩余定理(EXCRT)

1、 中国剩余定理 (n条同余式子, 前提是 m[1] ~ m[n] 两两互质)
x=r[1](mod m[1])
x=r[1](mod m[2])

x=r[n](mod m[n])
2、 扩展中国剩余定理 (n条同余式子, m[1] ~ m[n] 不一定两两互质)
x=r[1](mod m[1])
x=r[1](mod m[2])

x=r[n](mod m[n])

	考虑签名两条方程, x=r[1](mod m[1]), x=r[1](mod m[2])
	转化为不定方程  x = m[1] * p + r[1] = m[2] * q + r2
	那么 m[1] * p - m[2] * q = r[2] - r[1]
	gcd(m[1], m[2]) 整除 r[2] - r[1] 说明有解,否则无解
	由扩展欧几里得算法求出方程 m[1] * p - m[2] * q = r[2] - r[1] 一组特解,
	再得到通解,最后这两条方程等价合并为一个新的方程式
	x = r(mod m),  再用这新的方程 去和第 3条方程联解,以此类推, 
	n个同余方程,合并 n - 1 

3、如果数据用 long long ,最后一个测试点会爆 long long .
直接用 __int128, 输出的时候,用强转

#include 
using namespace std;
typedef __int128 ll;
const int MaxN = 1e5 + 10;
int n;
ll r[MaxN], m[MaxN];

ll exgcd(ll a, ll b, ll& x, ll& y)

你可能感兴趣的:(noi,后缀数组,扩展中国剩余定理)