OFFSET
0,2
COMMENTS
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Reinhard Zumkeller, Enumerations of Divisors
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1).
FORMULA
a(n) = C(n,0) + C(n,1) + 8*C(n,2) - 6*C(n,3).
G.f.: -(-1+2*x-9*x^2+14*x^3)/(-1+x)^4. - R. J. Mathar, Jun 18 2009
a(n) = 4*a(n-1)-6*a(n-2)+4*a(n-3)-a(n-4) with a(0)=1, a(1)=2, a(2)=11, a(3)=22. - Harvey P. Dale, Nov 12 2013
E.g.f.: (-x^3 + 4*x^2 + x + 1)*exp(x). - G. C. Greubel, Jul 16 2017
EXAMPLE
Differences of divisors of 22 to compute the coefficients of their interpolating polynomial, see formula:
1 2 11 22
1 9 11
8 2
-6
MATHEMATICA
Table[-n^3+7n^2-5n+1, {n, 0, 40}] (* or *) LinearRecurrence[{4, -6, 4, -1}, {1, 2, 11, 22}, 40] (* Harvey P. Dale, Nov 12 2013 *)
PROG
(Magma) [-n^3 + 7*n^2 - 5*n + 1: n in [0..40]]; // Vincenzo Librandi, Jul 17 2011
(PARI) a(n)=-n^3+7*n^2-5*n+1 \\ Charles R Greathouse IV, Sep 24 2015
CROSSREFS
KEYWORD
sign,easy
AUTHOR
Reinhard Zumkeller, Jun 17 2009
STATUS
approved