欧几里得算法

最新推荐文章于 2024-08-05 22:05:22 发布

「已注销」

最新推荐文章于 2024-08-05 22:05:22 发布

阅读量961

点赞数

分类专栏：数论****************** soj 文章标签： 1092

本文链接：https://blog.csdn.net/bo_jwolf/article/details/18923515

版权

数论****************** 同时被 2 个专栏收录

117 篇文章 0 订阅

订阅专栏

soj

51 篇文章 0 订阅

订阅专栏

本文介绍了一种用于解决同余方程 x*a=1(mod n) 的扩展欧几里得算法，并通过将该方程转化为 x*a+n*y=1 的形式来求解未知数 x 和 y。通过具体的 C++ 实现代码展示了算法的工作原理。

摘要由CSDN通过智能技术生成

x*a=1(mod n) , 转换成x*a + n * y = 1的形式，进行求解即可

// File Name: 1092.cpp
// Author: bo_jwolf
// Created Time: 2014年02月04日 星期二 15时30分26秒

#include<vector>
#include<list>
#include<map>
#include<set>
#include<deque>
#include<stack>
#include<bitset>
#include<algorithm>
#include<functional>
#include<numeric>
#include<utility>
#include<sstream>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<ctime>

using namespace std;
typedef long long LL;

LL ex_gcd( LL a, LL b, LL &x, LL &y ){
	if( b == 0 ){
		x = 1;
		y = 0;
		return a;
	}
	else{
		 LL temp = ex_gcd( b, a % b, x, y );
		 LL t = x;
		 x = y;
		 y = t - ( a / b ) * y;
		 return temp;
	}
}

int main(){
	LL n, m, a, b, r, s, t, x, y;
	while( cin >> a >> n ){
		if( !a && !n )
			break;
		a = a, b = n;
		r = ex_gcd( a, b, x, y );
		if( x < 0 ){
			x += b;
		}
		//LL res = ( x * a ) % y;
		cout << x << endl;
	}
return 0;
}