数论基础知识小结

窝来大致说一下数论的一些东西：
1.欧几里得相关的两个定理：
欧几里得原理：$(a,b) = (b,a \ mod\ b)$
证明：

$$设r = a \ mod \ b,then\ a = kb + r$$
$$设c为a,b的任意公因数$$
则： $$c | a;c | b;c | (a - kb)\ <=> \ c | r$$
$$所以c也是b,r的因数$$
$$同理，我们设d为b,r的任意因数，则可知d也是a,b的公因数。$$
$$设A = \{ d | a \ and\ d | b \} ,B = \{ d | r \ and \ d | b \},\ then A = B$$
$$所以我们可以知道(a,b) = (b,r)$$
证明完毕！

代码如下：

int gcd(int x,int y){return y ? (y,x % y) : x;}

---

扩展欧几里得算法：

$$贝祖等式x,y\in Z,则ax + by = (a,b)有解。$$
证明：
我们可以直接用扩展欧几里得算法步骤来证明：
根据 $(a,b) = (b,a \ mod \ b)$
则设:
$$ax_1 + by_1 = (a,b)$$
$$bx_2 + (a\ mod\ b)y_2 = (b,a\ mod\ b)$$
$$所以：ax_1 + by_1 = bx_2 + (a\ mod\ b)y_2$$
根据取模的定义，我们知道 $$a\ mod \ b = a - \left \lfloor \frac{a}{b} \right \rfloor \cdot b$$
那么可以把等式右边替换：
$$ax_1 + by_1 = bx_2 + ( a - \left \lfloor \frac{a}{b} \right \rfloor \cdot b)y_2$$
我们知道，如果这样一直递归下去，那么肯定会使得y的系数为0。
这样子就能证明有解了。
那么怎么计算出方程的值呢？
我们根据递归出口：
$$when \ b == 0,then\ x = 1,y = 0$$
理由如下：
$$ax + by = (a,b) = a，则x = 1,y = 0必然是一组解。$$
回到刚刚得问题上：
$$ax_1 + by_1 = bx_2 + ( a - \left \lfloor \frac{a}{b} \right \rfloor \cdot b)y_2$$
我们现在已经知道了右边等式的系数都是多少，现在怎么知道 $x_1和y_1的值呢?$
移项：
$$a(x_1-y_2) + b(y_1 - x_2 + \left\lfloor \frac{a}{b} \right \rfloor \cdot y_2) = 0$$
那么显然有：
$$x_1 = y_2 ; y_1 = x_2 - \left\lfloor \frac{a}{b} \right \rfloor\cdot y_2$$
那么我们可以得出算法代码：

void exgcd(int a,int b,int &x,int &y){b ? (exgcd(b,a % b,y,x),y -= a / b * x): (x = 1,y = 0);}

---

中场休息时间QAQ

---

中国剩余定理：
内容如下：
假如有n个同余方程：

$$x \equiv a_i (mod \ m_i)，且m_i两两互素。$$
设 $M = \prod m_i$
求其在模 $M$意义下的解。
我们假设 $$M_i = M / m_i$$
再设 $t_i为M_i在模m_i下的逆元$，即 $t_i M_i = 1(mod \ m_i)$
那么我们就可以方程的解：
$$x = \sum a_it_iM_i$$
证明：
先证明x合法。
$$即证明 \forall i \in[1,n],x \equiv a_i (mod \ m_i)。$$
$$x = \sum _{j \not = i}a_jt_jM_j + a_it_iM_i$$
则在 $mod \ m_i 的意义下：$
$$x = \sum 0 + a_i = a_i (mod \ m_i)$$
这样x是合法的。
实际上，我们可以肯定x是小于M的。
因为如果 $x = \sum a_it_im_i > M合法，那么显然可以把M减去对吧。$
我们这样就得到 $模M$意义下的一个解了.
但是……
注意一下实际上我们是用CRT去求解同余方程组的。
所以我们应该没有必要那么麻烦。
我们考虑把同余方程合并：
对于 $$x\equiv r_1 \mod m_1$$
$$x\equiv r_2 \mod m_2$$
那么我们想把它合成为1个方程。
$$k_1m_1+r_1 = k_2m_2+r_2$$
$$k_1m_1-k_2m_2 = r_2 - r_1$$
可以扩欧来求，对吧。
当 $d = gcd(m_1,m_2) \not | (r_2 - r_1)$时，无解。
我们想让这个 $k_1$尽量小，并且是正整数。
那么显然的恒等式是：
$$k_1m_1+k_2m_2 + lcm(k_1,k_2) - lcm(k_1,k_2)= r_2 - r_1$$
$$k_1(m_1\% \frac {k_2}{d})+k_2(m_2 + p) = r_2 - r_1$$
这个 $p$显然是补足$m_1$少了的那部分的，具体值懒得算了。
所以我们可以注意到，这个 $m_1$显然可以在模意义下去搞。
这样我们就求出了一个 $x$,注意到，$x + lcm(m_1,m_2)$显然是通解，那么实际上我们的方程就变成了：
$$x = r_2 - r_1(mod \ lcm(m_1,m_2))$$
这样一个一个合并就好了。

#include <cstdio>
#include <algorithm>
#include <cstring>
#define Rep(i,n) for(int i = 1;i <= n;++ i)
using namespace std;
typedef long long LL;
LL m[1002],r[1002];
int n;
LL gcd(LL a,LL b){return b ? gcd(b,a % b) : a;}
void exgcd(LL a,LL b,LL &x,LL &y){b ? (exgcd(b,a % b,y,x),y -= a / b * x) : (x = 1,y = 0);}
LL crt()
{
    LL R,M,k1,k2;
    R = r[1],M = m[1];
    for(int i = 2;i <= n;++ i)
    {
        LL dt = r[i] - R,d = gcd(m[i],M);
        if(dt % d)return -1;
        exgcd(M / d,m[i] / d,k1,k2);
        k1 = (k1 * dt / d) % (m[i] / d);
        R = R + k1 * M;
        M = M / d * m[i];
        R %= M;
    }
    return R < 0 ? R + M : R;
}
int main()
{
    while(~scanf("%d",&n))
    {
        Rep(i,n)scanf("%lld%lld",&m[i],&r[i]);
        printf("%lld\n",crt());
    }
    return 0;
}