OFFSET
1,1
COMMENTS
Subsequence of A191860. [Reinhard Zumkeller, Jun 18 2011]
MATHEMATICA
f[n_] := Block[{y = 1}, While[ !IntegerQ[ Sqrt[n*y^2 - 1]], y++ ]; Sqrt[n*y^2 - 1]]; lst = {}; Do[p = Prime@ n; If[ Mod[p, 4] == 1, AppendTo[lst, f@p]; Print[{n, Prime@n, f@p}]], {n, 66}]; lst
PROG
(Haskell)
a094048 n = head [m | m <- map (a037213 . subtract 1 . (* a002144 n))
(tail a000290_list), m > 0]
-- Reinhard Zumkeller, Jun 13 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Matthijs Coster, Apr 29 2004
EXTENSIONS
Edited by Don Reble, Apr 30 2004
STATUS
approved