|
|
A006578
|
|
Triangular numbers plus quarter squares: n*(n+1)/2 + floor(n^2/4) (i.e., A000217(n) + A002620(n)).
(Formerly M3329)
|
|
36
|
|
|
0, 1, 4, 8, 14, 21, 30, 40, 52, 65, 80, 96, 114, 133, 154, 176, 200, 225, 252, 280, 310, 341, 374, 408, 444, 481, 520, 560, 602, 645, 690, 736, 784, 833, 884, 936, 990, 1045, 1102, 1160, 1220, 1281, 1344, 1408, 1474, 1541, 1610, 1680, 1752, 1825, 1900, 1976, 2054
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
Equals (1, 2, 3, 4, ...) convolved with (1, 2, 1, 2, ...). a(4) = 14 = (1, 2, 3, 4) dot (2, 1, 2, 1) = (2 + 2 + 6 + 4). - Gary W. Adamson, May 01 2009
We observe that is the transform of A032766 by the following transform T: T(u_0,u_1,u_2,u_3,...) = (u_0, u_0+u_1, u_0+u_1+u_2, u_0+u_1+u_2+u_3+u_4,...). In other words, v_p = Sum_{k=0..p} u_k and the g.f. phi_v of is given by phi_v = phi_u/(1-z). - Richard Choulet, Jan 28 2010
Equals row sums of a triangle with (1, 4, 7, 10, ...) in every column, shifted down twice for columns > 1. - Gary W. Adamson, Mar 03 2010
Number of pairs (x,y) with x in {0,...,n}, y odd in {0,...,2n}, and x < y. - Clark Kimberling, Jul 02 2012
|
|
REFERENCES
|
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
|
|
LINKS
|
|
|
FORMULA
|
Expansion of x*(1+2*x) / ((1-x)^2*(1-x^2)). - Simon Plouffe in his 1992 dissertation
a(n) = (6*n^2 + 4*n - 1 + (-1)^n)/8. - Paul Barry, May 30 2003
a(n) = Sum_{i=1..n} (n - i + 1) * 2^( (i+1) mod 2 ). - Wesley Ivan Hurt, Mar 30 2014
Sum_{n>=1} 1/a(n) = 3 - Pi/(4*sqrt(3)) - 3*log(3)/4. - Amiram Eldar, May 28 2022
E.g.f.: (x*(5 + 3*x)*cosh(x) - (1 - 5*x - 3*x^2)*sinh(x))/4. - Stefano Spezia, Aug 22 2023
|
|
EXAMPLE
|
G.f. = x + 4*x^2 + 8*x^3 + 14*x^4 + 21*x^5 + 30*x^6 + 40*x^7 + 52*x^8 + 65*x^9 + ...
|
|
MAPLE
|
with (combinat): seq(count(Partition((3*n+1)), size=3), n=0..52); # Zerinvary Lajos, Mar 28 2008
# 2nd program
(6*n^2 + 4*n - 1 + (-1)^n)/8 ;
|
|
MATHEMATICA
|
Accumulate[LinearRecurrence[{1, 1, -1}, {0, 1, 3}, 100]] (* Harvey P. Dale, Sep 29 2013 *)
a[ n_] := Quotient[n + 1, 2] (Quotient[n, 2] 3 + 1); (* Michael Somos, Jun 09 2014 *)
a[ n_] := Quotient[3 (n + 1)^2 + 1, 4] - (n + 1); (* Michael Somos, Jun 10 2015 *)
LinearRecurrence[{2, 0, -2, 1}, {0, 1, 4, 8}, 53] (* Ray Chandler, Aug 03 2015 *)
|
|
PROG
|
(PARI) {a(n) = (3*(n+1)^2 + 1)\4 - n - 1}; /* Michael Somos, Mar 10 2006 */
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|