扩展中国剩余定理模板

这就是神奇的非互质版CRT，同余方程组中各个模数可能不互质，如果再用以前互质版的做法就会出错（这个的原因我也没有深究，如果有兴趣的小伙伴懂的话欢迎在评论区里指点本蒟蒻）。
对于模数不互质的情况，需要逐个合并方程求解，具体证明如下图：

还是FZU 1402的代码, 非互质版算法同样适用于模数互质的题目，有兴趣的同学可以做一下hdu 1573
FZU 1402：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<climits>
using namespace std;
typedef long long ll;
const int maxn=14;
int n;
ll a[maxn],m[maxn];
ll gcd(ll a,ll b) {return !b?a:gcd(b,a%b);
}
ll exgcd(ll a,ll b,ll &x,ll &y) {if (!b) {x=1,y=0;return a;}ll d=exgcd(b,a%b,y,x);y-=(a/b)*x;return d;
}
ll inv(ll a,ll m) {ll x,y;ll d=exgcd(a,m,x,y);if (d==-1) return -1;return (x%m+m)%m;
}
bool merge(ll a1,ll m1,ll a2,ll m2,ll &a3,ll &m3) {ll d=gcd(m1,m2),c=a2-a1;if (c%d) return false;c=(c%m2+m2)%m2,c/=d,m1/=d,m2/=d,c*=inv(m1,m2),c=(c%m2+m2)%m2,c=(c*m1*d)+a1;m3=m1*m2*d;a3=(c%m3+m3)%m3;return true;
}
ll crt() {ll a1=a[1],m1=m[1];for (int i=2;i<=n;i++) {ll aa,mm;if (!merge(a1,m1,a[i],m[i],aa,mm)) return -1;a1=aa,m1=mm;}return (a1%m1+m1)%m1;
}
int main() {while(scanf("%d",&n)!=EOF) {for (int i=1;i<=n;i++) scanf("%I64d%I64d",&m[i],&a[i]);printf("%I64d\n",crt());}return 0;
}

扩展中国剩余定理模板相关推荐

中国剩余定理扩展中国剩余定理模板
中国剩余定理解线性同余线性方程 /*long long gcd(LL a,LL b) {return b==0?a:gcd(b,a%b); }*/ #include<bits/stdc++.h& ...
扩展中国剩余定理(模板)
证明:挂个大佬的博客:点击查看挂上代码供以后白嫖: LL n,m[N],a[N];///模数为m,余数为a, X % m = aLL extend_gcd(LL a,LL b,LL &x,L ...
中国剩余定理扩展中国剩余定理（模板）
中国剩余定理 && 扩展中国剩余定理一个整数除以三余二,除以五余三,除以七余二,求这个整数. 例题: 一个正整数K,给出K Mod 一些质数的结果,求符合条件的最小的K.例如,K % ...
中国剩余定理（CRT）扩展中国剩余定理(exCRT)
前言中国剩余定理(也叫孙子定理)并不是很复杂,由于最近用到了,以前学的时候还不写博客,所以现在补一下中国剩余定理(CRT) 问题给出nnn个同余方程 x≡a1(modp1)x≡a2(modp2) ...
[数论]-----中国剩余定理（扩展中国剩余定理）
中国剩余定理中国剩余定理(CRT)用于求形如: { x ≡ a 1 ( m o d m 1 ) x ≡ a 2 ( m o d m 2 ) ⋯ ⋯ x ≡ a k ( m o d m k ) \be ...
中国剩余定理（CRT）和扩展中国剩余定理（EXCRT）
Tip:建议读者不要太着急后翻,按照顺序阅读有助于理解中国剩余定理(CRT) 问题引出 "有物不知其数,三三数之剩二,五五数之剩三,七七数之剩二.问物几何?"即,一个整数除以三余 ...
数学--数论--中国剩余定理+扩展中国剩余定理（孙子定理）
中国剩余定理问题求解同余方程组其中m1,m2,m3...mkm_1,m_2,m_3...m_km1,m2,m3...mk为两两互质的整数求x的最小非负整数解定理令M=∏i=1kmi ...
模线性方程(中国剩余定理+扩展中国剩余定理)
已知一系列除数和模数,求最小的满足条件的数我们先考虑一般的情况,即模数不互质.(扩展中国剩余定理) 我们考虑两个方程的情况 x%M=R x=k1∗M+Rx=k1 * M+Rx=k1∗M+R x%m= ...
中国剩余定理扩展中国剩余定理入门详解
中国剩余定理例题已知以下 n n n同余方程(所有 m i m_i mi互质): x ≡ a 1 ( m o d m 1 ) x≡a_1(\mod m_1) x≡a1(modm1) x ≡ ...

扩展中国剩余定理模板

扩展中国剩余定理模板相关推荐

最新文章

热门文章