数学基础III 扩展中国剩余定理 扩展Lucas定理 扩展BSGS 原根 指标

扩展中国剩余定理

一般的中国剩余定理只能解决模数互质的情况。扩展中国剩余定理可以解决模数不互质的情况。
考虑依次将$n$个方程合并起来。设前面合并起来的方程是$x=m_u{k}_u+a_u$，当前枚举到的方程是$x=m_v{k}_v+a_v$，则联立得到（$x$的符号可忽略）$$m_u{k}_u+m_v{k}_v=a_u-a_v$$
使用扩展欧几里得，求出最小正整数解$k_u^{\prime }$和$k_v^{\prime }$。那么我们合并得到的方程就是$$x=\mathrm{l}\mathrm{c}\mathrm{m}({\mathrm{m}}_{\mathrm{u}},{\mathrm{m}}_{\mathrm{v}})\mathrm{k}+{\mathrm{m}}_{\mathrm{u}}{\mathrm{k}}_{\mathrm{u}}^{\prime }+{\mathrm{a}}_{\mathrm{u}}$$
这样合并到最后就能得到答案。

long long cm = m[0], ca =a[0];
for(int i = 1; i < n; i++) {
	long long d, x, y, dela = (a[i]-ca%m[i]+m[i])%m[i];
	exgcd(cm, m[i], d, x, y);
	if(dela % d) {
		printf("You've got no electricity, I come to your house!\n");
		return 0;
	}
	long long sol = qmul(x, dela/d, m[i]/d);
	ca += sol * cm;
	cm *= m[i] / d;
	ca = (ca%cm+cm)%cm;
}
return ca;
//2019.3.14 // mod = 苯酚醛

20180929UPD：由$$m_1{k}_1=a_2-a_1+m_2{k}_2$$
两边同时除以$\gcd (m_1,m_2)$，令其为$G$（注意判断是否$a_2-a_1|G$，不成立则方程无解）。令$A=\frac{a_2-a_1}{G}$。
$$\frac{m_1}{G}{k}_1=A+\frac{m_2}{G}{k}_2$$
$$\frac{m_1}{G}{k}_1=A(\mathrm{m}\mathrm{o}\mathrm{d}\frac{m_2}{G})$$
$$k_1=A\cdot inv(\frac{m_1}{G},\frac{m_2}{G})(\mathrm{m}\mathrm{o}\mathrm{d}\frac{m_2}{G})$$
$$k_1=A\cdot inv(\frac{m_1}{G},\frac{m_2}{G})+\frac{m_2}{G}\cdot k^{\prime }$$
最后将$k_1$代为$\frac{x-a_1}{m_1}$
$$x=A\cdot inv(\frac{m_1}{G},\frac{m_2}{G})\cdot m_1+a_1(\mathrm{m}\mathrm{o}\mathrm{d}\frac{m_1{m}_2}{G})$$
这就是合并后的方程。

扩展Lucas定理

Lucas定理只能解决模数是质数的情况。事实上很多扩展定理都是将质数推广到非质数而来的。扩展Lucas定理可以解决模数为任意正整数的问题，复杂度为$O(p^t)$，其中$p^t$是模数分解为的质因数的幂。

首先学习求$n!\mathrm{m}\mathrm{o}\mathrm{d}{p}^t$。将乘数中将$p$的倍数提出来，除以倍数后依然是阶乘形式，因此这一部分可以递归计算。而其它部分按值域每$(p-1)p^{t-1}$分为一组，它们对$\mathrm{m}\mathrm{o}\mathrm{d}{p}^t$是同余的，因此计算一组即可，复杂度$O(p^t)$。
举例说明，模数为$3^2时，$$25!\equiv 3^8\times 8!\times (1\times 2\times 4\times 5\times 7\times 8{)}^3\times 25$。
然后就可以计算$C_n^m\mathrm{m}\mathrm{o}\mathrm{d}{p}^t$了。
回到原问题，我们将模数质因数分解为$\prod p_i^{t_i}$，设答案为$ans$，则我们得到若干个同余方程$ans\equiv C_n^m(\mathrm{m}\mathrm{o}\mathrm{d}{p}_i^{t_i})$。使用上述方法算出系数，然后中国剩余定理合并即可。

扩展BSGS

求$a^x\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$，$m$为任意正整数。
因为要求$a^{-1}\mathrm{m}\mathrm{o}\mathrm{d}m$，所以在$\gcd (a,m)=1$时才能用BSGS解。考虑转化成这个形式。
原方程化为$a^{x-1}\frac{a}{\gcd (a,m)}\equiv \frac{b}{\gcd (a,m)}(\mathrm{m}\mathrm{o}\mathrm{d}\frac{m}{\gcd (a,m)})$。将$\frac{a}{\gcd (a,m)}$等看作系数，递归计算直到$\gcd (a,m)=1$。然后反带回去即可。答案$=x+递归次数$。

upd on 2019.1.13
设$d=\gcd (a,p)$，则当$d|p$时有解。每次$p$会除以一个数，因此最多递归$O(\log p)$次。每次$a$的次数会减少$1$，因此最后求出的是$x-递归次数$。注意答案有可能小于递归次数，因此还要枚举小于等于递归次数的答案判断是否合法。

原根 指标

在模$p$意义下的原根$g$有以下性质：
原根个数$=\phi (\phi (p))$；
$g^i\equiv g^j(\mathrm{m}\mathrm{o}\mathrm{d}p)\iff i\equiv j(\mathrm{m}\mathrm{o}\mathrm{d}\phi (p))$；
指标的定义：$g^{ind(x)}\equiv x$；
由此，$x^k(\mathrm{m}\mathrm{o}\mathrm{d}p)$可转化为$ind(x)\cdot k(\mathrm{m}\mathrm{o}\mathrm{d}\phi (p))$。
求指标用BSGS，将指标还原用快速幂。

原根的求法

设模数为$p$，将$\phi (p)$分解为$\prod P_i^{t_i}$。从$2$开始枚举$g$，如果有$g^{\frac{\phi (p)}{P_i}}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$满足则不是原根，否则找到一个原根。
如果广义黎曼猜想成立，$g$取到$\log \log \log \log \log \log n$规模时即可找到一个原根。由一个数的质因数个数是$\frac{\log n}{\log \log n}$，总复杂度为$\frac{\log n\times \log \log \log \log \log \log n}{\log \log n}$，该算法的瓶颈在于分解质因数。

二次剩余

本季度内
from **999
已隐藏