OFFSET
0,2
COMMENTS
Odd hexagonal numbers. Bisection of A000384. - Omar E. Pol, Apr 06 2008
Sequence found by reading the line from 1, in the direction 1, 15, ..., in the square spiral whose vertices are the triangular numbers A000217. - Omar E. Pol, Sep 03 2011
a(n) is also the sum of natural numbers which can be placed in a center box and expanded ones on 4 arms on N, S, E, W or NE, NW, SW, SE directions. See illustration in links. - Kival Ngaokrajang, Jul 08 2014
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..5000
Kival Ngaokrajang, Illustration of initial terms.
Leo Tavares, Illustration: Dual Square Stars
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).
FORMULA
a(n) = A157870(n)/2. - Vladimir Joseph Stephan Orlovsky, Mar 10 2009
a(n) = a(n-1) + 16*n-2 (with a(0)=1). - Vincenzo Librandi, Nov 20 2010
G.f.: (1+12*x+3*x^2)/(1-x)^3. - Colin Barker, Jan 08 2012
E.g.f.: (8*x^2 +14*x + 1)*exp(x). - G. C. Greubel, Jul 18 2017
From Amiram Eldar, Feb 28 2022: (Start)
Sum_{n>=0} 1/a(n) = Pi/4 + log(2)/2.
Sum_{n>=0} (-1)^n/a(n) = Pi*(sqrt(2)-1)/4 + log(sqrt(2)+1)/sqrt(2). (End)
MAPLE
MATHEMATICA
lst={}; Do[a=(2*n+1)*(4*n+1); AppendTo[lst, a], {n, 0, 5!}]; lst (* Vladimir Joseph Stephan Orlovsky, Mar 10 2009 *)
Table[(2 n + 1) (4 n + 1), {n, 0, 50}] (* Wesley Ivan Hurt, Jul 09 2014 *)
LinearRecurrence[{3, -3, 1}, {1, 15, 45}, 50] (* Harvey P. Dale, Aug 30 2021 *)
PROG
(Magma) [(2*n+1)*(4*n+1) : n in [0..50]]; // Wesley Ivan Hurt, Jul 09 2014
(PARI) a(n)=(2*n+1)*(4*n+1) \\ Charles R Greathouse IV, Sep 24 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
More terms from Wesley Ivan Hurt, Jul 09 2014
STATUS
approved