cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 15 results. Next

A317614 a(n) = (1/2)*(n^3 + n*(n mod 2)).

Original entry on oeis.org

1, 4, 15, 32, 65, 108, 175, 256, 369, 500, 671, 864, 1105, 1372, 1695, 2048, 2465, 2916, 3439, 4000, 4641, 5324, 6095, 6912, 7825, 8788, 9855, 10976, 12209, 13500, 14911, 16384, 17985, 19652, 21455, 23328, 25345, 27436, 29679, 32000, 34481, 37044, 39775, 42592
Offset: 1

Views

Author

Stefano Spezia, Aug 01 2018

Keywords

Comments

Terms are obtained as partial sums in an algorithm for the generation of the sequence of the fourth powers (A000583). Starting with the sequence of the positive integers (A000027), it is necessary to delete every 4th term and to consider the partial sums of the obtained sequence, then to delete every 3rd term, and lastly to consider again the partial sums (see References).
a(n) is the trace of an n X n square matrix M(n) formed by writing the numbers 1, ..., n^2 successively forward and backward along the rows in zig-zag pattern as shown in the examples below. Specifically, M(n) is defined as M[i,j,n] = j + n*(i-1) if i is odd and M[i,j,n] = n*i - j + 1 if i is even, and it has det(M(n)) = 0 for n > 2 (proved).
From Saeed Barari, Oct 31 2021: (Start)
Also the sum of the entries in an n X n matrix whose elements start from 1 and increase as they approach the center. For instance, in case of n=5, the entries of the following matrix sum to 65:
1 2 3 2 1
2 3 4 3 2
3 4 5 4 3
2 3 4 3 2
1 2 3 2 1. (End)
The n X n square matrix of the preceding comment is defined as: A[i,j,n] = n - abs((n + 1)/2 - j) - abs((n + 1)/2 - i). - Stefano Spezia, Nov 05 2021

Examples

			For n = 1 the matrix M(1) is
  1
with trace Tr(M(1)) = a(1) = 1.
For n = 2 the matrix M(2) is
  1, 2
  4, 3
with Tr(M(2)) = a(2) = 4.
For n = 3 the matrix M(3) is
  1, 2, 3
  6, 5, 4
  7, 8, 9
with Tr(M(3)) = a(3) = 15.
		

References

  • Edward A. Ashcroft, Anthony A. Faustini, Rangaswami Jagannathan, and William W. Wadge, Multidimensional Programming, Oxford University Press 1995, p. 12.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 64.
  • G. Polya, Mathematics and Plausible Reasoning: Induction and analogy in mathematics, Princeton University Press 1990, p. 118.
  • Shailesh Shirali, A Primer on Number Sequences, Universities Press (India) 2004, p. 106.
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, Exercise 3.7.3 on pages 122-123.

Crossrefs

Cf. A000583, A000027, A186424 (first differences).
Cf. related to the M matrices: A074147 (antidiagonals), A130130 (rank), A241016 (row sums), A317617 (column sums), A322277 (permanent), A323723 (subdiagonal sums), A323724 (superdiagonal sums).

Programs

  • GAP
    a_n:=List([1..nmax], n->(1/2)*(n^3 + n*RemInt(n, 2)));
    
  • GAP
    List([1..50],n->(1/2)*(n^3+n*(n mod 2))); # Muniru A Asiru, Aug 24 2018
  • Magma
    [IsEven(n) select n^3/2 else (n^3+n)/2: n in [1..50]]; // Vincenzo Librandi, Aug 07 2018
    
  • Maple
    a:=n->(1/2)*(n^3+n*modp(n,2)): seq(a(n),n=1..50); # Muniru A Asiru, Aug 24 2018
  • Mathematica
    CoefficientList[Series[1/4 E^-x (1 + 3 E^(2 x) + 6 E^(2 x) x + 2 E^(2 x) x^2), {x, 0, 45}], x]*Table[(k + 1)!, {k, 0, 45}]
    CoefficientList[Series[-(1 + x^2)/((-1 + x)*(1 + x)^3), {x, 0, 45}], x]*Table[(k + 1)*(-1)^k, {k, 0, 45}]
    CoefficientList[Series[-(1 + x^2)/((-1 + x)^3*(1 + x)), {x, 0, 45}], x]*Table[(k + 1), {k, 0, 45}]
    From Robert G. Wilson v, Aug 01 2018: (Start)
    a[i_, j_, n_] := If[OddQ@ i, j + n (i - 1), n*i - j + 1]; f[n_] := Tr[Table[a[i, j, n], {i, n}, {j, n}]]; Array[f, 45]
    CoefficientList[Series[(x^4 + 2x^3 + 6x^2 + 2x + 1)/((x - 1)^4 (x + 1)^2), {x, 0, 45}], x]
    LinearRecurrence[{2, 1, -4, 1, 2, -1}, {1, 4, 15, 32, 65, 108}, 45]
    (End)
  • Maxima
    a(n):=(1/2)*(n^3 + n*mod(n,2))$ makelist(a(n), n, 1, nmax);
    
  • PARI
    Vec(x*(1 + 2*x + 6*x^2 + 2*x^3 + x^4) / ((1 - x)^4*(1 + x)^2) + O(x^40)) \\ Colin Barker, Aug 02 2018
    
  • PARI
    M(i, j, n) = if (i % 2, j + n*(i-1), n*i - j + 1);
    a(n) = sum(k=1, n, M(k, k, n)); \\ Michel Marcus, Aug 07 2018
    
  • R
    for (n in 1:nmax){
       a <- (n^3+n*n%%2)/2
       output <- c(n, a)
       cat(output, "\n")
    }
    (MATLAB and FreeMat)
    for(n=1:nmax); a=(n^3+n*mod(n,2))/2; fprintf('%d\t%0.f\n',n,a); end
    

Formula

a(n) = (1/2)*(A000578(n) + n*A000035(n)).
a(n) = A006003(n) - (n/2)*(1 - (n mod 2)).
a(n) = Sum_{k=1..n} T(n,k), where T(n,k) = ((n + 1)*k - n)*(n mod 2) + ((n - 1)*k + 1)*(1 - (n mod 2)).
E.g.f.: E(x) = (1/4)*exp(-x)*x*(1 + 3*exp(2*x) + 6*exp(2*x)*x + 2*exp(2*x)*x^2).
L.g.f.: L(x) = -x*(1 + x^2)/((-1 + x)*(1 + x)^3).
H.l.g.f.: LH(x) = -x*(1 + x^2)/((-1 + x)^3*(1 + x)).
Dirichlet g.f.: (1/2)*(Zeta(-3 + s) + 2^(-s)*(-2 + 2^s)*Zeta(-1 + s)).
From Colin Barker, Aug 02 2018: (Start)
G.f.: x*(1 + 2*x + 6*x^2 + 2*x^3 + x^4) / ((1 - x)^4*(1 + x)^2).
a(n) = 2*a(n-1) + a(n-2) - 4*a(n-3) + a(n-4) + 2*a(n-5) - a(n-6) for n>6.
a(n) = n^3/2 for n even.
a(n) = (n^3+n)/2 for n odd. (End)
a(2*n) = A317297(n+1) + A001489(n). - Stefano Spezia, Dec 28 2018
Sum_{n>0} 1/a(n) = (1/2)*(-2*polygamma(0, 1/2) + polygamma(0, (1-i)/2)+ polygamma(0, (1+i)/2)) + zeta(3)/4 approximately equal to 1.3959168891658447368440622669882813003351669... - Stefano Spezia, Feb 11 2019
a(n) = (A000578(n) + A193356(n))/2. - Stefano Spezia, Jun 27 2022
a(n) = A210378(n-1)/n. - Stefano Spezia, Jul 15 2024

A024000 a(n) = 1 - n.

Original entry on oeis.org

1, 0, -1, -2, -3, -4, -5, -6, -7, -8, -9, -10, -11, -12, -13, -14, -15, -16, -17, -18, -19, -20, -21, -22, -23, -24, -25, -26, -27, -28, -29, -30, -31, -32, -33, -34, -35, -36, -37, -38, -39, -40, -41, -42, -43, -44, -45, -46, -47, -48, -49, -50, -51, -52, -53, -54, -55, -56
Offset: 0

Views

Author

Keywords

Comments

a(n) is the weighted sum over all derangements (permutations with no fixed points) of n elements where each permutation with an odd number of cycles has weight +1 and each with an even number of cycles has weight -1. [Michael Somos, Jan 19 2011]

Examples

			a(4) = -3 because there are 6 derangements with one 4-cycle with weight -1 and 3 derangements with two 2-cycles with weight +1. - _Michael Somos_, Jan 19 2011
		

Crossrefs

A022958 shifted left.

Programs

  • Magma
    [1-n: n in [0..50]]; // Vincenzo Librandi, Apr 29 2011
  • Maple
    A024000:=n->1-n: seq(A024000(n), n=0..100); # Wesley Ivan Hurt, Mar 02 2016
  • Mathematica
    CoefficientList[Series[(1 - 2 x)/(1 - x)^2, {x, 0, 60}], x] Range[0, 60]!
    CoefficientList[Series[Exp[x] (1 - x), {x, 0, 60}], x]
    1-Range[0,60] (* Harvey P. Dale, Sep 18 2013 *)
    Flatten[NestList[(#/.x_/;x>1->Sequence[x,2x])-1&,{1},60]]
    (* Robert G. Wilson v, Mar 02 2016 *)
  • PARI
    {a(n) = 1 - n} /* Michael Somos, Jan 19 2011 */
    

Formula

E.g.f.: (1-x)*exp(x).
a(n) = Sum_{k=0..n} A094816(n,k)*(-1)^k (alternating row sums of Poisson-Charlier coefficient matrix).
O.g.f.: (1-2*x)/(1-x)^2. a(n+1) = A001489(n). - R. J. Mathar, May 28 2008
a(n) = 2*a(n-1)-a(n-2) for n>1. - Wesley Ivan Hurt, Mar 02 2016

A269044 a(n) = 13*n + 7.

Original entry on oeis.org

7, 20, 33, 46, 59, 72, 85, 98, 111, 124, 137, 150, 163, 176, 189, 202, 215, 228, 241, 254, 267, 280, 293, 306, 319, 332, 345, 358, 371, 384, 397, 410, 423, 436, 449, 462, 475, 488, 501, 514, 527, 540, 553, 566, 579, 592, 605, 618, 631, 644, 657, 670, 683, 696, 709, 722, 735
Offset: 0

Views

Author

Bruno Berselli, Feb 18 2016

Keywords

Comments

After 7 (which corresponds to n=0), all terms belong to A090767 because a(n) = 3*n*2*1 + 2*(n*2+2*1+n*1) + (n+2+1).
This sequence is related to A152741 by the recurrence A152741(n+1) = (n+1)*a(n+1) - Sum_{k = 0..n} a(k).
Any square mod 13 is one of 0, 1, 3, 4, 9, 10 or 12 (A010376) but not 7, and for this reason there are no squares in the sequence. Likewise, any cube mod 13 is one of 0, 1, 5, 8 or 12, therefore no a(k) is a cube.
The sum of the squares of any two terms of the sequence is also a term of the sequence, that is: a(h)^2 + a(k)^2 = a(h*(13*h+14) + k*(13*k+14) + 7). Therefore: a(h)^2 + a(k)^2 > a(a( h*(h+1) + k*(k+1) )) for h+k > 0.
The primes of the sequence are listed in A140371.

Crossrefs

Cf. A010376, A022271 (partial sums), A088227, A090767, A140371, A152741.
Similar sequences with closed form (2*k-1)*n+k: A001489 (k=0), A000027 (k=1), A016789 (k=2), A016885 (k=3), A017029 (k=4), A017221 (k=5), A017461 (k=6), this sequence (k=7), A164284 (k=8).
Sequences of the form 13*n+q: A008595 (q=0), A190991 (q=1), A153080 (q=2), A127547 (q=4), A154609 (q=5), A186113 (q=6), this sequence (q=7), A269100 (q=11).

Programs

  • Magma
    [13*n+7: n in [0..60]];
    
  • Mathematica
    13 Range[0, 60] + 7 (* or *) Range[7, 800, 13] (* or *) Table[13 n + 7, {n, 0, 60}]
    LinearRecurrence[{2, -1}, {7, 20}, 60] (* Vincenzo Librandi, Feb 19 2016 *)
  • Maxima
    makelist(13*n+7, n, 0, 60);
    
  • PARI
    vector(60, n, n--; 13*n+7)
    
  • Sage
    [13*n+7 for n in (0..60)]

Formula

G.f.: (7 + 6*x)/(1 - x)^2.
a(n) = A088227(4*n+3).
a(n) = -A186113(-n-1).
Sum_{i=h..h+13*k} a(i) = a(h*(13*k + 1) + k*(169*k + 27)/2).
Sum_{i>=0} 1/a(i)^2 = 0.0257568950542502716970... = polygamma(1, 7/13)/13^2.
E.g.f.: exp(x)*(7 + 13*x). - Stefano Spezia, Aug 02 2021

A001478 The negative integers.

Original entry on oeis.org

-1, -2, -3, -4, -5, -6, -7, -8, -9, -10, -11, -12, -13, -14, -15, -16, -17, -18, -19, -20, -21, -22, -23, -24, -25, -26, -27, -28, -29, -30, -31, -32, -33, -34, -35, -36, -37, -38, -39, -40, -41, -42, -43, -44, -45, -46, -47, -48, -49, -50, -51, -52, -53, -54, -55, -56, -57, -58, -59, -60, -61, -62, -63, -64, -65
Offset: 1

Views

Author

Keywords

Crossrefs

Cf. A001477 (complement in A001057), A000027, A001489, A130472.

Programs

Formula

a(n) = -n.
G.f.: -x/(1-x)^2.
G.f. A(x) satisfies A(x) + A(-x) = 4A(x^2), A(x)A(-x) = A(x^2), A(x)^2 + A(x^2) = 4A(x)A(x^2). - Michael Somos, Mar 23 2004
a(n) = -A000027(n). - Felix Fröhlich, Jul 23 2021

A286509 Square array A(n,k), n>=0, k>=0, read by antidiagonals, where column k is the expansion of k-th power of continued fraction 1/(1 + x/(1 + x^2/(1 + x^3/(1 + x^4/(1 + x^5/(1 + ...)))))).

Original entry on oeis.org

1, 1, 0, 1, -1, 0, 1, -2, 1, 0, 1, -3, 3, 0, 0, 1, -4, 6, -2, -1, 0, 1, -5, 10, -7, -1, 1, 0, 1, -6, 15, -16, 3, 4, -1, 0, 1, -7, 21, -30, 15, 6, -6, 1, 0, 1, -8, 28, -50, 40, 0, -17, 6, 0, 0, 1, -9, 36, -77, 84, -26, -30, 24, -3, -1, 0, 1, -10, 45, -112, 154, -90, -30, 64, -21, -2, 2, 0, 1, -11, 55, -156, 258, -217, 15, 125, -81, 6, 9, -3, 0
Offset: 0

Views

Author

Ilya Gutkovskiy, May 10 2017

Keywords

Examples

			Square array begins:
1,  1,  1,  1,   1,   1,  ...
0, -1, -2, -3,  -4,  -5,  ...
0,  1,  3,  6,  10,  15,  ...
0,  0, -2, -7, -16, -30,  ...
0, -1, -1,  3,  15,  40,  ...
0,  1,  4,  6,   0, -26,  ...
		

Crossrefs

Columns k=0-5 give: A000007, A007325, A055101, A055102, A055103, A078905 (with offset 0).
Rows n=0-2 give: A000012, A001489, A000217.
Main diagonal gives A291651.
Antidiagonal sums give A302015.

Programs

  • Mathematica
    Table[Function[k, SeriesCoefficient[1/(1 + ContinuedFractionK[x^i, 1, {i, 1, n}])^k, {x, 0, n}]][j - n], {j, 0, 12}, {n, 0, j}] // Flatten
    Table[Function[k, SeriesCoefficient[Product[(1 - x^(5 i - 1)) (1 - x^(5 i - 4))/((1 - x^(5 i - 2)) (1 - x^(5 i - 3))), {i, n}]^k, {x, 0, n}]][j - n], {j, 0, 12},{n, 0, j}] // Flatten

Formula

G.f. of column k: Product_{j>=1} ((1 - x^(5*j-1))*(1 - x^(5*j-4)) / ((1 - x^(5*j-2))*(1 - x^(5*j-3))))^k.

A057428 Sign(-n): a(n) = 1 if -n > 0, = -1 if -n < 0, = 0 if n = 0.

Original entry on oeis.org

0, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1
Offset: 0

Views

Author

Henry Bottomley, Sep 05 2000

Keywords

Crossrefs

First differences of A001489.
Cf. A057427.

Programs

Formula

G.f.: A(x) = -x/(1-x).
a(n) = -A057427(n).
Series reversion of g.f. A(x) is A(A(x)). - Nikolaos Pantelidis, Apr 13 2022

A319933 A(n, k) = [x^k] DedekindEta(x)^n, square array read by descending antidiagonals, A(n, k) for n >= 0 and k >= 0.

Original entry on oeis.org

1, 0, 1, 0, -1, 1, 0, -1, -2, 1, 0, 0, -1, -3, 1, 0, 0, 2, 0, -4, 1, 0, 1, 1, 5, 2, -5, 1, 0, 0, 2, 0, 8, 5, -6, 1, 0, 1, -2, 0, -5, 10, 9, -7, 1, 0, 0, 0, -7, -4, -15, 10, 14, -8, 1, 0, 0, -2, 0, -10, -6, -30, 7, 20, -9, 1, 0, 0, -2, 0, 8, -5, 0, -49, 0, 27, -10, 1
Offset: 0

Views

Author

Peter Luschny, Oct 02 2018

Keywords

Comments

The columns are generated by polynomials whose coefficients constitute the triangle of signed D'Arcais numbers A078521 when multiplied with n!.

Examples

			[ 0] 1,   0,   0,    0,     0,    0,     0,     0,     0,     0, ... A000007
[ 1] 1,  -1,  -1,    0,     0,    1,     0,     1,     0,     0, ... A010815
[ 2] 1,  -2,  -1,    2,     1,    2,    -2,     0,    -2,    -2, ... A002107
[ 3] 1,  -3,   0,    5,     0,    0,    -7,     0,     0,     0, ... A010816
[ 4] 1,  -4,   2,    8,    -5,   -4,   -10,     8,     9,     0, ... A000727
[ 5] 1,  -5,   5,   10,   -15,   -6,    -5,    25,    15,   -20, ... A000728
[ 6] 1,  -6,   9,   10,   -30,    0,    11,    42,     0,   -70, ... A000729
[ 7] 1,  -7,  14,    7,   -49,   21,    35,    41,   -49,  -133, ... A000730
[ 8] 1,  -8,  20,    0,   -70,   64,    56,     0,  -125,  -160, ... A000731
[ 9] 1,  -9,  27,  -12,   -90,  135,    54,   -99,  -189,   -85, ... A010817
[10] 1, -10,  35,  -30,  -105,  238,     0,  -260,  -165,   140, ... A010818
    A001489,  v , A167541, v , A319931,  v ,         diagonal: A008705
           A080956       A319930      A319932
		

References

  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. Fifth ed., Clarendon Press, Oxford, 2003.

Crossrefs

Transpose of A286354.
Cf. A078521, A319574 (JacobiTheta3).

Programs

  • Julia
    # DedekindEta is defined in A000594
    for n in 0:10
        DedekindEta(10, n) |> println
    end
  • Maple
    DedekindEta := (x, n) -> mul(1-x^j, j=1..n):
    A319933row := proc(n, len) series(DedekindEta(x, len)^n, x, len+1):
    seq(coeff(%, x, j), j=0..len-1) end:
    seq(print([n], A319933row(n, 10)), n=0..10);
  • Mathematica
    eta[x_, n_] := Product[1 - x^j, {j, 1, n}];
    A[n_, k_] := SeriesCoefficient[eta[x, k]^n, {x, 0, k}];
    Table[A[n - k, k], {n, 0, 11}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Nov 10 2018 *)
  • Sage
    from sage.modular.etaproducts import qexp_eta
    def A319933row(n, len):
        return (qexp_eta(ZZ['q'], len+4)^n).list()[:len]
    for n in (0..10):
        print(A319933row(n, 10))
    

A305548 a(n) = 27*n.

Original entry on oeis.org

0, 27, 54, 81, 108, 135, 162, 189, 216, 243, 270, 297, 324, 351, 378, 405, 432, 459, 486, 513, 540, 567, 594, 621, 648, 675, 702, 729, 756, 783, 810, 837, 864, 891, 918, 945, 972, 999, 1026, 1053, 1080, 1107, 1134, 1161, 1188, 1215, 1242, 1269, 1296, 1323, 1350, 1377, 1404, 1431, 1458, 1485, 1512
Offset: 0

Views

Author

Eric Chen, Jun 05 2018

Keywords

Crossrefs

For a(n) = k*n: A001489 (k=-1), A000004 (k=0), A001477 (k=1), A005843 (k=2), A008585 (k=3), A008591 (k=9), A008607 (k=25), A252994 (k=26), this sequence (k=27), A135628 (k=28), A195819 (k=29), A249674 (k=30), A135631 (k=31), A174312 (k=32), A044102 (k=36), A085959 (k=37), A169823 (k=60), A152691 (k=64).

Programs

  • Mathematica
    Range[0,2000,27]
  • PARI
    a(n)=27*n

Formula

a(n) = 27*n.
a(n) = A008585(A008591(n)) = A008591(A008585(n)).
G.f.: 27*x/(x-1)^2.
From Elmo R. Oliveira, Apr 10 2025: (Start)
E.g.f.: 27*x*exp(x).
a(n) = 2*a(n-1) - a(n-2). (End)

A319930 a(n) = (1/24)*n*(n - 1)*(n - 3)*(n - 14).

Original entry on oeis.org

0, 0, 1, 0, -5, -15, -30, -49, -70, -90, -105, -110, -99, -65, 0, 105, 260, 476, 765, 1140, 1615, 2205, 2926, 3795, 4830, 6050, 7475, 9126, 11025, 13195, 15660, 18445, 21576, 25080, 28985, 33320, 38115, 43401, 49210, 55575, 62530, 70110, 78351, 87290, 96965
Offset: 0

Views

Author

Peter Luschny, Oct 02 2018

Keywords

Crossrefs

Cf. A000012 (m=0), A001489 (m=1), A080956 (m=2), A167541 (m=3), this sequence (m=4), A319931 (m=5), A319932 (m=6).
Cf. A319933.

Programs

  • Maple
    a := n -> (1/24)*n*(n-1)*(n-3)*(n-14):
    seq(a(n), n=0..44);
  • Mathematica
    Table[(n(n-1)(n-3)(n-14))/24,{n,0,70}] (* Harvey P. Dale, Apr 29 2022 *)

Formula

a(n) = [x^4] DedekindEta(x)^n.
a(n) = A319933(n, 4).

A319931 a(n) = -(1/120)*n*(n - 3)*(n - 6)*(n^2 - 21*n + 8).

Original entry on oeis.org

0, 1, 2, 0, -4, -6, 0, 21, 64, 135, 238, 374, 540, 728, 924, 1107, 1248, 1309, 1242, 988, 476, -378, -1672, -3519, -6048, -9405, -13754, -19278, -26180, -34684, -45036, -57505, -72384, -89991, -110670, -134792, -162756, -194990, -231952, -274131, -322048, -376257
Offset: 0

Views

Author

Peter Luschny, Oct 02 2018

Keywords

Crossrefs

Cf. A000012 (m=0), A001489 (m=1), A080956 (m=2), A167541 (m=3), A319930 (m=4), this sequence (m=5), A319932 (m=6).
Cf. A319933.

Programs

  • Maple
    a := n -> -(1/120)*n*(n-3)*(n-6)*(n^2-21*n+8):
    seq(a(n), n=0..41);
  • PARI
    a(n)=-n*(n-3)*(n-6)*(n^2-21*n+8)/120 \\ Charles R Greathouse IV, Oct 21 2022

Formula

a(n) = [x^5] DedekindEta(x)^n.
a(n) = A319933(n, 5).
From Chai Wah Wu, Jul 27 2022: (Start)
a(n) = 6*a(n-1) - 15*a(n-2) + 20*a(n-3) - 15*a(n-4) + 6*a(n-5) - a(n-6) for n > 5.
G.f.: x*(-7*x^4 + 6*x^3 + 3*x^2 - 4*x + 1)/(x - 1)^6. (End)
Showing 1-10 of 15 results. Next