OFFSET
1,1
MATHEMATICA
a = {}; For[i = 1, i < 140, i++, k = 1; b = (2*i*(i + 1) + 1); While[b^(1/k) > 4, If[PrimeQ[b^(1/k)], If[Mod[b^(1/k), 4] == 1, AppendTo[a, b]]]; k++ ]]; a (* Stefan Steinerberger, Oct 21 2007 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Mats Granvik, Oct 18 2007, Oct 19 2007
EXTENSIONS
More terms from Stefan Steinerberger, Oct 21 2007
STATUS
approved