OFFSET
1,1
MATHEMATICA
f[n_, z_] := FindInstance[x > 0 && y > 0 && 3*x^2 - n*y^2 == z, {x, y}, Integers, 1]; a[n_] := For[z = 0, True, z = z + GCD[n, 3], fz = f[n, z]; If[fz =!= {}, Print["a(", n, ") = ", z, " {x, y} = ", {x, y} /. fz[[1]]]; Return[z]]]; a[32] := (Print["a(32) = ", 3, " {x, y} = ", {4801, 1470}]; 3); a[64] := (Print["a(64) = ", 3, " {x, y} = ", {97, 21}]; 3); a[101] := (Print["a(101) = ", 3, " {x, y} = ", {2524, 435}]; 3); a[103] := (Print["a(103) = ", 3, " {x, y} = ", {64202725495, 10957096332}]; 3); Table[TimeConstrained[a[n], 300], {n, 1, 105}] (* Jean-François Alcover, Oct 09 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved