OFFSET
0,1
REFERENCES
I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, Wiley, N.Y., 1983, (2.3.14).
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (6,-15,20,-15,6,-1).
FORMULA
G.f.: (x+2)*(2*x^4 + 8*x^3 + 36*x^2 + 31*x + 4)/(1-x)^6. Generally, g.f. for k-th column of A046741 is coefficient of y^k in expansion of (1-y)/((1 - y - y^2)*(1-y) - (1+y)*x).
From G. C. Greubel, Jan 31 2019: (Start)
a(n) = (320 + 1114*n + 1515*n^2 + 1125*n^3 + 405*n^4 + 81*n^5)/40.
E.g.f.: (320 + 4240*x + 8940*x^2 + 5580*x^3 + 1215*x^4 + 81*x^5)*exp(x)/40. (End)
MATHEMATICA
Table[(320+1114*n+1515*n^2+1125*n^3+405*n^4+81*n^5)/40, {n, 0, 40}] (* G. C. Greubel, Jan 31 2019 *)
PROG
(PARI) vector(40, n, n--; (320+1114*n+1515*n^2+1125*n^3+405*n^4 + 81*n^5)/40) \\ G. C. Greubel, Jan 31 2019
(Magma) [(320+1114*n+1515*n^2+1125*n^3+405*n^4+81*n^5)/40: n in [0..40]]; // G. C. Greubel, Jan 31 2019
(Sage) [(320+1114*n+1515*n^2+1125*n^3+405*n^4+81*n^5)/40 for n in range(40)] # G. C. Greubel, Jan 31 2019
(GAP) List([0..40], n -> (320+1114*n+1515*n^2+1125*n^3+405*n^4+81*n^5 )/40); # G. C. Greubel, Jan 31 2019
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Vladeta Jovovic, Jun 04 2001
EXTENSIONS
More terms from Larry Reeves (larryr(AT)acm.org), Jun 06 2001
STATUS
approved