金华集训Day7——宇-CSDN博客

本文链接：https://blog.csdn.net/weihongyu_2023/article/details/98312778

Contents)

一、质因子分解

1.素性测试

①两种常见的素性测试

$\qquad$ 试除法，配合质数筛法可以做到 $O(\sqrt n)$ 的复杂度。
$\qquad$ $M i l l e r - R a b i n$ 素性测试，可以做到 $O(log\ n)$ 的复杂度，属于 $R P$ 算法。
$\qquad$ 其它还有印度人的 $A K S$ 算法可以做到正确率 $100\%$ 的 $O(log\ n)$ 复杂度。
$\qquad$ 绝大多数情况 $M i l l e r - R a b i n$ 都够优秀了。
$\qquad$ 可以查阅 $w i k i$ 以确定 $M i l l e r - R a b i n$ 应该试到多少。

②Miller-Rabin

$\qquad$ 基本原理是飞马小定理：若 $p$ 是质数， $a, p$ 互质，则 $a^{p-1}\equiv1\pmod p$ 。
$\qquad$ 于是对于某个 $p$ ，若能找到与它互质的 $a$ 使得 $a^{p-1}\ne 1\pmod p$ ，则 $p$ 必定不是质数。
$\qquad$ 然而有些合数 $p$ ，满足所有它互质的 $a$ 都有 $a^{p-1}\equiv 1\pmod p$ ，这种数称为 $C a r n i c h a e l$ 数（如 $561 = 3 * 11 * 17$ ），这样的数是用上面的方法检验不出来的。
$\qquad$ 所以还需要奇素数判定。对于奇素数 $p$ ，如果 $a^{p-1}\equiv1\pmod p$ 即 $(a^\frac{p-1}2+1)(a^\frac{p-1}2-1)\equiv0\pmod p$ ，由于 $F_p$ 是整环，所以 $a^\frac{p-1}2\equiv1$ 或 $p - 1$ 。如果 $\frac{p-1}2$ 还是偶数则可以继续往下检验……
$\qquad$ 用原根的一些理伦可以证明这样就能保证对于任意合数至少存在一个 $a$ 可以判定它是合数。

2.质因子分解

①两种常见的质因子分解方法

$\qquad$ 试除法，复杂度 $O(\sqrt n)$ ，太慢了。
$\qquad$ $P o l l a r d^{'} s R h o$ ，期望复杂度 $O(\sqrt[4] n \ log\ n)$ ，有名启发式分解。

②为什么叫启发式分解？

$\qquad$ $n$ 是素数的时候用 $M i l l e r - R a b i n$ 素性测试。
$\qquad$ $n$ 不是素数的时候，复杂度只和最小的质因子有关。所以称为启发式。

③Pollard’s Rho

$\qquad$ 假如要分解一个数 $n$ ：
$\qquad$ 首先进行素性测试，是素数直接返回。
$\qquad$ 否则就要生成一些随机的 $x_i$ ，去求 $gcd(|x_i-x_j|,n)$ ，如果这个 $\in(1,n)$ 则找到了 $n$ 的一个因子，递归分解。
$\qquad$ 一个挺靠谱的随机方法就是 $x←x^2+c$ ， $c$ 是个随机数。
$\qquad$ 这样随机出来的 $x$ 可能会进入循环，假如进入循环了我们还没找到因子，就重新随个 $x$ 和 $c$ ，重新做。

④如何判已经进入循环？

$\qquad$ 可以证明 $x←x^2+c$ ，形成的一定是一个 $p$ 形结构。
$\qquad$ 每次当 i 为 2 的幂次的时候就令 $y←x_i$ ，如果某时刻 $x_i=y$ 了则说明已经在环上绕了一圈了。
$\qquad$ 即：看 $x_{(2,4]}$ 是否 $x_2$ ，看 $x_{(4,8]}$ 是否 $= x 4$ ，看 $x_{(8,16]}$ 是否 $x_8……$
$\qquad$ 这样“浪费”的步数仅仅是 $O (环长)$ 级别的。

二、数论

1.欧几里得算法

①欧几里得求gcd

$\qquad$ 若 $x ∣ a, x ∣ b$ ，则 $x ∣ a + b$
$\qquad$ 于是 $gcd(a,b)=gcd(b,a\%b)$ 。

②扩展欧几里德算法

$\qquad$ 已知 $a, b$ ，求出 $x, y$ 满足 $a x + b y = g c d (a, b)$ 。
$\qquad$ 在欧几里德算法中递归地求：若已有 $b = 0$ ，则 $g c d = a$ ，令 $x = 1, y = 0$ 。
$\qquad$ 否则求出 $x', y'$ 满足 $bx′+(a-a/b∗b)∗y′=gcd(b,a\%b)=gcd(a,b)$ 。
$\qquad$ 于是 $a * y' + b * (x' - a / b * y') = g c d (a, b)$ 。

③辗转相减（除）的其它用处

$\qquad$ 只要一个环定义了带余除法，就可以在上面辗转相减（除）。比如模合数的环、多项式环等等。
$\qquad$ 比如求行列式模一个合数，根据行列式的性质可以把一行的倍数加到另一行上。辗转相减把其中一个位置消成 $0$ 即可。
$\qquad$ 比如求两个多项式的最大公约式，或者说多项式取模，可以不停地把一个的倍数加到另一个上，把其中一个多项式变成 $0$ 。

④类欧几里德算法

$\qquad$ 最基本的模型： $solve(n,A,B,C)=\sum^n_{i=1}⌊\frac{A_i+B}C⌋$ 。
如果 $A \geq C$ ， $ans+=\frac{n(n+1)}2∗(A/C)$ ，然后 $A\%=C$ 。
$\qquad$ 如果 $∣ B ∣ \geq C$ ，讨论下正负号算下，然后把 $B$ 搞到 $[0, C)$ 之间。
$\qquad$ 设 $m=⌊\frac{A_n+B}C⌋$ ，则要算 $\sum^n_{i=1}\sum^m_{j=1}[C_j≤A_i+B]$ 。
$\qquad$ 即 $nm-\sum^m_{j=1}\sum^n_{i=1}[A_i≤C_j-B-1]$ 。
$\qquad$ 即 $n m - s o l v e (m, C, - B - 1, A)$ 。
$\qquad$ 每次 $A$ 对 $C$ 取模后互换位置，复杂度同欧几里德算法，为 $O(log\ C)$ 。

2.中国剩余定理

①中国剩余定理

$\qquad$ 有 $n$ 个方程 $x\equiv a_i\pmod p_i$ ， $p_i$ 两两互质，求 $x$ 。
$\qquad$ 设 $w_i=\prod_{j\ne i}p_j$ ，则答案 $\equiv\sum^n_{i=1}ai×w_i×inv(w_i, p_i)$ 。

②Ex 中国剩余定理

$\qquad$ 如果模数不互质，只要两两合并方程即可。
$\qquad$ 比如 $x\equiv a_1\pmod {p_1},x\equiv a_2\pmod {p_2}$ ，设 $d=gcd(p_1, p_2)$ ，那么必须 $a_1\equiv a_2 \pmod d$ 。然后答案一定可以表示成 $w×d+(a_1\ mod\ d)$ 。
$\qquad$ 用普通中国剩余定理求出 $w$ 即可。
$\qquad w\equiv(a_1/d)\pmod {p1/d},w\equiv(a_2/d)\pmod {p2/d}$ 。
$\qquad$ 每次把两个方程合并成一个模数是它们 $l c m$ 的方程。

③费马小定理

$\qquad$ $p$ 是质数，则 $a^p\equiv a\pmod p$ 。

⑤欧拉定理

$\qquad$ $p > 1$ ， $a, p$ 互质，则 $a^{ϕ(p)}\equiv 1\pmod p$ 。
$\qquad$ 证明： $a, p$ 互质，于是对于缩系内所有 $x$ ， $a x$ 同样构成缩系。即 $\prod_{gcd(x,p)=1}ax\equiv\prod_{gcd(x,p)=1}x$ 。于是 $a^{ϕ(p)}\equiv 1$ 。