Pagini recente » Cod sursa (job #2758906) | Cod sursa (job #787166) | Cod sursa (job #680570) | Cod sursa (job #430738) | Cod sursa (job #2986343)
#include <fstream>
using namespace std;
ifstream fin("inversmodular.in");
ofstream fout("inversmodular.out");
int A, N;
// Oricare ar fi A si N numere intregi exista doua numere X si Y de asemenea intregi astfel incat A * X + N * Y = cmmdc(A, N)
// Cum in problema determinarii inversului modular avem cmmdc(A, N) = 1, exista X si Y astfel incat A * X + N * Y = 1
// X este inversul modular al lui A. Cand X este negativ, adunam N pana devine pozitiv.
void euclidExtins(long long int &x, long long int &y, int a, int b)
{
if (b == 0)
x = 1, y = 0;
else
{
euclidExtins(x, y, b, a % b);
long long int aux = x;
x = y;
y = aux - y * (a / b);
}
}
int main()
{
long long x = 0, y;
fin >> A >> N;
euclidExtins(x, y, A, N);
if (x <= 0)
x += N;
fout << x << "\n";
return 0;
}