辗转相除法（gcd、exgcd）

欧几里得算法

用于求两个数的最大公约数。复杂度为$O(logn)$。

#include <cstdio>
#include <iostream>
using namespace std;

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

int main(void)
{
	int a, b;
	cin >> a >> b;
	cout << gcd(a, b) << endl; 
	return 0;
}

求n个数的最大公约数

$gcd(a,b,c)=gcd(gcd(a,b),c)$

#include <cstdio>
#include <iostream>
using namespace std;
const int N = 1e5 + 5;

int arr[N];

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

int main(void)
{
	int n, a, b, ans;
	cin >> n;
	for (int i = 0; i < n; i++){
		cin >> arr[i];
	}
	ans = arr[0];
	for (int i = 0; i < n - 1; i++){
		ans = gcd(ans, arr[i + 1]);
	}
	cout << ans << endl; 
	return 0;
}

最小公倍数

$lcm(a,b)=(a\ast b)/gcd(a,b)$
$lcm(a,b,c)=lcm(a,b)\ast c/gcd(lcm(a,b),c)$

扩展欧几里德算法

常用于求$ax+by=gcd(a,b)$的一组可行解。返回值为$gcd(a,b)$，同时计算出$x$和$y$。

可以用于求解逆元，根据逆元的定义：$$ab\equiv 1(modp)$$

可以得出不定方程：$$ab-pk=1$$

当$gcd(a,p)=1$时，上述方程一定有解，所以一定存在一个$b$为$a$的乘法逆元。

将$a、p$的值代入方程$ax+by=gcd(a,b)$中的$a、b$，就可以利用扩展欧几里得算法解出这个方程，解出的$x$即为$a$的乘法逆元$b$，$y$为$k$。

这个算法的时间复杂度为$O(logn)$

#include <cstdio>
#include <iostream>
using namespace std;

int exgcd(int a, int b, int& x, int& y)
{
	int d = a;
	if (b != 0){
		d = exgcd(b, a % b, y, x);
		y -= (a / b) * x;
	}
	else{
		x = 1;
		y = 0;
	}
	return d;
}

int main(void)
{
	int a, b, x, y;
	cin >> a >> b;
	cout << exgcd(a, b, x, y) << endl;
	cout << x << " " << y << endl;
	return 0;
}

对于int& x

说明这是传递给这个参数的实参的引用
对这个形参修改 就会对实参修改
是一种传址操作