家喻户晓的数学小技巧——欧几里得算法

背景1

• 给定两个数$x$与$y$。
• 求出这两个数的最大公因数。

暴力1

• 直接从1到$min(x,y)$枚举，找出最大的能分别整除它们两个的数即可。
• 时间复杂度：$O(min(x,y))$。

暴力2

• 求出$x$的所有约数，再求出$y$的所有约数。
• 从中查找出公共的数，在所有公共的数中找出最大值即可。
• 时间复杂度$O(\sqrt{max(x,y)})$。

引入——更相减损术

• 更相减损术是出自《九章算术》的一种求最大公约数的算法，它原本是为约分而设计的，但它适用于任何需要求最大公约数的场合。
• 原文：可半者半之，不可半者，副置分母、子之数，以少减多，更相减损，求其等也。以等数约之。
• 说人话：
  • 对于两个数，假如两个数都是偶数，先分别对每一个数除2，直到至少有一个数变成奇数为止。
  • 若这两个数相同，那么最大公约数就是这个数乘上$2^n$，$n$是刚才除2的次数。
  • 否则，设这两个数分别为$x,y$，不妨设$x<y$，则新的$x^{\prime }=y-x$,$y^{\prime }=x$。再重新判断是否相等。
• 举例子：
  • 对于$60$和$100$，求出它们的最大公约数。
  • 首先不断除2，最终得到$15$和$25$。
  • 发现$15\ne 25$，所以$15$变成$25-15=10$，$25$变成$15$。
  • 发现$10\ne 15$，所以$10$变成$15-10=5$，$15$变成$10$。
  • 发现$5\ne 10$，所以$10$变成$10-5=5$，$5$变成$5$。
  • 这时$5=5$，所以最大公约数就是$5\ast 2^2=20$。
• 其实我们发现，之前不用除2也一样可以求出最大公约数。
• 正确性的证明
  • 假若$a$是$x,y$的公因数，$x<y$。则$a$一定也是$y-x,x$的公因数。
  • 这时十分显然的，设$y=a\ast y0,x=a\ast x0$，则$y-x=a\ast (y0-x0)$，$a$也是其因数。
  • 反过来也是一样的，若$a$是$y-x,x$的公因数，$x<y$。则$a$一定也是$x,y$的公因数。
  • 所以$y-x,x$的公因数与$x,y$的公因数是等价的。
  • 对于两个相同的数，它们的最大公因数显然是它们本身，有这个边界条件往上推即可得到原本两个数的最大公因数。
• 但是，假若$x$特别小，$y$特别大，时间复杂度就会退化到$O(y)$，没能完全优化时间复杂度。
• 不过，它给我们提供了一个非常好的思想：

欧几里得算法

• 我们发现，重复的减法其实相当于取模的过程，所以我们直接取模就好了。
• 有$x,y(x<y)$推到$y\mathrm{m}\mathrm{o}\mathrm{d}x,x$。
• 发现这两组数的最大公因数也是等价的。
• 但注意此时的边界条件变为其中一个数为0，则最大公因数为另一个数。
• 时间复杂度：
  • 若$x>\frac{y}{2}$，则$y\mathrm{m}\mathrm{o}\mathrm{d}x<\frac{y}{2}$
  • 若$x<=\frac{y}{2}$，则$y\mathrm{m}\mathrm{o}\mathrm{d}x<x<=\frac{y}{2}$。
  • 即每一次，较大的那个数都会缩减至少一半。
  • 所以极限的时间复杂度是$O(lo{g}_2(max(x,y)))$的。

代码实现

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