2021-02-10

中国剩余定理例题

（其实和标题没啥关系，重点在于分析

题意

https://www.acwing.com/problem/content/206/

给定 2n 个整数a1,a2,…,an和m1,m2,…,mn,求一个最小的非负整数 x，满足∀i∈[1,n], x≡mi(mod ai)。

输入格式
第1 行包含整数 n。

第 2…n+1行：每 i+1 行包含两个整数ai和mi,数之间用空格隔开。

输出格式
输出最小非负整数 x，如果 x 不存在，则输出 −1。
如果存在 x，则数据保证 x 一定在64位整数范围内。

数据范围
1≤ai≤2^31−1,
0≤mi<ai
1≤n≤25

题解

$取两个等式，x\equiv m_1(mod{a}_1)和x\equiv m_2(mod{a}_2),$
$即x=k_1\times a_1+m_1和x=k_2\times a_2+m_2，$
$即k_1\times a_1+m_1=k_2\times a_2+m_2，$
$即k_1\times a_1-k_2\times a_2=m_2-m_1,$
$即a\times x+b\times y=d形式$

$两式子有解，即m_2-m_1是d的倍数（d=gcd(a,b)=gcd(a_1,a_2)）$

$k1先随着d扩大到m_2-m_1,即k1\times \frac{m_2-m_1}{d}$
$解出k_1的通解为k_1=k_1+\frac{a_2}{d}\times kk,kk为整数$
$所求x=a_i\times k_i+m_i=a_1\times k_1+m_1,$
$通解为x=a_1\ast k_1+m_1+\frac{a_1\ast a_2}{d}\times kk，$
$（a_1、m_1给出，k_1求出，kk未知(由于k_1有题解所以有kk)）$

$要使得x为最小非负整数,那么k1要小，且保证为非负整数，$
$即k1=(k1\%abs(\frac{a_2}{d})+abs(\frac{a_2}{d}))\%abs(\frac{a_2}{d})$

概括一下：（思路）
$有n个等式，先任意取两个，然后算出新的等式，看成一个，$
$再与后面的新的式子，即x相等，前两个等式化成的等式$
$x=a_1\ast k_1+m_1+\frac{a_1\ast a_2}{d}\times kk应该写成x=a_i\times k_i+m_i的格式$
$在上面的式子中，只有kk是未知数，$
$所以写成x=\frac{a_1\ast a_2}{d}\times kk+(a_1\ast k_1+m_1)$
$该式子为两个式子合并后的式子，由此与新的式子联立求解$

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int INF = 1e9 + 7;
const int MAX = 1e5 + 10;
int a[MAX];
int m[MAX];

ll exgcd(ll a, ll b, ll& x, ll& y) {
   if(b == 0) {
       x = 1, y = 0;
       return a;
   }
   ll d = exgcd(b, a % b, y, x);
   y -= a / b * x;
   return d;
}

int main() {
   int n;
   scanf("%d", &n);
   ll a1, m1, a2, m2;
   scanf("%lld%lld", &a1, &m1);
   bool flag = 1;
   ll x = 0;
   for(int i = 1; i < n; i++) {
       scanf("%lld%lld", &a2, &m2);
       ll k1, k2;
       ll d = exgcd(a1, -a2, k1, k2);
       if((m2 - m1) % d != 0) { flag = 0; }
       ll k11 = k1 * (m2 - m1) / d ;
       k1 = (k11 % abs(a2 / d) + abs(a2 / d)) % abs(a2 / d);
       x = a1 * k1 + m1;
       m1 = a1 * k1 + m1;
       a1 = abs(a1 / d * a2);
   }
   if(flag == 0) printf("-1\n");
   else printf("%lld\n", x);
}