金华集训Day7——宇

Contents)

一、质因子分解

1.素性测试

①两种常见的素性测试

\qquad 试除法，配合质数筛法可以做到$O(\sqrt{n})$的复杂度。
\qquad $Miller-Rabin$素性测试，可以做到$O(logn)$的复杂度，属于$RP$算法。
\qquad 其它还有印度人的$AKS$算法可以做到正确率$100\%$的$O(logn)$复杂度。
\qquad 绝大多数情况$Miller-Rabin$都够优秀了。
\qquad 可以查阅$wiki$以确定$Miller-Rabin$应该试到多少。

②Miller-Rabin

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

2.质因子分解

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

\qquad 试除法，复杂度$O(\sqrt{n})$，太慢了。
\qquad $Pollar{d}^{\prime }sRho$，期望复杂度$O(\sqrt[4]{n}logn)$，有名启发式分解。

②为什么叫启发式分解？

\qquad $n$是素数的时候用$Miller-Rabin$素性测试。
\qquad $n$不是素数的时候，复杂度只和最小的质因子有关。所以称为启发式。

③Pollard’s Rho

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

④如何判已经进入循环？

\qquad 可以证明$x\leftarrow x^2+c$，形成的一定是一个$p$形结构。
\qquad 每次当 i 为 2 的幂次的时候就令$y\leftarrow x_i$，如果某时刻$x_i=y$了则说明已经在环上绕了一圈了。
\qquad 即：看$x_{(2,4]}$是否$=x_2$，看$x_{(4,8]}$是否$=x4$，看$x_{(8,16]}$是否$=x_8\dots \dots$
\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$满足$ax+by=gcd(a,b)$。
\qquad 在欧几里德算法中递归地求：若已有$b=0$，则$gcd=a$，令$x=1,y=0$。
\qquad 否则求出$x\prime ,y\prime$满足$bx\prime +(a-a/b\ast b)\ast y\prime =gcd(b,a\%b)=gcd(a,b)$。
\qquad 于是$a\ast y\prime +b\ast (x\prime -a/b\ast y\prime )=gcd(a,b)$。

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

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

④类欧几里德算法

\qquad 最基本的模型：$solve(n,A,B,C)={\sum }_{i=1}^n\lfloor \frac{A_i+B}{C}\rfloor$。
如果$A\ge C$，$ans+=\frac{n(n+1)}{2}\ast (A/C)$，然后$A\%=C$。
\qquad 如果$|B|\ge C$，讨论下正负号算下，然后把$B$搞到$[0,C)$之间。
\qquad 设$m=\lfloor \frac{A_n+B}{C}\rfloor$，则要算${\sum }_{i=1}^n{\sum }_{j=1}^m[C_j\le A_i+B]$。
\qquad 即$nm-{\sum }_{j=1}^m{\sum }_{i=1}^n[A_i\le C_j-B-1]$。
\qquad 即$nm-solve(m,C,-B-1,A)$。
\qquad 每次$A$对$C$取模后互换位置，复杂度同欧几里德算法，为$O(logC)$。

2.中国剩余定理

①中国剩余定理

\qquad 有$n$个方程$x\equiv a_i(\mathrm{m}\mathrm{o}\mathrm{d}p{)}_i$，$p_i$两两互质，求$x$。
\qquad 设$w_i={\prod }_{j̸=i}{p}_j$，则答案$\equiv {\sum }_{i=1}^{n}ai\times w_i\times inv(w_i,p_i)$。

②Ex 中国剩余定理

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

③费马小定理

\qquad $p$是质数，则$a^p\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}p)$。

⑤欧拉定理

\qquad $p>1$，$a,p$互质，则$a^{\varphi (p)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$。
\qquad 证明：$a,p$互质，于是对于缩系内所有$x$，$ax$ 同样构成缩系。即${\prod }_{gcd(x,p)=1}ax\equiv {\prod }_{gcd(x,p)=1}x$。于是$a^{\varphi (p)}\equiv 1$。

⑥扩展欧拉定理

\qquad

离散对数问题

原根相关

二次剩余

一些例题

数论函数

定义

计算积性函数

一些例题