最大公约数简谈

根据教师课件，写成本文。未经允许，禁止转载。

$\gcd$，即最大公约数，指的是几个整数中公有的约数中最大的一个。

问题

已知 $a,b$，求 $\gcd (a,b)$。

方法一：枚举法

从 $\min (a,b)$ 到 $1$ 枚举，找到第一个 $x$ 符合题意，就退出循环。

时间复杂度 $\mathcal{O}(\min (a,b))$。

方法二：分解质因数

令 $a=p_1^{x_1}{p}_2^{x_2}\dots p_n^{x_n}$，$b=p_1^{y_1}{p}_2^{y_2}\dots p_n^{y_n}$，其中$x_i,y_i\ge 0$ 且不同时为 0。

则 $\gcd (a,b)=p_1^{\min (x_1,y_1)}{p}_2^{\min (x_2,y_2)}\dots p_3^{\min (x_3,y_3)}$。

时间复杂度 $\mathcal{O}(\sqrt{\min (a,b)})$。

void Gcd()
{
	for(int x=2;x*x<=min(a,b);x++)
	{
		while (a % x==0 && b % x==0) {a/=x;b/=x;ans*=x;}
		while (a % x==0)a/=x;
    	while (b % x==0)b/=x;
	} 
	if (a % b==0)ans*=b;
    else if (b % a==0)ans*=a;
    printf("%d",ans);
}

方法三：辗转相除（欧几里得算法）

定理：$\gcd (a,b)=\gcd (b,a\%b)$。

证明：

设 $\gcd (a,b)=p$，则有 $a=a^{\prime }\ast p,b=b^{\prime }\ast p,\gcd (a^{\prime },b^{\prime })=1$。

$a\%b=a-\lfloor \frac{a}{b}\rfloor \ast b=a^{\prime }\ast p-\lfloor \frac{a}{b}\rfloor \ast b^{\prime }\ast p=p\ast (a^{\prime }-\lfloor \frac{a}{b}\rfloor \ast b^{\prime })$。

$\gcd (b,a\%b)=\gcd (b^{\prime }\ast p,p\ast (a^{\prime }-\lfloor \frac{a}{b}\rfloor \ast b^{\prime }))=p\ast \gcd (b^{\prime },a^{\prime }-\lfloor \frac{a}{b}\rfloor \ast b^{\prime })$。

现证明 $\gcd (b^{\prime },a^{\prime }-\lfloor \frac{a}{b}\rfloor \ast b^{\prime })=1$，使用反证法。假设 $gcd(b^{\prime },a^{\prime }-\lfloor \frac{a}{b}\rfloor \ast b^{\prime })=t(t>1)$。

$b^{\prime }=b^{\prime \prime }\ast t$，$a^{\prime }-\lfloor \frac{a}{b}\rfloor \ast b^{\prime }=c^{\prime }\ast t$。

$a^{\prime }=c^{\prime }\ast t+\lfloor \frac{a}{b}\rfloor \ast b^{\prime }=c^{\prime }\ast t+\lfloor \frac{a}{b}\rfloor \ast b^{\prime \prime }\ast t=t\ast (c^{\prime }+\lfloor \frac{a}{b}\rfloor \ast b^{\prime \prime })$。

则 $\gcd (a^{\prime },b^{\prime })\ge p$，与 $\gcd (a^{\prime },b^{\prime })=1$ 矛盾，因此 $\gcd (b^{\prime },a^{\prime }-\lfloor \frac{a}{b}\rfloor \ast b^{\prime })=1$。

$\gcd (b,a\%b)=p=\gcd (a,b)$。

---

时间复杂度 $\mathcal{O}(\log (\max (a,b)))$。

分析：

1. 设 $a>b$。

   若 $a>2b$，则 $b\le \frac{a}{2}$，规模减小一半。反之 $a<2b$，则 $a\%b<\frac{a}{2}$。

   因此时间复杂度是 $\log$ 级别。

2. 斐波那契分析：传送门。

---

int Gcd(int a,int b)
{
    if (b==0) return a;
    else return Gcd(b,a % b);
}

方法四：二进制法

当 $a<b$ 时， $\gcd (a,b)=\gcd (b,a)$。

当 $a=b$ 时 $\gcd (a,b)=a$。

当 $a,b$ 同为偶数时，$\gcd (a,b)=2\ast \gcd (\frac{a}{2},\frac{b}{2})$。

当 $a$ 为偶数，$b$ 为奇数时，$\gcd (a,b)=\gcd (\frac{a}{2},b)$。

当 $a$ 为奇数，$b$ 为偶数时，$\gcd (a,b)=\gcd (a,\frac{b}{2})$。

当 $a,b$ 为奇数时，$\gcd (a,b)=\gcd (a-b,b)$。

注：此法适合高精度求最大公约数。

int Gcd(int m,int n)
{
    if (m==n) return m;
    if (m<n) return Gcd(n,m);
    if (m & 1==0) return (n & 1==0)? 2*Gcd(m/2,n/2):Gcd(m/2,n);
    return (n & 1==0)? Gcd(m,n/2): Gcd(n,m-n);
}