|
|
A008748
|
|
Expansion of (1 + x^5) / ((1-x) * (1-x^2) * (1-x^3)) in powers of x.
|
|
3
|
|
|
1, 1, 2, 3, 4, 6, 8, 10, 13, 16, 19, 23, 27, 31, 36, 41, 46, 52, 58, 64, 71, 78, 85, 93, 101, 109, 118, 127, 136, 146, 156, 166, 177, 188, 199, 211, 223, 235, 248, 261, 274, 288, 302, 316, 331, 346, 361, 377, 393, 409, 426, 443, 460, 478, 496, 514, 533, 552, 571
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
LINKS
|
|
|
FORMULA
|
a(n) = 1 + a(n-1) + a(n-3) - a(n-4) if n>4; a(n) = n if n=1..4. - Michael Somos, Jun 16 1999
|
|
EXAMPLE
|
G.f. = 1 + x + 2*x^2 + 3*x^3 + 4*x^4 + 6*x^5 + 8*x^6 + 10*x^7 + 13*x^8 + ...
|
|
MAPLE
|
A061347 := proc(n) op(1+(n mod 3), [-2, 1, 1]) ; end proc:
|
|
MATHEMATICA
|
CoefficientList[Series[(1+x^5)/((1-x)(1-x^2)(1-x^3)), {x, 0, 60}], x] (* Vincenzo Librandi, Jun 11 2013 *)
LinearRecurrence[{2, -1, 1, -2, 1}, {1, 1, 2, 3, 4}, 60] (* Harvey P. Dale, Apr 08 2019 *)
|
|
PROG
|
(Magma) [1 + Floor(n*(n+1)/6): n in [0..60]]; // G. C. Greubel, Aug 03 2019
(Sage) ((1 + x^5)/((1-x)*(1-x^2)*(1-x^3))).series(x, 30).coefficients(x, sparse=False) # G. C. Greubel, Aug 03 2019
(GAP) List([0..60], n-> 1 + Int(n*(n+1)/6)); # G. C. Greubel, Aug 03 2019
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|