我感觉中国剩余定理我并不能顺推出来,这是一个很巧妙的构造,这个构造思想值得学习,但是扩展的可以顺推。
问题:
有一
x
x
x,满足:
x
≡
c
1
(
m
o
d
m
1
)
x≡c1(mod~m1)
x≡c1(mod m1)
x
≡
c
2
(
m
o
d
m
2
)
x≡c2(mod~m2)
x≡c2(mod m2)
…
x
≡
c
k
(
m
o
d
m
k
)
x≡ck(mod~mk)
x≡ck(mod mk)
∀ i ≠ j , ( m i , m j ) = 1 ∀i≠j,(mi,mj)=1 ∀i=j,(mi,mj)=1
求 x x x的通解形式。
设
M
=
∏
m
i
,
M
i
=
M
m
i
,
t
i
=
M
i
−
1
(
m
o
d
m
i
)
M=\prod mi,Mi={M \over mi},ti=Mi^{-1}(mod~mi)
M=∏mi,Mi=miM,ti=Mi−1(mod mi)
则
x
≡
∑
M
i
∗
t
i
∗
c
i
(
m
o
d
M
)
x≡\sum Mi*ti*ci(mod~M)
x≡∑Mi∗ti∗ci(mod M)
这个怎么理解?
对一个确定的 M i ∗ t i Mi*ti Mi∗ti,我们让它 m o d mod mod一下 m i mi mi,肯定等于1, m o d m j ( j ≠ i ) mod~mj(j≠i) mod mj(j=i),肯定为0,类似于拉格朗日插值。
扩展中国剩余定理比中国剩余定理应用广的多,也很好写。
问题为:
x
≡
c
1
(
m
o
d
m
1
)
x≡c1(mod~m1)
x≡c1(mod m1)
x
≡
c
2
(
m
o
d
m
2
)
x≡c2(mod~m2)
x≡c2(mod m2)
m 1 , m 2 m1,m2 m1,m2不一定互质,求是否有解?有解的话,通解形式是什么?
设 M = l c m ( m 1 , m 2 ) M=lcm(m1,m2) M=lcm(m1,m2),如果能够合并,那么显然是合并成: x ≡ c 3 ( m o d M ) x ≡ c3 (mod~M) x≡c3(mod M)的形式。
设 t = g c d ( m 1 , m 2 ) t=gcd(m1, m2) t=gcd(m1,m2)
设:
c
3
=
c
1
+
a
∗
m
1
(
0
≤
a
<
(
M
/
m
1
=
m
2
/
t
)
)
c3=c1+a*m1(0 \le a < (M/m1=m2/t))
c3=c1+a∗m1(0≤a<(M/m1=m2/t))
c
3
=
c
2
+
b
∗
m
2
(
0
≤
a
<
(
M
/
m
2
=
m
1
/
t
)
)
c3=c2+b*m2(0 \le a < (M/m2=m1/t))
c3=c2+b∗m2(0≤a<(M/m2=m1/t))
现在的目标是求出 a a a或者 b b b,这样就可以求出 c 3 c3 c3了。
联立得:
c
1
+
a
∗
m
1
=
c
2
+
b
∗
m
2
c1+a*m1=c2+b*m2
c1+a∗m1=c2+b∗m2
c
2
−
c
1
=
a
∗
m
1
−
b
∗
m
2
c2-c1=a*m1-b*m2
c2−c1=a∗m1−b∗m2
根据扩展欧几里得那一套, a ∗ m 1 − b ∗ m 2 a*m1-b*m2 a∗m1−b∗m2只能表示 t t t的倍数,所以 t ∣ ( c 2 − c 1 ) t|(c2-c1) t∣(c2−c1)才会有解。
方程两边同时除以
t
t
t,得:
(
c
2
−
c
1
)
/
t
=
a
∗
m
1
/
t
−
b
∗
m
2
/
t
(c2-c1)/t=a*m1/t-b*m2/t
(c2−c1)/t=a∗m1/t−b∗m2/t
假设求 a a a,注意到 0 ≤ a < ( M / m 1 = m 2 / t ) 0 \le a < (M/m1=m2/t) 0≤a<(M/m1=m2/t),不放对整个方程 m o d ( M / m 1 = m 2 / t ) mod ~(M/m1=m2/t) mod (M/m1=m2/t), b ∗ m 2 / t b*m2/t b∗m2/t就没有了。
得:
(
c
2
−
c
1
)
/
t
=
a
∗
m
1
/
t
(
m
o
d
m
2
/
t
)
(c2-c1)/t=a*m1/t~(mod~{ m2/t})
(c2−c1)/t=a∗m1/t (mod m2/t)
a
=
(
i
n
v
(
m
1
/
t
,
m
2
/
t
)
∗
(
c
2
−
c
1
)
/
t
m
o
d
(
m
2
/
t
)
)
a=(inv(m1/t,m2/t)*(c2-c1)/t~~mod~~(m2/t))
a=(inv(m1/t,m2/t)∗(c2−c1)/t mod (m2/t))
求逆元用 e x g c d exgcd exgcd,就可以得到 a a a,进而得到 c 3 c3 c3。
M ≤ 1 0 18 M \le 10^{18} M≤1018时小心爆 l o n g l o n g long long longlong。
#include<cstdio>
#define pp printf
#define ll long long
#define fo(i, x, y) for(int i = x; i <= y; i ++)
using namespace std;
const int N = 1e5 + 5;
int n;
ll m1, m2, c1, c2;
ll mul(ll x, ll y, ll mo) {
x %= mo; y %= mo;
ll z = (long double) x * y / mo;
z = x * y - z * mo;
if(z < 0) z += mo; else
if(z > mo) z -= mo;
return z;
}
ll gcd(ll x, ll y) {
return !y ? x : gcd(y, x % y);
}
void eg(ll a, ll b, ll &x, ll &y) {
if(!b) {x = a, y = 0; return;}
eg(b, a % b, y, x); y = y - (a / b) * x;
}
ll inv(ll v, ll p) {
ll x, y;
eg(v, -p, x, y);
x = (x % p + p) % p;
return x;
}
int main() {
scanf("%d", &n);
c1 = 0, m1 = 1;
fo(i, 1, n) {
scanf("%lld %lld", &m2, &c2);
ll t = gcd(m1, m2); ll m3 = m2 / t, M = m3 * m1;
c1 = (mul(mul(inv(m1 / t, m3), (c2 - c1) / t, m3), m1, M) + c1) % M;
m1 = M; if(c1 < 0) c1 += m1;
}
pp("%lld\n", c1);
}