本文搬运自本人初中博客园博客,若图片加载不出来,可到原文查看:https://www.cnblogs.com/zhangtingxi/p/16531857.html

相关文章:

这篇文章内容之前已经记过一次了,但用的时候又忘了,再记一下

之前的这篇会详细很多

拓展欧几里得复习

ax+by=gcd(a,b)\Large ax+by=\gcd(a,b)

其中 a,ba,b 已知,求 x,yx,y

正常的推导应该都会,拆开后合并同类项最终化为:

ay+b(xaby)=gcd(b,amodb)\Large ay'+b(x'-\frac a by')=\gcd(b,a\bmod b)

所以 x=y,y=abyx=y', y=\dfrac a b y',代码实现时我们是知道当前这一步的 x,yx',y',要推回上一步的 x,yx,y,直接套上面式子

边界条件和打普通gcd一样,若 b=0b=0 就结束,此时 x=1,y=0x=1,y=0

拓欧求逆元

前提:gcd(a,p)=1\gcd(a,p)=1

aa 在模 pp 意义下的逆元,可以表示为求 ax1(modp)ax\equiv 1\pmod p 中的 xx

怎么用拓欧做?回到拓欧式子:

ax+by=gcd(a,b)\Large ax+by=\gcd(a,b)

已知的是 a,pa,p,直接套进去:

ax+py=gcd(a,p)\Large ax+py=\gcd(a,p)

gcd(a,p)\gcd(a,p) 已知为1:

ax+py=1\Large ax+py=1

那模 pp 后必然相同,即:

ax+py1(modp)\Large ax+py\equiv 1\pmod p

约去 pypy

ax1(modb)\Large ax\equiv 1\pmod b

xx 就是答案