0,2

For n>0, 30*a(n) is the sum of the ten distinct products of 2*n-1, 2*n+1, and 2*n+3. For example, when n = 1, we sum the ten distinct products of 1, 3, and 5: 1*1*1 + 1*1*3 + 1*1*5 + 1*3*3 + 1*3*5 + 1*5*5 + 3*3*3 + 3*3*5 + 3*5*5 + 5*5*5 = 330 = 30*11 = 30*a(1). -J. M. Bergot,Apr 06 2014

Vincenzo Librandi,Table of n, a(n) for n = 0..1000

T. P. Martin,Shells of atoms,Phys. Reports, 273 (1996), 199-241, eq. (10).

Index entries for linear recurrences with constant coefficients,signature (4,-6,4,-1).

a(n) = 2*A050533(n) + 1. -N. J. A. Sloane,Sep 22 2004

G.f.: (1+7*x+7*x^2+x^3)/(1-x)^4. -Colin Barker,Mar 01 2012

G.f. for sequence with interpolated zeros: 1/(8*x)*sinh(8*arctanh(x)) = 1/(16*x)*( ((1 + x)/(1 - x))^4 - ((1 - x)/(1 + x))^4 ) = 1 + 11*x^2 + 45*x^4 + 119*x^6 +.... Cf.A019560.-Peter Bala,Apr 07 2017

E.g.f.: (3 + 30*x + 36*x^2 + 8*x^3)*exp(x)/3. -G. C. Greubel,Dec 01 2017

FromPeter Bala,Mar 26 2024: (Start)

12*a(n) = (2*n + 1)*(a(n + 1) - a(n - 1)).

Sum_{n >= 0} (-1)^n/(a(n)*a(n+1)) = 3*Pi/16 - 1/2. Cf.A016754andA336266.(End)

A057813:=n->(2*n + 1)*(4*n^2 + 4*n + 3)/3; seq(A057813(n), n=0..50); #Wesley Ivan Hurt,Apr 06 2014

Table[(2*n + 1)*(4*n^2 + 4*n + 3)/3, {n, 0, 50}] (*David Nacin,Mar 01 2012 *)

(PARI) P(x, y, z) = x^3 + x^2*y + x^2*z + x*y^2 + x*y*z + x*z^2 + y^3 + y^2*z + y*z^2 + z^3;

a(n) = P(2*n-1, 2*n+1, 2*n+3)/30; \\Michel Marcus,Apr 22 2014

(Magma) [(2*n+1)*(4*n^2+4*n+3)/3: n in [0..50]] //Wesley Ivan Hurt,Apr 22 2014

1/12*t*(2*n^3-3*n^2+n)+2*n-1 for t = 2, 4, 6,... givesA049480,A005894,A063488,A001845,A063489,A005898,A063490,A057813,A063491,A005902,A063492,A005917,A063493,A063494,A063495,A063496.

Cf.A019560,A336266.

Sequence in context:A154106A232613A357736*A051740A370534A263227

Adjacent sequences:A057810A057811A057812*A057814A057815A057816

nonn,easy

N. J. A. Sloane,Nov 07 2000

approved