乘法逆元

给出2个数M和N(M < N)，且M与N互质，找出一个数K满足0 < K < N且K * M % N = 1，如果有多个满足条件的，输出最小的。

Input

输入2个数M, N中间用空格分隔（1 <= M < N <= 10^9)

Output

输出一个数K，满足0 < K < N且K * M % N = 1，如果有多个满足条件的，输出最小的。

Sample Input

2 3

Sample Output

2

逆元：设c是b的逆元，则有b*c≡1(mod m)；

求逆元：扩展欧几里得；

算法作用：求解a,b 最大公约数m 和a*x+b*y=m 的一个解；

证明：   a*x1+b*y1=gcd(a,b),b*x2 +(a%b)* y2=gcd(b,a%b)
             由辗转相除法可得到gcd(a,b)= gcd(b,a%b)
             从而得到a*x1+b*y1=b*x2 +(a%b)* y2
             即a*x1+b*y1=b*x2 +(a -(a/b)*b)* y2=a*y2+b*(x2 -(a/b)* y2)
             所以x1=y2 ,y1=x2 -(a/b)* y2.

代码1：提交c++，不然好像会编译错误。

#include<stdio.h>
int exgcd(int a,int b,int &x,int &y)
{
	if(!b)
	{
		x=1;
		y=0;
		return a;
	}
	int r=exgcd(b,a%b,x,y);
	int t=x;
	x=y;
	y=t-(a/b)*y;
	return r;
}
int main()
{
	int a,b,x,y;
	scanf("%d %d",&a,&b);
	int c=exgcd(a,b,x,y);
	printf("%d",(x+b)%b);
	return 0;
}

代码2：

#include<stdio.h>
typedef long long ll;
ll extend_gcd(ll a,ll b,ll &x,ll &y)
{
	if(b==0)
	{
		x=1,y=0;
		return a;
	}
	else
	{
		ll r=extend_gcd(b,a%b,y,x);
		y-=x*(a/b);
		return r;
    }	
}
ll inv(ll a,ll n)
{
	ll x,y;
	extend_gcd(a,n,x,y);
	x=(x%n+n)%n;
	return x;
}
int main()
{
	ll m,n;
	scanf("%lld %lld",&m,&n);
	printf("%lld\n",inv(m,n));
	return 0;
}

代码大致一样，格式不同，相互对比一下，以供学习，可以当模板记。