A197169 Values of gcd(n,y) for successive y = A197128(n).
2, 2, 2, 2, 2, 2, 2, 46, 2, 2, 2, 2, 3, 10, 6, 2, 3, 3, 2, 2, 2, 6, 5, 2, 2, 5, 2, 2, 2, 2, 2, 2, 22, 5, 2, 2, 3, 2, 2, 3, 2, 2, 3, 46, 5, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 14, 2, 4, 3, 2, 2, 2, 4, 14, 3, 2, 3, 2, 5, 2, 2, 2, 5, 2, 2, 2, 3, 6, 29, 3, 2, 2, 3
Offset: 1
Keywords
Programs
-
Mathematica
cr = {}; Do[If[IntegerQ[Sqrt[n]], , kk = NumberFieldFundamentalUnits[Sqrt[n]]; d1 = kk[[1]][[2]][[1]]; d2 = kk[[1]][[1]] kk[[1]][[2]][[2]]; d4 = Numerator[d2/Sqrt[n]]; If[GCD[d4, n] == 1, , AppendTo[cr, GCD[d4, n]]]], {n, 2, 20000}]; cr