#include<cstdio>
using namespace std;
int x,y;
int gcd(int a,int b,int &x,int &y){
    if(!b){
        x=1,y=0;
        return a;
    }
    int ret=gcd(b,a%b,x,y);
    int t=x;
    x=y;
    y=t-a/b*y;
    return ret;
} 
int a,b;
int main(){
    scanf("%d%d",&a,&b);
    gcd(a,b,x,y
 
);
    printf("%d",(x+b)%b);


}
複製程式碼
歐幾里得gcd(a,b)=gcd(b,a%b)=......

因為ax≡1（mod b） -》ax%b=1%b=1

所以ax+by=1,因為y是整數所以加個by就相當於%b(因為%b的本質是+上y個b)，所以兩個式子等價

由擴充套件歐幾里得得 ax+by=gcd(a，b)=gcd(b,a%b)=bx+(a%b=a-a/b*b)*y=....(x,0)=x。(x是最大公約數,/是整除)

最後bx=b,(a%b)*y=0

所以得 x=1，y=0

因此遞迴就行了

x2=y1,y2=x1-a/b*y1(原x-整除的數*a/b=模數)

最後得到的是無數解中的一種，有可能是負的，只需要加上模數就行了