Macaulay2 Engine
Loading...
Searching...
No Matches

◆ gcd_extended()

long gcd_extended ( long a,
long b,
long & u,
long & v )

Definition at line 14 of file dpoly.cpp.

15{
16 long g;
17 long q;
18 long u1, v1, g1;
19 long utemp, vtemp, gtemp;
20
21 g1 = b;
22 u1 = 0;
23 v1 = 1;
24 g = a;
25 u = 1;
26 v = 0;
27 while (g1 != 0)
28 {
29 q = g / g1;
30 gtemp = g - q * g1;
31 utemp = u - q * u1;
32 vtemp = v - q * v1;
33 g = g1;
34 u = u1;
35 v = v1;
36 g1 = gtemp;
37 u1 = utemp;
38 v1 = vtemp;
39 }
40 return g;
41}

Referenced by ZZp_INVERT().