欧几里德算法

欧几里德算法又叫辗转相除法，作用是求最大公约数

//其核心方法是如果 a=q*b+r; 那么gcd(a,b) = gcd(b,r);
int gcd(int a, int b) {
    if (b == 0)return a;
    return gcd(b, a % b);
}

扩展欧几里德算法：

根据gcd(a,b) = xa + yb; 可以求出最大公约数gcd(a,b)，x，y；

根据辗转相除法: r_i = q_i+1 * r_i+1 + r_i+2;

　　　　　　　　r_k-1 = q_k * r_k;

把上式改写 r_i+2 = r_i - q_i+1 * r_i+1;

设r_k = x_i * r_i-1 + y_i * r_i;

将r_i 带入以后得到r_k = y_i * r_i-2 + (x_i - q_i-1 * y_i) * r_i-1;

x_i-1 =y_i,   y_i-1 = x_i - q_i-1 * y_i;

eg:已知4=1*12+(-1)×8；就可以得到4=(-1)*32+3*12;

int q,x,y;//q=gcd(a,b);

void exgcd(int a, int b) {
    if (b == 0) {
        x = 1;
        y = 0;
        q = a;
    } else {  //gcd(a,b)=gcd(b,a%b);
        exgcd(b, a%b);
        int temp = x;
        x = y;
        y = temp - a / b * y;
    }
}

 

 

转载于:https://www.cnblogs.com/shijianupc/p/3203190.html