扩展欧几里得算法求逆元

//写在前面
//当有ax=1(mod n)---1 时,x称做a的乘法逆元
//一个数有逆元的充分必要条件是gcd(a,n)=1
// 1式等价于 ax+ny=1 ---2
//2式可用扩展欧几里得算法求出x0的值
//其中最小逆元 x=(x0%n+n)%n 
//求出逆元后 (a/b)%n=(a*inv(a))%n 其中a*inv(a)=1(mod n) 
#include 
#define LL long long
using namespace std;
int gcd(int a,int b)//欧几里得算法判断是否有逆元 
{
	if(b==0)return a;
	else gcd(b,a%b);
}
void Ex_gcd(int a,int b,int &x,int &y)//扩展欧几里得算法 
{
	if(b==0)
	{
		x=1;
		y=0;
		return;
	}
	int x1,y1;
	Ex_gcd(b,a%b,x1,y1);
	x=y1;
	y=x1-(a/b)*y1;
}
int main()
{
	LL a,b;cin>>a>>b;
	int x,y;
	Ex_gcd(a,b,x,y);
    
    x=(x%b+b)%b;//最小的逆元 
    cout<

 

你可能感兴趣的:(算法模板,数论)