首页 > 其他分享> 文章详细

正睿金华Day6数论&杂题选讲

2019-08-12 16:44:05 阅读：374 来源： 互联网

标签：right gcd Day6 bmod 选讲 VeryThinSpace 正睿 equiv mod

<前言>

好久没写博客了，从Day5开始，那么我就从Day6开始补吧。

等等让我找找day6讲什么的、。。。

偶，是任轩笛讲的，上午讲数论和数论函数，下午杂题选讲。

<正文>

质因数

一开始讲的是质因数的素性测试、质因子分解之类的，听着还挺正常。讲到线性筛的时候感觉还行，就去上了个厕所回来。

素性测试：
- 试除法，配合质数筛法可以 $O(\sqrt n)$ O(n)解决
  
  这个只要筛出 $\sqrt n$ n之内的素数。然后一直试除，能除就除，只要能出就不是素数
- Miller-Rabin素性测试， $O(log\ n)$ O(log n)完成但可能错误。
  
  这个我没听啊，但大致讲一下吧。
  - 基本原理是费马小定理：若 p 是质数，a, p 互质，则 $a^{p-1} \equiv 1(\bmod\ p)$ ap−1≡1(mod p)
  - 于是对于某个 p，若能找到与它互质的 a 使得 $a^{p-1} \neq 1(mod\ p)$ ap−1̸=1(mod p)，则p不是质数。
  - 然而有一些合数 p，满足所有与它互质的 a 都有 $a^{p-1} \equiv 1(\bmod\ p)$ ap−1≡1(mod p)，这种数称为 Carmichael 数（如 561 = 3 ∗ 11 ∗ 17），这些数是用上面的方法检验不出来的。
  - 所以还需要奇素数判定。对于奇素数 p， $\left(a^{\frac{p-1}{2}}+1\right)\left(a^{\frac{p-1}{2}}-1\right) \equiv 0\ (\bmod \ p)$ (a2p−1+1)(a2p−1−1)≡0 (mod p)，由于 $F_p$ Fp是整环，所以 $a^{\frac{p-1}{2}} \equiv 1$ a2p−1≡1 或 $p-1$ p−1.
  - 如果 $\frac{p-1}{2}$ 2p−1还是偶数则可以继续往下检验。
- 印度人的AKS就不说了，是保证正确的log n算法:)
质因数分解：
- 试除法：能除就除，不能除就加除数，直到大于原数。f复杂度就算用线性筛优化寻找素数，也是 $O(\sqrt n)$ O(n)，太慢了。
- Pollard’s Rho算法，期望复杂度 $O(\sqrt[4]{n} \log n)$ O(4nlogn)
  - 假如要分解一个数 n，首先进行素性测试，是素数直接返回。
  - 否则就要生成一些随机的 $x_i$ xi，去求 $gcd(|x_i − x_j|, n)$ gcd(∣xi−xj∣,n)，如果这个 $∈ (1, n)$ ∈(1,n) 则找到了 $n$ n的一个因子，递归分解。
  - 一个挺靠谱的随机方法就是 $x \leftarrow x^{2}+c$ x←x2+c,c为随机数。
  - 这样随机出来的 x 可能会进入循环，假如进入循环了我们还没找到因子，就重新随个 x 和 c，重新做。
  - 如何判定是否进入环中：
    - 可以证明 $x ← x^2 + c$ x←x2+c，形成的一定是一个 ρ 形结构,从一头进入循环就出不来了。我们可以采取标记的方法
    - 每次当 i 为 2 的幂次的时候就令 $y ← x_i$ y←xi，如果某时刻 $x_i = y$ xi=y
      了则说明已经在环上绕了一圈了。
      即：看 x(2,4] 是否 = x2，看 x(4,8] 是否 = x4，看 x(8,16] 是否
      = x8……
      
      这样“浪费”的步数仅仅是 O( 环长 ) 级别的。

数论

老师的ppt足够详细，但是我还是有不明白。

但就着重讲几个我会的吧：

欧几里得算法

在求两数gcd得时候，我们用的多是辗转相除法或者辗转相减法。欧几里得算法就是辗转相除法了。

如果 $x|a,x|b$ x∣a,x∣b,那么 $x|(a-b)(a>b)$ x∣(a−b)(a>b),所以gcd(a,b)=gcd(a,a%b).这样不断转换直到a%b变成零，此时a就是gcd了。、

扩展欧几里得：

已知 a, b，求出 x, y 满足 ax + by = gcd(a, b)。

在欧几里德算法中递归地求，若已有 b = 0，则 gcd = a，令x = 1, y = 0。

否则求出x‘, y′满足
$b x^{\prime}+(a-a / b * b) * y=g c d(b, a \% b)=g c d(a, b)$ bx′+(a−a/b∗b)∗y=gcd(b,a%b)=gcd(a,b)

于是：
$a * y+b *(x-a / b * y)=g c d(a, b)$ a∗y+b∗(x−a/b∗y)=gcd(a,b)

有 $，x1=y2; y1=x2-a/b*y2;$ ，x1=y2;y1=x2−a/b∗y2;然后对x1，y1做，递归解决，回溯时实行逆运算，最后找到最小解。

其实扩欧我们老早学过了，不过忘得差不多了，这次讲的时候听不进去，已经被任轩笛的骚操作秀到了。

类欧几里得：

额，抱歉，这个是真的挂了。在这个暑假之前我根本不知道有这个玩意。先是在hl集训中出现，然后在这里又讲了，可惜我完全看不懂公式啊。

$\operatorname{solve}(n, A, B, C)=\sum_{i=1}^{n}\left\lfloor\frac{A i+B}{C}\right\rfloor$ solve(n,A,B,C)=i=1∑n⌊CAi+B⌋

这是我全程看得懂的一个公式（也就是定义）

中国剩余定理

有 n 个方程 x ≡ ai (mod pi)，pi 两两互质，求 x。

mb题养猪

设 $w_{i}=\prod_{j \neq i} p_{j}$ wi=∏j̸=ipj，则答案 $\equiv \sum_{i=1}^{n} a_{i} * w_{i} * i n v\left(w_{i}, p_{i}\right)$ ≡∑i=1nai∗wi∗inv(wi,pi)

inv为逆元

已知 $x \equiv a_{1}\left(\bmod p_{1}\right), x \equiv a_{2}\left(\bmod p_{2}\right)$ x≡a1(modp1),x≡a2(modp2)，若 $d=\operatorname{gcd}\left(p_{1}, p_{2}\right)$ d=gcd(p1,p2)
则 $a_{1} \equiv a_{2}(\bmod d)$ a1≡a2(modd)成立

所以答案可以表示为 $w * d+\left(a_{1} \bmod d\right)$ w∗d+(a1modd)

求得：

$w \equiv\left(a_{1} / d\right)\left(\bmod p_{1} / d\right), w \equiv\left(a_{2} / d\right)\left(\bmod p_{2} / d\right)$ w≡(a1/d)(modp1/d),w≡(a2/d)(modp2/d)

每次把两个方程合并成一个模数是它们 lcm 的方程

费马小定理

时间不多长话短说：

费马小定理：p 是质数，则 $a^{p} \equiv a(\bmod p)$ ap≡a(modp)

欧拉定理：

p>1,a,p互质，则 $a^{\phi(p)} \equiv 1(\bmod p)$ aϕ(p)≡1(modp)

扩展欧拉定理

$p>1, \quad n \geq \phi(p)$ p>1,n≥ϕ(p)情况下，存在

$a^{n} \equiv a^{n \bmod \phi(p)+\phi(p)}(\bmod p)$ an≡anmodϕ(p)+ϕ(p)(modp)

还有一个数论函数，不写了不写了，反正横竖也全是题目。

算了，不能指望我听懂这些，果然听数论就是个错误的选择吗。

更可恶的是，任轩笛

标签：right,gcd,Day6,bmod,选讲,VeryThinSpace,正睿,equiv,mod
来源： https://blog.csdn.net/qq_40900472/article/details/99184087

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9