A060999 Nearest integer to (n+1)^3/9.

0, 1, 3, 7, 14, 24, 38, 57, 81, 111, 148, 192, 244, 305, 375, 455, 546, 648, 762, 889, 1029, 1183, 1352, 1536, 1736, 1953, 2187, 2439, 2710, 3000, 3310, 3641, 3993, 4367, 4764, 5184, 5628, 6097, 6591, 7111, 7658, 8232, 8834, 9465, 10125, 10815, 11536

Offset

0,3

Links

Harry J. Smith, Table of n, a(n) for n=0..1000

Index entries for linear recurrences with constant coefficients, signature (3,-3,2,-3,3,-1).

Formula

G.f.: x*(1+x^2)/((1-x)^3*(1-x^3)).

Expansion of x * (1 - x^4) / ((1 - x)^3 * (1 - x^2) * (1 - x^3)) in powers of x.

Euler transform of length 4 sequence [ 3, 1, 1, -1]. - Michael Somos, Aug 12 2009

a(-2-n) = -a(n). - Michael Somos, Aug 12 2009

G.f.: ( (1 + 4*x + x^2) / (1 - x)^4 - 1 / (1 + x + x^2) ) / 9.

Example

x + 3*x^2 + 7*x^3 + 14*x^4 + 24*x^5 + 38*x^6 + 57*x^7 + 81*x^8 + ...

Mathematica

Table[Floor[(n+1)^3/9+1/2], {n, 0, 50}] (* Harvey P. Dale, Jan 20 2013 *) (* or *)
LinearRecurrence[{3, -3, 2, -3, 3, -1}, {0, 1, 3, 7, 14, 24}, 47] (* Georg Fischer, Oct 13 2020 *)

Prog

(PARI) { default(realprecision, 100); for (n=0, 1000, write("b060999.txt", n, " ", round((n + 1)^3/9)) ) } \\ Harry J. Smith, Jul 16 2009
(PARI) {a(n) = n++; (n^3 - kronecker(-3, n)) / 9} /* Michael Somos, Aug 12 2009 */

Crossrefs

Sequence in context: A093523 A173247 A123386 * A089187 A333980 A316319

Adjacent sequences: A060996 A060997 A060998 * A061000 A061001 A061002

Keyword

nonn

Author

N. J. A. Sloane, May 14 2001

Status

approved