乘法逆元應(yīng)用在組合數(shù)學(xué)取模問(wèn)題中，這里給出的實(shí)現(xiàn)不見(jiàn)得好用

給出拓展GCD算法：

擴(kuò)展歐幾里得算法是指對(duì)于兩個(gè)數(shù)a,b 一定能找到x,y（均為整數(shù)，但不滿足一定是正數(shù)） 滿足x*a+y*b=gcd(a,b) gcd（x,y)是指x 與 y的最大公約數(shù)

有啥用呢？求解形如 a*x +b*y = c 的通解

然后我們先介紹同余方程，再介紹乘法逆元

同余方程 a≡b(mod m) 等價(jià)于小學(xué)的運(yùn)算式 b÷m 余數(shù)為a 也就是a mod m=b

其實(shí)介紹這個(gè)就是看怎么把≡拿掉

乘法逆元 ax ≡ 1 (mod m) 我們稱 x 是 a 關(guān)于 m 的乘法逆元 可以等價(jià)于這樣的表達(dá)式： a*x + m*y = 1

當(dāng)滿足這個(gè)式子的時(shí)候：a*x + b*y = c 有解的充要條件： c % gcd(a , b) == 0

一般，我們能夠找到無(wú)數(shù)組解滿足條件，但是一般是讓你求解出最小的那組解

我們求解出來(lái)了一個(gè)特殊的解 x0 ，我們用 x0 % m其實(shí)就得到了最小的解了

1 #include<cstdio> 2 using namespace std; 3 inline long long read() 4 { 5 long long x=0,f=1;char ch=getchar(); 6 while(ch<'0'||ch>'9') {if(ch=='-')f=-1;ch=getchar();} 7 while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} 8 return x*f; 9 } 10 int a,b; 11 void exgcd(int a,int b,int &x,int &y) 12 { 13 if(b==0) {x=1;y=0;return;} 14 exgcd(b,a%b,x,y); 15 int t=x;x=y;y=t-a/b*y; 16 } 17 //ax ≡ 1 (mod b) 18 //-> a*x + b*y = 1 19 //->求出x和y后讓x%b就是最小解了 20 int main() 21 { 22 a=read();b=read(); 23 int x,y; 24 exgcd(a,b,x,y); 25 x=(x%b+b)%b; 26 printf("%d",x); 27 return 0; 28 }

?

轉(zhuǎn)載于:https://www.cnblogs.com/aininot260/p/9480161.html

創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎(jiǎng)勵(lì)來(lái)咯，堅(jiān)持創(chuàng)作打卡瓜分現(xiàn)金大獎(jiǎng)

總結(jié)

以上是生活随笔為你收集整理的数学：乘法逆元-拓展GCD的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。