Расширенный алгоритм Euclid С++
У меня проблема с расширенным алгоритмом Евклида. (ax + by = gcd (a, b)) Я пытаюсь определить как GCD, так и x и y. GCD не проблема, но с использованием метода цикла что-то идет не так с x и y. Обычно одно число появляется как 0, а другое - аномально большое отрицательное число. Код следует:
#include <iostream>
using namespace std;
main ()
{
int a,b,q,x,lastx,y,lasty,temp,temp1,temp2,temp3;
cout << "Please input a" << endl;
cin >> a;
cout << "Please input b" << endl;
cin >> b;
if (b>a) {//we switch them
temp=a; a=b; b=temp;
}
//begin function
x=0;
y=1;
lastx=1;
lasty=0;
while (b!=0) {
q= a/b;
temp1= a%b;
a=b;
b=temp1;
temp2=x-q*x;
x=lastx-q*x;
lastx=temp2;
temp3=y-q*y;
y=lasty-q*y;
lasty=temp3;
}
cout << "gcd" << a << endl;
cout << "x=" << lastx << endl;
cout << "y=" << lasty << endl;
return 0;
}
Ответы
Ответ 1
Два из ваших назначений ошибочны:
temp2 = x;
x=lastx-q*x;
lastx = temp2;
temp3 = y;
y = lasty-q*y;
lasty=temp3;
Пример вывода с приведенными выше исправлениями:
Please input a
54
Please input b
24
gcd6
x=1
y=-2
Ответ 2
Хотя вопрос задан давно, но ответ поможет кому-то, кто нашел реализацию расширенного евклидова алгоритма на С++.
Рекурсивная реализация С++:
int xGCD(int a, int b, int &x, int &y) {
if(b == 0) {
x = 1;
y = 0;
return a;
}
int x1, y1, gcd = xGCD(b, a % b, x1, y1);
x = y1;
y = x1 - (a / b) * y1;
return gcd;
}
Пример с кодом:
#include <iostream>
int main()
{
int a = 99, b = 78, x, y, gcd;
if(a < b) std::swap(a, b);
gcd = xGCD(a, b, x, y);
std::cout << "GCD: " << gcd << ", x = " << x << ", y = " << y << std::endl;
return 0;
}
Ввод:
a = 99, b = 78
Вывод:
GCD: 3, x = -11, y = 14