1. 程式人生 > >擴充套件的歐幾里德演算法求乘法逆元

擴充套件的歐幾里德演算法求乘法逆元

計算乘法逆元,比如3mod8的乘法逆元為3

是如何用歐幾里得演算法計算的呢???

數對 x,y ,使得 gcd(a,b)=ax+by。
c++語言實現:

#include <iostream>
#include<algorithm>
#include<cstdio>
using namespace std;
int x,y,q;
void ex_Eulid(int a,int  b){
	if(b==0){
		x=1;y=0;q=a;
	}
	else{
		ex_Eulid(b,a%b);
		double temp=x;
		x=y;y=temp-a/b*y;
	}
}
int main() {
	int a,b;
	cin>>a>>b;
	if(a<b)swap(a,b);
	ex_Eulid(a,b);
	printf("%d=(%d)*%d+(%d)*%d\n",q,x,a,y,b);
	return 0;
}
題目實際上就是: 給定 a 和b。
a 要有逆元 , 那麼gcd( a , b ) = 1
假設a的逆元 為x , 那麼就有 a * x mod b = 1
也就是 a * x + b * y = 1
上面的程式中輸入a和b就可以求出對應的x和y。
其中 x 就是 a的逆元。



相關推薦

擴充套件演算法乘法(C語言版)

#include <stdio.h>   int ExtendedEuclid( int f,int d ,int *result);   int main()   {   int x,y,z;   z = 0;   printf("輸入

擴充套件演算法乘法

計算乘法逆元,比如3mod8的乘法逆元為3 是如何用歐幾里得演算法計算的呢??? 數對 x,y ,使得 gcd(a,b)=ax+by。 c++語言實現: #include <iostream&

演算法擴充套件演算法乘法

最近看了一本書《程式設計師》裡面說的一個面試題: 求兩個數的最大公約數: SoEasy的題目看過C 的人都知道怎麼寫這個程式 1.傳統方法:窮舉 #include <math.h>int main(){int m=1970,n=1066,p=0;p=m<n?m:n;for(;p>=1

演算法最大公因子)及擴充套件得(乘法

一、歐幾里得演算法歐幾里得演算法又稱輾轉相除法,是指用於計算兩個正整數a,b的最大公約數。gcd(a,b)=gcd(b,a mod b)。演算法描述:1. 輸入:兩個非負整數a,b,且a≥b。2. 輸出

擴充套件演算法乘法 最小正整數解 直線上的整數點)

參考資料: 本文證明過程來自百度百科和劉汝佳的演算法入門經典。 擴充套件歐幾里得演算法介紹: 前置知識:歐幾里得演算法(其實就是輾轉相除法,用於計算兩個整數a,b的最大公約數。) 歐幾里得演算法: 在開始之前,我們先說明幾個定理: gcd(a,b)=gcd(b,a

擴充套件演算法模版題(+分析+題目)HDU1576 A/B

首先給大家普及一下什麼是擴充套件歐幾里德演算法,它是由歐幾里德演算法演變的,即我們常說的輾轉相除法。 程式碼如下: int gcd(int a,int b){ return b?gcd(b,a%b):a; } 那麼對於不完全為0的非負整數,a,b,gcd(a,b

同餘方程(擴充套件演算法

同餘方程 時間限制: 1 Sec  記憶體限制: 128 MB 題目描述 求關於 x 的同餘方程 ax ≡ 1 (mod b)的最小正整數解。 輸入 輸入只有一行,包含兩個正整數 a, b,用一個空格隔開。 輸出 輸出只有一行,包含

caioj 1153 擴充套件演算法(解不定方程)

模板題 注意exgcd函式要稍微記一下 #include<cstdio> #include<cctype> #include<algorithm> #define

數論-擴充套件演算法

找出一對整數(x,y),使得ax+by=gcd(a,b)。注意,這裡的x和y不一定是正數,也可能是負數或者0.例如,gcd(6,15)=3,6*3-15*1=3,其中,x=3,y=-1.這個方程還有其他解,如x=-2,y=1。 用數學歸納法並不難證明演算法的正確性。此處略去

擴充套件演算法求解線性同餘方程

 歐幾里德演算法  歐幾里德演算法又稱輾轉相除法,用於計算兩個整數a,b的最大公約數。其計算原理依賴於下面的定理:  定理:gcd(a,b) = gcd(b,a mod b)  證明:a可以表示成a = kb + r,則r = a mod b  假設d是a,b的一個公約數,則

演算法的推導與證明 || 擴充套件演算法的解釋說明

序言:      當博主第一次見到歐幾里德演算法時,我是不屑一顧的,由於模板比較好背,所以也沒有仔細研究過其中的數學原理.這段時間突然喜歡上了數學,碰巧同學講了一下基礎數論,就去聽了一聽. 由於博主數學基礎和學習能力都比較差,沒有立即消化其中的知識,於是研究

擴充套件演算法 遞迴和非遞迴實現及證明

關於歐幾里得演算法,貝祖等式,擴充套件歐幾里得演算法,Wikipedia的解釋非常非常詳細了。 另外,看了好多別人優秀的總結,我認為最詳盡的就是ACM之家的總結。 這裡自己再總結一次…實際上就是把別人總結的,我認為有助於自己理解的內容copy過來,再加上

gcd演算法/extgcd擴充套件演算法以及在解不定方程中的應用

這個應該是我在noip前就應該會的東西 ,但是當時也許只是記下了程式碼吧 ,現在有諸多的不理解。後來藉著書和幾篇部落格弄懂了並小證了一下,鑑於網上有些部落格關於這個的寫的真的不好看,所以自己來總結一下,順帶以後也能看。 順帶一提,gcd(a,b)表示a,b的最

#數論# 演算法擴充套件演算法 、費馬小、求解(ing)

歐幾里德求gcd(輾轉相除法): 定理: gcd(a, b) = gcd(b, a % b) 兩個正整數a和b(a>b),它們的最大公約數等於a除以b的餘數c和b之間的最大公約數 證明: a可以表示成a = kb + r,則r = a %

擴充套件演算法Java實現和青蛙相遇

問題銜接: http://poj.org/problem?id=1061 Description 兩隻青蛙在網上相識了,它們聊得很開心,於是覺得很有必要見一面。它們很高興地發現它們住在同一條緯度線上,於是它們約定各自朝西跳,直到碰面為止。可是它們出發之前忘記了一件很重要

演算法最大公約數

歐幾里德演算法又稱輾轉相除法,用於計算兩個整數a,b的最大公約數。 定理:gcd(a,b) = gcd(b,a mod b)。 第一種可以寫成: int Gcd(int a, int b) { while(b != 0) { int r =

HDU 1576 A/B 擴充套件演算法

解決該題的關鍵是: 1、瞭解擴充套件歐幾里德演算法,可以運用其解出gcd(a,b)=ax1+by1中的x1、y1的值 2、由題可得以下內容: n=A%9973,則n=A-A/9973*9973。又A/B=x,則A=Bx。所以Bx-A/9973*9973=n。即Bx-997

演算法(兩數最大公因數)

     兩個整數的最大公因數(gcd)是同時整除兩個大最大整數。即gcd(50,15)=5.      演算法連續計算餘數直到除數為0,最後的非0餘數就是最大公因數。因此若M=1989,N=1590

擴充套件演算法解二元一次不定方程

擴充套件歐幾里德演算法: 已知兩個不完全為 0 的非負整數 a,b,必然存在整數對 x,y ,使它們滿足貝祖等式: 解一定存在,根據數論中的相關定理。下面給出程式碼: int extgcd(int a, int b, int& x, int& y)

輾轉相除法 費馬小定理 尤拉定理 擴充套件演算法簡介

證明:對於(2),p是素數,a是整數且gcd(a,p)=1即他們的最大公約數是1。由於a, 2a, 3a, ……,(p-1)a 模p的餘數都不相同。否則若(i*a) mod p=(j*a) mod p其中 1 =< i < j <= p-1 則p|(j-i)*a, 而gcd(a,p)=1,那