OFFSET
1,1
MAPLE
A000265 := proc(n) local nshft ; if n mod 2 = 1 then RETURN(n); else nshft := n ; while nshft mod 2 = 0 do nshft := nshft/2 ; od ; RETURN(nshft) ; fi ; end: A135013 := proc(n) option remember ; if n = 1 then RETURN(1) ; else RETURN( A135013(n-1)+A000265(n)) ; fi ; end: isA002378 := proc(n) RETURN( issqr(1+4*n) ) ; end: for n from 1 do a135013 := A135013(n) ; if isA002378(a135013) then printf("%d, ", a135013) ; fi ; od: # R. J. Mathar, Feb 11 2008
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Feb 10 2008
EXTENSIONS
More terms from R. J. Mathar, Feb 11 2008
a(15)-a(19) from Donovan Johnson, Aug 08 2010
a(20)-a(21) from Donovan Johnson, Jul 24 2011
STATUS
approved