OFFSET
1,2
COMMENTS
A001221(a(n)) and a(n) are of opposite parity.
If m is in the sequence, then also 2*m.
Conjecture : a(n) is asymptotic to c*n where c is around 4 - Benoit Cloitre, Jan 06 2003
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
EXAMPLE
n = 99 = [(4*0+3)^2]*[(4*1+1)], therefore 99 is a term.
MATHEMATICA
fQ[n_]:=Plus@@((Mod[#[[1]], 4]-2)&/@If[n==1, {}, FactorInteger[n]])==0; Select[Range[240], fQ] (* Ray Chandler, Dec 18 2011*)
PROG
(Haskell)
a078613 n = a078613_list !! (n-1)
a078613_list = filter ((== 0) . a005094) [1..]
-- Reinhard Zumkeller, Jan 07 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Dec 10 2002
EXTENSIONS
Edited by Ray Chandler, Dec 18 2011
STATUS
approved