Найти обратную величину по модулю

[yurets17]

Здравствуйте! Подскажите пожалуйста ошибку?
Функция должна находить такие х и у, что a * x + b * y = gcd(a, b) = d! но х и у находит неправильные!

#include <stdio.h>
#include <iostream>
#include <cmath>

#define DEBUG
using namespace std;

void extended_euclid(long a, long b, long *x, long *y, long *d)

/* calculates a * *x + b * *y = gcd(a, b) = *d */

{
    
    long q, r, x1, x2, y1, y2;
    
    if (b == 0) {
        
        *d = a, *x = 1, *y = 0;
        
        return;
        
    }
    
    x2 = 1, x1 = 0, y2 = 0, y1 = 1;
    
#ifdef DEBUG
    printf("------------------------------");
    printf("-------------------\n");
    printf("   q    r    x    y    a    b   ");
    printf("x2   x1   y2   y1 \n");
    printf("------------------------------");
    printf("-------------------\n");
#endif
    
    while (b > 0) {
        
        q = a / b, r = a - q * b;
        
        *x = x2 - q * x1, *y = y2 - q * y1;
        
        a = b, b = r;
        
        x2 = x1, x1 = *x, y2 = y1, y1 = *y;
        
#ifdef DEBUG
        printf("%4ld %4ld %4ld %4ld ", q, r, *x, *y);
        printf("%4ld %4ld %4ld %4ld ", a, b, x2, x1);
        printf("%4ld %4ld\n", y2, y1);
#endif
        
    }
    
    *d = a, *x = x2, *y = y2;
    
#ifdef DEBUG
    printf("------------------------------");
    printf("-------------------\n");
#endif
    
}

int main(void)
{
    long a, b, d, x, y;
    printf("Введите значение а: ");
    cin>>a;
    printf("Введите значение b: ");
    cin>>b;
    cout<<"Введите модуль: "; cin>>d;
    extended_euclid(a, b, &x, &y, &d);
    
    printf("x = %ld y = %ld d = %ld\n", x, y, d);
    
    return 0;
}