最 大 公 约 数 [Greatest Common Divisor]

$最大公约数$

$欧几里得算法$

$不妨设a>b$
$1.$ $能整除$
$b是a的约数，那么b就是二者的最大公约数$

$2.$ $不能整除$
$①$ $gcd(a,b)=>gcd(b,a$ $mod$ $b)$

$设a=bk+c(c<b)，显然有c=a$ $mod$ $b。设d$ $|$ $a，d$ $|$ $b，则c=a-bk,\frac{c}{d}=\frac{a}{d}-\frac{b}{d}k，\frac{c}{d}为整数，即d$ $|$ $c，\therefore gcd(a,b)=>gcd(b,a$ $mod$ $b)$

$②$ $gcd(b,a$ $mod$ $b)=>gcd(a,b)$

$设d$ $|$ $b，d$ $|$ $(a$ $mod$ $b)$
$=>\frac{amodb}{d}=\frac{a}{d}-\frac{b}{d}k=>\frac{a}{d}=\frac{b}{d}k+\frac{amodb}{d}，\frac{a}{d}为整数，即d$ $|$ $a，\therefore gcd(b,a$ $mod$ $b)=>gcd(a,b)$

$综上，gcd(a,b)=gcd(b,a$ $mod$ $b)$

实现
$①递归求法$

//version 1
int gcd(int a, int b){
	if(!b) return a;//!b<=>b == 0
	return gcd(b, a % b);
}

//version 2
int gcd(int a, int b){return (!b) ? a : gcd(b, a % b);}//!b<=>b == 0

$②迭代求法$

//version 1
int gcd(int a, int b){
	while(b){//b<=>b != 0
		int tmp = a;
		a = b;
		b = tmp % b;
	}
	return a;
}

异或使人疑惑

//version 2
int gcd(int a, int b){
	while(b){//b<=>b != 0
		a %= b;
		b ^= a;
		a ^= b;
		b ^= a;
	}
	return a;
}

//version 3
int gcd(int a, int b){
    while(b ^= a ^= b ^= a %= b);
    return a;
}

$c++14$

__gcd(a, b);

$c++17$

std::gcd(a, b);

$时间复杂度为O(log$ $max(a,b))$

• $a\ge b，gcd(a,b)=gcd(b,a$ $mod$ $b)，a$ $mod$ $b会让a至少减半$
  $证明：$
  $①当b\le \frac{a}{2}时，a$ $mod$ $b的最大值为b-1<\frac{a}{2}$
  $②当b>\frac{a}{2}时，有a$ $mod$ $b=a-\lfloor \frac{a}{b}\rfloor \times b$
  $\because b>\frac{a}{2}$
  $\therefore \frac{a}{b}<2$
  $\therefore \lfloor \frac{a}{b}\rfloor =1$
  $\therefore a$ $mod$ $b=a-b$
  $\because b>\frac{a}{2}$
  $\therefore a-b<\frac{a}{2}$
  $即a$ $mod$ $b<\frac{a}{2}$
  $综上，a$ $mod$ $b<\frac{a}{2}$
  $证毕.$
  $时间复杂度为O(log$ $a)=O(n)$

• $a<b，gcd(a,b)=gcd(b,a)=gcd(a,b$ $mod$ $a)$
  $时间复杂度为O(log$ $b)=O(n)$

$欧几里得算法求斐波那契数列相邻两项的最大公约数，会让该算法达到最坏复杂度$

$F_n=F_{n-1}+F_{n-2}$
$F_n$ $mod$ $F_{n-1}=(F_{n-1}+F_{n-2})$ $mod$ $F_{n-1}$
$\because F_{n-1}$ $mod$ $F_{n-1}=0$
$\therefore F_n$ $mod$ $F_{n-1}=F_{n-2}$ $mod$ $F_{n-1}=F_{n-2}$

$gcd(F_n,F_{n-1})=gcd(F_{n-1},F_{n-2})=\dots =gcd(F_2,F_1)$

$迭代次数为n次，时间复杂度为O(n)$