歐幾里得：

int gcd(int a, int b)
{
    return !b ? a : gcd(b, a%b);
}
int lcm(int a, int b)//最小公倍數
{
    return a / gcd(a, b) * b;//先除後乘避免溢位
}

擴充套件歐幾里得：

存在整數對(x,y)使得ax+by=gcd(a,b)
推導過程：
用遞迴求解擴充套件歐幾里得，設已經求出了下一層遞迴的解，即：ax1+by1=gcd(a,b)的解(x1,y1)
又a%b=a−(a/b)∗b
將(x1,y1)代入到bx1+(a%b)y1=gcd(a,b)中，得
bx1+(a−

int extgcd(int a, int b, int &x, int &y)
{
    int d = a;
    if(b != 0)
    {
        d = extgcd(b, a%b, y, x);
        y -= (a / b) * x;
    }
    else x = 1, y = 0;
    return d;
}

求解不定方程：

若c%gcd(a,b

證明：

首先 a∗x0+a∗k∗b/gcd(a,b)+b∗y0−a∗k∗b/gcd(a,b)=gcd(a,b)
即 a∗(x0+k∗b/gcd(a,b))+b∗(y0−k∗a/gcd(a,b))=gcd(a,b)
通解為x=x0+k∗b/gcd(a,b)，y=y0−k∗a/gcd(a,b)，其中k=...−2,−1,0,1,2...
在所有解中最小的正整數為(x0+b/gcd(a,b))，
所以對於方程a∗x+b∗y=c，最小正整數解（以x

int cal(int a, int b, int c)
{
    int x, y;
    int gcd = extgcd(a, b, x, y);
    if(c % gcd != 0) return -1;
    x *= c/gcd;
    b /= gcd;
    if(b < 0) b = -b;
    int ans = x % b;
    if(ans <= 0) ans += b;
    return ans;
}

同餘方程：

根據上面的內容，我們可以得到：

a∗x≡b(modn)，