OFFSET
1,2
MATHEMATICA
f[n_, z_] := FindInstance[x > 0 && y > 0 && 5*x^2 - n*y^2 == z, {x, y}, Integers, 1]; a[n_] := For[z = 0, z <= 5, z = z + GCD[n, 5], fz = f[n, z]; If[fz =!= {}, Print["a(", n, ") = ", z, " {x, y} = ", {x, y} /. fz[[1]]]; Return[z]]];
a[49] := (Print["a(49) = ", 5, " {x, y} = ", {51841, 16560}]; 5);
a[81] := (Print["a(81) = ", 5, " {x, y} = ", {161, 40}]; 5);
a[83] := (Print["a(83) = ", 5, " {x, y} = ", {1919, 471}]; 5);
a[98] := (Print["a(98) = ", 3, " {x, y} = ", {31, 7}]; 3);
Table[TimeConstrained[a[n], 300], {n, 1, 105}] (* Jean-François Alcover, Oct 10 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved