Pagini recente » Cod sursa (job #2976733) | Cod sursa (job #914580) | Cod sursa (job #2600165) | Cod sursa (job #813638) | Cod sursa (job #1096057)
/*Invers modular - log(N)
daca gcd(A,N)=1
A*X=1 (mod N ) X e inversul lui A
din gcd extins avem oricare are fi A,N
exista X,Y a.i. AX+BY=gcd(A,B)
deci exista X,y a.i. AX+NY=1
AX+NY=1 (mod N ) <=> AX=1 (mod N)
*/
#include <fstream>
using namespace std;
ifstream f("inversmodular.in");
ofstream g("inversmodular.out");
int N,A,X,inv;
int GCD(int A,int B,int &x,int &y)
{
if(!B)
{
x=1;
y=0;
return A;
}
int x0,y0,D;
D=GCD(B, A % B , x0,y0);
x=y0;
y=x0-(A/B)*y0;
return D;
}
int main()
{
f>>N>>A;
GCD(A,N,X,inv);
while(inv<=0)inv+=N;
g<<inv<<'\n';
//g<<X<<'\n';
return 0;
}