有必要重新學一下擴展GCD emmmm。

主要是擴展GCD求解線性同余方程$ax≡b (mod p)$。

1.方程有解的充分必要條件：b%gcd（a，p）=0。

　　證明：

• $ax-py=b$
• 由于求解整數解，ax是gcd（a，p）的整數倍，py也是，所以b是gcd（a，p）的整數倍。

2.擴展GCD模板

int exgcd(int a,int b,int &x,int &y) {if(b==0){x=1,y=0;return a;}//注意x，y的賦值。int gcd=exgcd(b,a%b,x,y),t=x;x=y;y=t-a/b*y;return gcd; }

?3.求解線性同余方程：

擴展歐幾里得可以求解形如$ax-py=b$的解。

方程可化為$ax≡b (mod p)$,注意b和p的位置。

令t=gcd（a，p）。方程可化為$\frac {a}{t}x-\frac{p}{t}y=\frac{b}{t}$。exgcd求出$\frac {a}{t}x-\frac{p}{t}y=1$的一組特解x，y。$x*=b/t,y*=b/t$即可求出一組解。

而要求最小整數解，可以發現如果x減p，y加a等式仍然成立，所以最小整數解為（x%p+p）%p.

int GCD(int a,int b){return !b?a:GCD(b,a%b);} int exgcd(int a,int b,int &x,int &y) {if(b==0){x=1,y=0;return a;}int gcd=exgcd(b,a%b,x,y),t=x;x=y;y=t-a/b*y;return gcd; } int fcc(int a,int b,int p) { int x,y,t=GCD(a,p);if(b%t)return -1;int tem=b/t;a/=t,p/=t;exgcd(a,p,x,y);x*=tem,y*=tem;return (x%p+p)%p; }

?

轉載于:https://www.cnblogs.com/Al-Ca/p/11353711.html

總結

以上是生活随笔為你收集整理的模板—扩展GCD*2的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。