数论 —— 最大公约数与最小公倍数

【概念】

1.公约数:有 k 个非零整数 a_{1},a_{2},...,a_{k},若 \exists \: d \in Z^+,s.t. d|a_{1},d|a_{2},...,d|a_{k},则称 d 为 a_{1},a_{2},...,a_{k} 的公约数。

2.最大公约数:公约数中最大的一个数称为最大公约数,记为:GCD(a_1,a_2,...,a_k)

注:

  ① 最大公约数一定是存在的,其最小值为 1。

  ② 当 GCD=1 时,则称这些数是互质的。

  ③ 公约数一定是最大公约数的约数。

3.公倍数:有 k 个非零整数 a_{1},a_{2},...,a_{k},若 \exists \: d \in Z^+,s.t. a_{1}|d,a_{2}|d,...,a_{k}|d,则称 d 为 a_{1},a_{2},...,a_{k} 的公倍数。

4.最小公倍数:公倍数中最小的一个数称为最小公倍数,记为:LCM(a_1,a_2,...,a_k)

注:公倍数一定是最小公倍数的倍数。

【欧几里德算法】

欧几里德算法又称为辗转相除法,用于求两个数的最大公约数,其原理为:GCD(x,y)=GCD(x,y-x),x>y,即:GCD(x,y)=GCD(y,x\:\, mod\:\, y),x>y

1.实现

int GCD(int x,int y){
    if(y==0)
        return x;
    else
        return GCD(y,x%y);
}

2.二进制算法

1)原理

为提高效率,可将欧几里德算法进行优化,通过不断去除因子 2 来降低常数。

  • x=y 时,GCD(x,y)=x=y
  • 当 x、y 均为偶数时,GCD(x,y)=2*GCD(x/2,y/2)
  • 当 x 为偶数,y 为奇数时,GCD(x,y)=GCD(x/2,y)
  • 当 x 为奇数,x 为偶数时,GCD(x,y)=GCD(x,y/2)

2)实现

int GCD(int x,int y){
    if(x==0)
        return y;
    if(y==0)
        return x;

    int i,j;
    for(i=0;!(x&1);i++)//约去2
        x>>=1;
    for(j=0;!(y&1);j++)//约去2
        y>>=1;
    if(i>j)
        i=j;

    while(true){
        if(x>=1;
    }
}

【扩展欧几里德算法】

扩展欧几里德算法是在已知 x、y 时,求解一组 a、b,使得 a*x+b*y=GCD(x,y)

1.过程分析

x>y

① 当 y=0 时,GCD(x,y)=x,此时 a=1,b=0

② 当 x>y>0 时

     设 a_1*x+b_1*y=GCD(x,y)a_2*y+b_2*(x \: \, mod\: \, y)=GCD(y,x \: \, mod\: \, y)

     由于 GCD(x,y)=GCD(y,x \: \, mod\: \, y)

     故 a_1*x+b_1*y=a_2*y+b_2*(x \: \, mod\: \, y)

     即 a_1*x+b_1*y=a_2*y+b_2*(x-[x/y]*y)=a_2*y+b_2*x- b_2*[x/y] *y

     因此 a_1*x+b_1*y=b_2*x+(a_2- [x/y] * b_2)*y

     易得  \left\{\begin{matrix}a1=b2 \\ b_1=a_2- [x/y] * b_2 \end{matrix}\right.

这样就得到了求解 a_1b_1 的方法:a_1b_1 的值基于 a_2b_2

由于 GCD 不断的递归求解,因此一定会在某时 b=0,结束递归,从而得出 a、b 的值。

注: x-[x/y]*y 即为 mod 运算,[x/y] 代表取小于 x/y 的最大整数。

2.实现

int Extended_GCD(int x,int y,int &a,int &b){
    if(y==0){
        a=1;
        b=0;
        return x;
    }

    int gcd=Extended_GCD(b,a%b,y,x);
    b-=a*(x/y);
    return gcd;
}
int main(){
    //形如 ax+by=GCD(x,y)
    int x,y,a,b;
    cin>>x>>y;
    int gcd=Extended_GCD(x,y,a,b);
    cout<<"GCD="<

在求 LCM 时,如果将 LCM 写成 a*b/GCD(a,b),a*b 可能会溢出,正确的方法应该是先除后乘,即:a/GCD(a,b)*b

int GCD(int x,int y){
    return !y?x:GCD(y,x%y);
}
int LCM(int x,int y){
    return x/GCD(x,y)*y;
}

【例题】

  1. 取石子游戏(信息学奥赛一本通-T1218)(欧几里得):点击这里
  2. 青蛙的约会(POJ-1061)(扩展欧几里得):点击这里
    同题:青蛙的约会(洛谷-P1516):点击这里
  3. 最大公约数和最小公倍数问题(洛谷-P1029)(LCM):点击这里
  4. LCM Walk(HDU-5584)(LCM+逆推):点击这里
  5. Happy 2006(POJ-2773)(扩展欧几里得求互素的数):点击这里
  6. LCM Challenge(CF-236C)(三个数的最大的LCM):点击这里
  7. Faraway(HDU-6639)(思维+LCM):点击这里

你可能感兴趣的:(——————数论——————,#,数论——GCD,与,LCM)