【NOIP2012 提高组 day2】同余方程

分析

数论
原式可以转化为ax+by=1
用拓展欧几里得算法求这个方程中x的最小值

代码如下

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll a,b,x,y;
ll read(){
	ll sum=0,f=1;
	char ch=getchar();
	while(ch>'9'||ch<'0')
	{
		if(ch=='-')f=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9')
	{
		sum=(sum<<3)+(sum<<1)+ch-'0';
		ch=getchar();
	}
	return sum*f;
}
void exgcd(ll a,ll b,ll&x,ll&y)
{
	if(b==0){
        x=1;y=0;
		return ;
	}
	exgcd(b,a%b,x,y);
	ll t;
	t=x; x=y; y=t-(a/b)*y;
}
int main(){
//	freopen("mod.in","r",stdin);
//	freopen("mod.out","w",stdout);
	a=read();
	b=read();;
    exgcd(a,b,x,y);
    cout<<(x%b+b)%b;
	return 0;
}