Movatterモバイル変換


[0]ホーム

URL:


login
A299276
Partial sums ofA008137.
51
1, 5, 14, 31, 59, 101, 161, 242, 347, 479, 641, 837, 1070, 1343, 1659, 2021, 2433, 2898, 3419, 3999, 4641, 5349, 6126, 6975, 7899, 8901, 9985, 11154, 12411, 13759, 15201, 16741, 18382, 20127, 21979, 23941, 26017, 28210, 30523, 32959, 35521, 38213, 41038
OFFSET
0,2
COMMENTS
Euler transform of length 6 sequence [5, -1, 1, -1, 1, -1]. -Michael Somos, Oct 03 2018
FORMULA
FromColin Barker, Feb 11 2018: (Start)
G.f.: (1 + x)^3*(1 - x + x^2)*(1 + x^2) / ((1 - x)^4*(1 + x + x^2 + x^3 + x^4)).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + a(n-5) - 3*a(n-6) + 3*a(n-7) - a(n-8) for n>7.
(End)
a(n) = -a(-1-n) for all n in Z.
EXAMPLE
G.f. = 1 + 5*x + 14*x^2 + 31*x^3 + 59*x^4 + 101*x^5 + 161*x^6 + ... -Michael Somos, Oct 03 2018
MATHEMATICA
a[ n_] := (8 n^3 + 12 n^2 + 40 n + 18 - {3, 3, 0, -3, -3, 3}[[Mod[n, 5] + 1]]) / 15; (*Michael Somos, Oct 03 2018 *)
PROG
(PARI) Vec((1 + x)^3*(1 - x + x^2)*(1 + x^2) / ((1 - x)^4*(1 + x + x^2 + x^3 + x^4)) + O(x^60)) \\Colin Barker, Feb 11 2018
(PARI) {a(n) = (8*n^3 + 12*n^2 + 40*n + 18 - 3*(n%5<2) + 3*(n%5>2)) / 15}; /*Michael Somos, Oct 03 2018 */
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Feb 10 2018
STATUS
approved


[8]ページ先頭

©2009-2025 Movatter.jp