欧几里德及其扩展

#include<bits/stdc++.h>
#define IOS ios::sync_with_stdio(false);cin.tie(0);

using namespace std;
typedef long long ll;

//朴素欧几里德
//求a,b的最大公约数
int gcd(int a, int b) {
	return b == 0 ? a : gcd(b, a % b);
}

//扩展欧几里德
//ax+by=m, m需是gcd(a,b)的倍数才有解
//当ax+by=1时，gcd(a.b)=1;
//d为gcd(a,b),递归回溯结束时的x是ax+by=m的一个特解
//通解 x=x0+b/gcd(a,b)  y=y0-a/gcd(a,b)
//求最小正整数解 (x%gcd(a,b)+gcd(a,b))%gcd(a,b)
//最终a,b的值不变，是原来cin的值；而&d,&x,&y都是新值
void exgcd(int a, int b, int& d, int& x, int& y) {
	if (!b) { d = a, x = 1, y = 0; }
	else { exgcd(b, a % b, d, y, x), y -= x * (a / b); }
}