CRT和扩展CRT学习笔记

CRT（中国剩余定理）（孙子定理）

一般的CRT只能用于解决模数彼此之间互质的情况。（以下假设$b$表示模数，$a$表示$x$在$modb$意义下的值）

大概想法是构造$M={\prod }_{i=1}^k{b}_i$，对于第$i$个方程，构造一组特殊的解$ti$满足$\frac{M}{b_i}\ast t_i=1(mod{b}_i)$。最后答案就是${\sum }_{i=1}^k\frac{M}{b_i}\ast t_i\ast a_i$。

正确性证明就考虑对于每个方程，不对应这个方程的项$mod{b}_i$为0，恰好对应这个方程的项为$a_i$。恰好构造出一组解。

ex_gcd寻找通解

对于$ax+by=c(modp)$，最后$b=0$的时候：$x=c/a,y=0；$
对于$ax+by=gcd(a,b)$，最后$b=0$的时候：$x=1,y=0.$

直接扩欧求出一组特值以后。假设$x$的变化量为$\delta$，那么$b|a\delta$。如果$gcd(a,b)=1$的话，$b|\delta$，直接$((modb)+b)modb$就好了。否则$a,b$同时除以$gcd(a,b)$。答案$x$的变化一定是$x-T\ast b^{\prime }$的形式了。

$luoguP3868$

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef long double LDB;
LL a[15], b[15], M = 1;
LL ans;
int n;
inline void ex_gcd(LL a, LL b, LL c, LL &x, LL &y) {
    if (!b) {
        x = c / a;
        y = 0;
        return ;
    }
    ex_gcd(b, a % b, c, y, x);
    y -= (a / b) * x;
}
inline LL q_c(LL x, LL y, LL P) {
    return (x * y - (LL)((LDB)x * y / P) * P);
}
int main() {
    scanf("%d", &n);
    for (int i = 1; i <= n; ++i) scanf("%lld", &a[i]);
    for (int i = 1; i <= n; ++i) scanf("%lld", &b[i]), M *= b[i];
    for (int i = 1; i <= n; ++i) {
        LL x, y;
        ex_gcd(M / b[i], b[i], 1, x, y);
        x = ((x % b[i]) + b[i]) % b[i];
        ans = (ans + q_c(q_c(M / b[i], x, M), a[i], M)) % M;
    }
    printf("%lld\n", ans);
    return 0;
}

扩展CRT和CRT没有什么关系。

扩展CRT是假设前$k-1$个方程的答案已经得出，通解是$x=x_0+t\ast M$的形式。$M$是前$k-1$个$b$的$lcm$。
现在要求一个$x$使得满足第$k$个方程，也就是找到一个$t_0$使得$x_0+t_0\ast M=a_i(mod{b}_i)$。 这个$t_0$可以用扩展欧几里得求出来，然后就求出了满足前$k$个方程的一组特解$x=x_0+t_0\ast M$。此时再把$M$更新就可以了。

具体实现的时候，当前计算的$x$是不可以取模的，往往都会要求答案不超过$longlong$的范围。但是中间计算可能会爆$longlong$。这个时候里面要写$O(1)$快速乘。

另，此题最好应采用$ax+by=gcd(a,b)$形式的扩欧。对于$ax+by=c$的不定方程的解就是$ax+by=gcd(a,b)$的解乘以$\frac{c}{gcd(a,b)}$，如果要调整的话，$x+=\frac{b}{gcd(a,b)},y-=\frac{a}{gcd(a,b)}。$而且一开始扩欧的系数都要是正的。两种CRT都要把答案对M取模！

一个巧合：$ax=1(modb)$的同余方程，即$ax+by=1(modb)(gcd(a,b)=1)$。如果扩欧解出了这个的一组解$(x0,y0)$的话。$x$的调整系数恰好是$b/1=b$。造成可以直接模$b$的效果。实际上，对于一个任意的不定方程，调整应该如上调整。

板子： $luoguP4777$

$luoguP2480$ 猪文

题意非常清楚，就是求$G^{{\sum }_{d|n}{C}_n^{\frac{n}{d}}}(mod999911659)$的值。

$n$的范围可以$O\sqrt{(}n)$枚举因子。对于每个因子可能我们会想用$lucas$算出这个组合数。但是有$3$个问题。第一,$lucas$要求模数必须是质数。 第二，$lucas$的模数不能太大，因为我们要暴力算$P$以内的组合数。第三，就算我们算出了G的幂次，也无法直接快速幂。因为$mod$运算是不能带指数的。

以下记$P=999911659$

学的定理多果然是有好处的。根据扩展欧拉定理，$a^b=a^{bmod\varphi (n)}(modn)$，要求$gcd(a,n)=1$。所以我们先判断$G$是不是$P$的倍数，如是的话就为$0$了。否则，由于$P$是一个质数，定有$gcd(a,P)=1$。就可以把上面取模了。

又有一个问题，$P-1$不是一个质数，我们可能要对上面暴力扩展$lucas$。但是这样$P-1$ 还是太大了不好做。 （mark一下没学过扩展$lucas$）

我们把$P-1$分解质因数。$999911658=2\ast 3\ast 4679\ast 35617$，所以我们把$modP-1$转化成$mod$这些互质的质因数。这样根据$CRT$，求出符合4个条件的最小正整数解就可以了(记得一边算一边要$mod(P-1)$)。最后一个快速幂就可以算出答案了。