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

qq_38232157

于 2022-11-01 22:50:30 发布

阅读量161

点赞数

分类专栏： noi 后缀数组文章标签：扩展中国剩余定理

本文链接：https://blog.csdn.net/qq_38232157/article/details/127642968

版权

本文介绍了扩展中国剩余定理（EXCRT）的基础知识，包括如何解决带有多个同余式子的问题。当数据规模较大，可能导致长整型溢出时，建议使用__int128类型来确保计算正确，并在输出时进行类型转换。

摘要由CSDN通过智能技术生成

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 <bits/stdc++.h>
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)