NOIp提高组 2012 同余方程————扩展欧几里得定理

题解：本题主要考查扩展欧几里得定理。
简要题意：求$ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$最小正整数解
1.扩展欧几里得定理：$ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$，等价与$a\ast x+m\ast y=1$ ，
因为裴蜀定理$ax+by=z,则gcd(a,b)|z$ 则$gcd(a,b)=1;$
又因为$gcd(a,b)=gcd(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$所以$b\ast x_1+(a\mathrm{m}\mathrm{o}\mathrm{d}b)\ast y_1=gcd;$
又因为$a\mathrm{m}\mathrm{o}\mathrm{d}b=a-(a/b)\ast b;$
所以$b\ast x1+(a-(a/b)\ast b)\ast y1$
$=b\ast x1+a\ast y1–(a/b)\ast b\ast y1$
$=a\ast y_1+b\ast (x1–a/b\ast y1)=gcd$
所以$x=y_1,y=x1–a/b\ast y1;$
因此到最后，如果$x$太小就不断加$b$直到大于等于$0$，太大则一直减$b$，直到最小正整数解。直接求模实现更快。
(其实暴力也可，亲测70分QwQ)
代码如下：

#include<iostream>
#include<cstdio>
using namespace std;
long long a,b,x,xx,y;
void exgcd(long long a,long long b)
{
	if(b==0)
	{
		x=1;y=0;
		return ;
	}
	exgcd(b,a%b);
	xx=x;x=y;y=xx-(a/b)*y;
	return ;
}
int main()
{
	ios::sync_with_stdio(false);
	cin>>a>>b;
	exgcd(a,b);
	cout<<(x+b)%b;
	return 0;
}