A036828 A036827/2.

0, 1, 17, 125, 637, 2637, 9549, 31501, 97037, 283661, 795661, 2158605, 5697549, 14696461, 37175309, 92471309, 226689037, 548667405, 1313079309, 3111125005, 7305429005, 17016291341, 39346765837, 90378862605, 206342979597, 468486979597, 1058239676429

Offset: 0

Author

N. J. A. Sloane

Keywords

• nonn
• easy

Comments

This sequence is related to A036826 by a(n) = n*A036826(n) - Sum_{i=0..n-1} A036826(i). - Bruno Berselli, Mar 06 2012

Links

• Bruno Berselli, Table of n, a(n) for n = 0..1000
• Index entries for linear recurrences with constant coefficients, signature (9,-32,56,-48,16).

Crossrefs

Cf. A036826, A209359.

Programs

Magma

m:=26; R:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!(-(4*x^2+8*x+1)/((x-1)*(2*x-1)^4))); // Bruno Berselli, Mar 06 2012

Mathematica

LinearRecurrence[{9, -32, 56, -48, 16}, {0, 1, 17, 125, 637}, 27] (* Bruno Berselli, Mar 06 2012 *)

PARI

a(n) = 2^n*(n^3-3*n^2+9*n-13)+13 \\ Bruno Berselli, Mar 06 2012

Formula

G.f.: -x*(4*x^2+8*x+1)/((x-1)*(2*x-1)^4). - Maksym Voznyy (voznyy(AT)mail.ru), Aug 13 2009

a(n) = 2^n*(n^3-3*n^2+9*n-13)+13. - Bruno Berselli, Mar 06 2012

Extensions

Typo in definition corrected by R. J. Mathar, Sep 16 2009