歐幾里得算法：

求a，b的最大公約數(shù)
gcd（a,b）= gcd(b,a%b)

ll gcd(ll a,ll b) {return b==0?a:gcd(b,a%b); }

擴(kuò)展歐幾里得算法：

如果a，b是整數(shù)，一定存在x和y使得ax+by=gcd（a，b）
也就是ax+by=m的話,m一定是gcd（a，b）的倍數(shù)
假設(shè)當(dāng)前我們在求的時a和b的最大公約數(shù)，而我們已經(jīng)求出了下一個狀態(tài)：b和a%b的最大公因數(shù)，并且求出了一組x1和y1使得 ： b * x1+(a%b) * y1 = gcd
結(jié)論： x = y1 , y = x1 – a / b * y1
得到的是特解，如何由特解推出其他的整數(shù)解：
對于關(guān)于x，y的方程a * x + b * y = g 來說，讓x增加b/g，讓y減少a/g，等式兩邊還相等。
對于ax+by=c的一組解（x0,y0）,則它的任意整數(shù)解都可以寫成（x0+k * b1,y0+k * a1）其中 b1 = b / gcd(a,b), a1 = a/gcd(a,b).
我們就可以用 (x0 % b1 + b1 ) % b1得到它的最小正整數(shù)解了
此處x0= x * c / gcd(a,b)

#include<iostream> #include<cstdio> #include<cmath>using namespace std;int exgcd(int a,int b,int &x,int &y)//擴(kuò)展歐幾里得算法 {if(b==0){x=1;y=0;return a; //到達(dá)遞歸邊界開始向上一層返回}int gcd=exgcd(b,a%b,x,y);int y1=y; //把x y變成上一層的int x1=x;y=x1-(a/b)*y1;x=y1;return gcd; //得到a b的最大公因數(shù) }

總結(jié)

以上是生活随笔為你收集整理的数论 欧几里得与扩展欧几里得的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

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