OFFSET
0,2
COMMENTS
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
R. Zumkeller, Enumerations of Divisors
Index entries for linear recurrences with constant coefficients, signature (6,-15,20,-15,6,-1).
FORMULA
a(n) = C(n,0) + C(n,1) + 2*C(n,3) - 4*C(n,4) + 12*C(n,5).
G.f.: ( 1-4*x+6*x^2-2*x^3-7*x^4+18*x^5 ) / (x-1)^6. - R. J. Mathar, Jul 12 2016
EXAMPLE
Differences of divisors of 18 to compute the coefficients of their interpolating polynomial, see formula:
1 2 3 6 9 18
1 1 3 3 9
0 2 0 6
2 -2 6
-4 8
12
MAPLE
A161704:=n->(3*n^5 - 35*n^4 + 145*n^3 - 235*n^2 + 152*n + 30)/30: seq(A161704(n), n=0..50); # Wesley Ivan Hurt, Jul 16 2017
MATHEMATICA
CoefficientList[Series[(1 - 4*x + 6*x^2 - 2*x^3 - 7*x^4 + 18*x^5)/(x - 1)^6, {x, 0, 50}], x] (* G. C. Greubel, Jul 16 2017 *)
PROG
(Magma) [(3*n^5 - 35*n^4 + 145*n^3 - 235*n^2 + 152*n + 30)/30: n in [0..50]]; // Vincenzo Librandi, Dec 27 2010
(PARI) a(n)=n*(3*n^4-35*n^3+145*n^2-235*n+152)/30+1
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Reinhard Zumkeller, Jun 17 2009
STATUS
approved