OFFSET
1,1
COMMENTS
LINKS
A. Tripathi, On Pythagorean triples containing a fixed integer, Fib. Q., 46/47 (2008/2009), 331-340. See Theorem 8.
MATHEMATICA
t={}; Do[If[(PrimeQ[n]&&Mod[n, 4] == 3)||(PrimeQ[n/2]&&Mod[n/2, 4] == 3), t=Join[t, {n}]], {n, 470}]; t
(* Positions of the ones in A046081, omitting position = 4; based on program by Jean-François Alcover *)
a[1] = 0; a[n_] := Module[{f}, f = Select[FactorInteger[n], Mod[#[[1]], 4] == 1 &][[All, 2]]; (DivisorSigma[0, If[OddQ[n], n, n/2]^2] - 1)/2 + (Times @@ (2*f + 1) - 1)/2]; arr = Array[a, nmax]; fl = Flatten[Position[arr, 1]]; Delete[fl, 2]
CROSSREFS
KEYWORD
nonn
AUTHOR
Manfred Boergens, Jul 22 2024
STATUS
approved