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 19 results. Next

A033996 8 times triangular numbers: a(n) = 4*n*(n+1).

Original entry on oeis.org

0, 8, 24, 48, 80, 120, 168, 224, 288, 360, 440, 528, 624, 728, 840, 960, 1088, 1224, 1368, 1520, 1680, 1848, 2024, 2208, 2400, 2600, 2808, 3024, 3248, 3480, 3720, 3968, 4224, 4488, 4760, 5040, 5328, 5624, 5928, 6240, 6560, 6888, 7224, 7568, 7920, 8280
Offset: 0

Views

Author

N. J. A. Sloane, Dec 11 1999

Keywords

Comments

Write 0, 1, 2, ... in a clockwise spiral; sequence gives numbers on one of 4 diagonals.
Also, least m > n such that T(m)*T(n) is a square and more precisely that of A055112(n). {T(n) = A000217(n)}. - Lekraj Beedassy, May 14 2004
Also sequence found by reading the line from 0, in the direction 0, 8, ... and the same line from 0, in the direction 0, 24, ..., in the square spiral whose vertices are the generalized decagonal numbers A074377. Axis perpendicular to A195146 in the same spiral. - Omar E. Pol, Sep 18 2011
Number of diagonals with length sqrt(5) in an (n+1) X (n+1) square grid. Every 1 X 2 rectangle has two such diagonals. - Wesley Ivan Hurt, Mar 25 2015
Imagine a board made of squares (like a chessboard), one of whose squares is completely surrounded by square-shaped layers made of adjacent squares. a(n) is the total number of squares in the first to n-th layer. a(1) = 8 because there are 8 neighbors to the unit square; adding them gives a 3 X 3 square. a(2) = 24 = 8 + 16 because we need 16 more squares in the next layer to get a 5 X 5 square: a(n) = (2*n+1)^2 - 1 counting the (2n+1) X (2n+1) square minus the central square. - R. J. Cano, Sep 26 2015
The three platonic solids (the simplex, hypercube, and cross-polytope) with unit side length in n dimensions all have rational volume if and only if n appears in this sequence, after 0. - Brian T Kuhns, Feb 26 2016
The number of active (ON, black) cells in the n-th stage of growth of the two-dimensional cellular automaton defined by "Rule 645", based on the 5-celled von Neumann neighborhood. - Robert Price, May 19 2016
The square root of a(n), n>0, has continued fraction [2n; {1,4n}] with whole number part 2n and periodic part {1,4n}. - Ron Knott, May 11 2017
Numbers k such that k+1 is a square and k is a multiple of 4. - Bruno Berselli, Sep 28 2017
a(n) is the number of vertices of the octagonal network O(n,n); O(m,n) is defined by Fig. 1 of the Siddiqui et al. reference. - Emeric Deutsch, May 13 2018
a(n) is the number of vertices in conjoined n X n octagons which are arranged into a square array, a.k.a. truncated square tiling. - Donghwi Park, Dec 20 2020
a(n-2) is the number of ways to place 3 adjacent marks in a diagonal, horizontal, or vertical row on an n X n tic-tac-toe grid. - Matej Veselovac, May 28 2021

Examples

			Spiral with 0, 8, 24, 48, ... along lower right diagonal:
.
  36--37--38--39--40--41--42
   |                       |
  35  16--17--18--19--20  43
   |   |               |   |
  34  15   4---5---6  21  44
   |   |   |       |   |   |
  33  14   3   0   7  22  45
   |   |   |   | \ |   |   |
  32  13   2---1   8  23  46
   |   |           | \ |   |
  31  12--11--10---9  24  47
   |                   | \ |
  30--29--28--27--26--25  48
                            \
[Reformatted by _Jon E. Schoenfield_, Dec 25 2016]
		

References

  • Stuart M. Ellerstein, J. Recreational Math. 29 (3) 188, 1998.
  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 2nd ed., 1994, p. 99.
  • Stephen Wolfram, A New Kind of Science, Wolfram Media, 2002; p. 170.

Crossrefs

Cf. A000217, A016754, A002378, A024966, A027468, A028895, A028896, A045943, A046092, A049598, A088538, A124080, A008590 (first differences), A130809 (partial sums).
Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.

Programs

  • Magma
    [ 4*n*(n+1) : n in [0..50] ]; // Wesley Ivan Hurt, Jun 09 2014
  • Maple
    seq(8*binomial(n+1, 2), n=0..46); # Zerinvary Lajos, Nov 24 2006
    [seq((2*n+1)^2-1, n=0..46)];
  • Mathematica
    Table[(2n - 1)^2 - 1, {n, 50}] (* Alonso del Arte, Mar 31 2013 *)
  • PARI
    nsqm1(n) = { forstep(x=1,n,2, y = x*x-1; print1(y, ", ") ) }
    

Formula

a(n) = 4*n^2 + 4*n = (2*n+1)^2 - 1.
G.f.: 8*x/(1-x)^3.
a(n) = A016754(n) - 1 = 2*A046092(n) = 4*A002378(n). - Lekraj Beedassy, May 25 2004
a(n) = A049598(n) - A046092(n); a(n) = A124080(n) - A002378(n). - Zerinvary Lajos, Mar 06 2007
a(n) = 8*A000217(n). - Omar E. Pol, Dec 12 2008
a(n) = A005843(n) * A163300(n). - Juri-Stepan Gerasimov, Jul 26 2009
a(n) = a(n-1) + 8*n (with a(0)=0). - Vincenzo Librandi, Nov 17 2010
For n > 0, a(n) = A058031(n+1) - A062938(n-1). - Charlie Marion, Apr 11 2013
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Wesley Ivan Hurt, Mar 25 2015
a(n) = A000578(n+1) - A152618(n). - Bui Quang Tuan, Apr 01 2015
a(n) - a(n-1) = A008590(n), n > 0. - Altug Alkan, Sep 26 2015
From Ilya Gutkovskiy, May 19 2016: (Start)
E.g.f.: 4*x*(2 + x)*exp(x).
Sum_{n>=1} 1/a(n) = 1/4. (End)
Product_{n>=1} a(n)/A016754(n) = Pi/4. - Daniel Suteu, Dec 25 2016
a(n) = A056220(n) + A056220(n+1). - Bruce J. Nicholson, May 29 2017
sqrt(a(n)+1) - sqrt(a(n)) = (sqrt(n+1) - sqrt(n))^2. - Seiichi Manyama, Dec 23 2018
a(n)*a(n+k) + 4*k^2 = m^2 where m = (a(n) + a(n+k))/2 - 2*k^2; for k=1, m = 4*n^2 + 8*n + 2 = A060626(n). - Ezhilarasu Velayutham, May 22 2019
Sum_{n>=1} (-1)^n/a(n) = 1/4 - log(2)/2. - Vaclav Kotesovec, Dec 21 2020
From Amiram Eldar, Feb 21 2023: (Start)
Product_{n>=1} (1 - 1/a(n)) = -(4/Pi)*cos(Pi/sqrt(2)).
Product_{n>=1} (1 + 1/a(n)) = 4/Pi (A088538). (End)

A122141 Array: T(d,n) = number of ways of writing n as a sum of d squares, read by ascending antidiagonals.

Original entry on oeis.org

1, 1, 2, 1, 4, 0, 1, 6, 4, 0, 1, 8, 12, 0, 2, 1, 10, 24, 8, 4, 0, 1, 12, 40, 32, 6, 8, 0, 1, 14, 60, 80, 24, 24, 0, 0, 1, 16, 84, 160, 90, 48, 24, 0, 0, 1, 18, 112, 280, 252, 112, 96, 0, 4, 2, 1, 20, 144, 448, 574, 312, 240, 64, 12, 4, 0, 1, 22, 180, 672, 1136, 840, 544, 320, 24, 30, 8, 0
Offset: 1

Views

Author

R. J. Mathar, Oct 29 2006

Keywords

Comments

This is the transpose of the array in A286815.
T(d,n) is divisible by 2d for any n != 0 iff d is a power of 2. - Jianing Song, Sep 05 2018

Examples

			Array T(d,n) with rows d = 1,2,3,... and columns n = 0,1,2,3,... reads
  1  2   0   0    2    0     0     0     0     2      0 ...
  1  4   4   0    4    8     0     0     4     4      8 ...
  1  6  12   8    6   24    24     0    12    30     24 ...
  1  8  24  32   24   48    96    64    24   104    144 ...
  1 10  40  80   90  112   240   320   200   250    560 ...
  1 12  60 160  252  312   544   960  1020   876   1560 ...
  1 14  84 280  574  840  1288  2368  3444  3542   4424 ...
  1 16 112 448 1136 2016  3136  5504  9328 12112  14112 ...
  1 18 144 672 2034 4320  7392 12672 22608 34802  44640 ...
  1 20 180 960 3380 8424 16320 28800 52020 88660 129064 ...
		

Crossrefs

Cf. A000122 (1st row), A004018 (2nd row), A005875 (3rd row), A000118 (4th row), A000132 (5th row), A000141 (6th row), A008451 (7th row), A000143 (8th row), A008452 (9th row), A000144 (10th row), A008453 (11th row), A000145 (12th row), A276285 (13th row), A276286 (14th row), A276287 (15th row), A000152 (16th row).
Cf. A005843 (2nd column), A046092 (3rd column), A130809 (4th column).
Cf. A010052 (1st row divides 2), A002654 (2nd row divides 4), A046897 (4th row divides 8), A008457 (8th row divides 16), A302855 (16th row divides 32), A302857 (32nd row divides 64).

Programs

  • Maple
    A122141 := proc(d,n) local i,cnts ; cnts := 0 ; for i from -trunc(sqrt(n)) to trunc(sqrt(n)) do if n-i^2 >= 0 then if d > 1 then cnts := cnts+procname(d-1,n-i^2) ; elif n-i^2 = 0 then cnts := cnts+1 ; fi ; fi ; od ; cnts ;
    end:
    for diag from 1 to 14 do for n from 0 to diag-1 do d := diag-n ; printf("%d,",A122141(d,n)) ; od ; od;
    # second Maple program:
    A:= proc(d, n) option remember; `if`(n=0, 1, `if`(n<0 or d<1, 0,
          A(d-1, n) +2*add(A(d-1, n-j^2), j=1..isqrt(n))))
        end:
    seq(seq(A(h-n, n), n=0..h-1), h=1..14); # Alois P. Heinz, Jul 16 2014
  • Mathematica
    Table[ SquaresR[d - n, n], {d, 1, 12}, {n, 0, d - 1}] // Flatten (* Jean-François Alcover, Jun 13 2013 *)
    A[d_, n_] := A[d, n] = If[n==0, 1, If[n<0 || d<1, 0, A[d-1, n] + 2*Sum[A[d-1, n-j^2], {j, 1, Sqrt[n]}]]]; Table[A[h-n, n], {h, 1, 14}, {n, 0, h-1}] // Flatten (* Jean-François Alcover, Feb 28 2018, after Alois P. Heinz *)
  • Python
    from sympy.core.power import isqrt
    from functools import cache
    @cache
    def T(d, n):
      if n == 0: return 1
      if n < 0 or d < 1: return 0
      return T(d-1, n) + sum(T(d-1, n-(j**2)) for j in range(1, isqrt(n)+1)) * 2  # Darío Clavijo, Feb 06 2024

Formula

T(n,n) = A066535(n). - Alois P. Heinz, Jul 16 2014

A319574 A(n, k) = [x^k] JacobiTheta3(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, 2, 1, 0, 0, 4, 1, 0, 0, 4, 6, 1, 0, 2, 0, 12, 8, 1, 0, 0, 4, 8, 24, 10, 1, 0, 0, 8, 6, 32, 40, 12, 1, 0, 0, 0, 24, 24, 80, 60, 14, 1, 0, 0, 0, 24, 48, 90, 160, 84, 16, 1, 0, 2, 4, 0, 96, 112, 252, 280, 112, 18, 1, 0, 0, 4, 12, 64, 240, 312, 574, 448, 144, 20, 1
Offset: 0

Views

Author

Peter Luschny, Oct 01 2018

Keywords

Comments

Number of ways of writing k as a sum of n squares.

Examples

			[ 0] 1,  0,    0,    0,     0,     0,     0      0,     0,     0, ... A000007
[ 1] 1,  2,    0,    0,     2,     0,     0,     0,     0,     2, ... A000122
[ 2] 1,  4,    4,    0,     4,     8,     0,     0,     4,     4, ... A004018
[ 3] 1,  6,   12,    8,     6,    24,    24,     0,    12,    30, ... A005875
[ 4] 1,  8,   24,   32,    24,    48,    96,    64,    24,   104, ... A000118
[ 5] 1, 10,   40,   80,    90,   112,   240,   320,   200,   250, ... A000132
[ 6] 1, 12,   60,  160,   252,   312,   544,   960,  1020,   876, ... A000141
[ 7] 1, 14,   84,  280,   574,   840,  1288,  2368,  3444,  3542, ... A008451
[ 8] 1, 16,  112,  448,  1136,  2016,  3136,  5504,  9328, 12112, ... A000143
[ 9] 1, 18,  144,  672,  2034,  4320,  7392, 12672, 22608, 34802, ... A008452
[10] 1, 20,  180,  960,  3380,  8424, 16320, 28800, 52020, 88660, ... A000144
   A005843,   v, A130809,  v,  A319576,  v ,   ...      diagonal: A066535
           A046092,    A319575,       A319577,     ...
		

References

  • E. Grosswald, Representations of Integers as Sums of Squares. Springer-Verlag, NY, 1985, p. 121.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954.
  • J. Carlos Moreno and Samuel S. Wagstaff Jr., Sums Of Squares Of Integers, Chapman & Hall/CRC, (2006).

Crossrefs

Variant starting with row 1 is A122141, transpose of A286815.

Programs

  • Maple
    A319574row := proc(n, len) series(JacobiTheta3(0, x)^n, x, len+1);
    [seq(coeff(%, x, j), j=0..len-1)] end:
    seq(print([n], A319574row(n, 10)), n=0..10);
    # Alternative, uses function PMatrix from A357368.
    PMatrix(10, n -> A000122(n-1)); # Peter Luschny, Oct 19 2022
  • Mathematica
    A[n_, k_] := If[n == k == 0, 1, SquaresR[n, k]];
    Table[A[n-k, k], {n, 0, 11}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Nov 03 2018 *)
  • Sage
    for n in (0..10):
        Q = DiagonalQuadraticForm(ZZ, [1]*n)
        print(Q.theta_series(10).list())

A130810 If X_1,...,X_n is a partition of a 2n-set X into 2-blocks then a(n) is equal to the number of 4-subsets of X containing none of X_i, (i=1,...,n).

Original entry on oeis.org

16, 80, 240, 560, 1120, 2016, 3360, 5280, 7920, 11440, 16016, 21840, 29120, 38080, 48960, 62016, 77520, 95760, 117040, 141680, 170016, 202400, 239200, 280800, 327600, 380016, 438480, 503440, 575360, 654720, 742016, 837760, 942480, 1056720
Offset: 4

Views

Author

Milan Janjic, Jul 16 2007

Keywords

Comments

Number of n permutations (n>=4) of 3 objects u,v,z, with repetition allowed, containing n-4 u's. Example: if n=4 then n-4 =(0) zero u, a(1)=16 because we have vvvv zzzz vvvz zzzv vvzv zzvz vzvv zvzz zvvv vzzz vvzz zzvv vzvz zvzv zvvz vzzv. - Zerinvary Lajos, Aug 05 2008
a(n) is the number of 3-dimensional elements in an n-cross polytope where n>=4. - Patrick J. McNab, Jul 06 2015

Crossrefs

Programs

  • Maple
    a:= n-> binomial(2*n,4) +binomial(n,2) -n*binomial(2*n-2,2);
    seq(binomial(n, n-4)*2^4, n=4..37); # Zerinvary Lajos, Dec 07 2007
  • Mathematica
    a[n_] := 16 * Binomial[n, 4]; Array[a, 34, 4] (* Amiram Eldar, Jul 25 2025 *)
  • PARI
    a(n) = 16 * binomial(n, 4); \\ Amiram Eldar, Jul 25 2025

Formula

a(n) = binomial(2*n,4) + binomial(n,2) - n*binomial(2*n-2,2).
a(n) = binomial(n,4)*16. - Zerinvary Lajos, Dec 07 2007
G.f.: 16*x^4/(1-x)^5. - Colin Barker, Apr 14 2012
a(n) = 2*n*(n-1)*(n-2)*(n-3)/3 = 2*A162668(n-3). - Robert Israel, Jul 06 2015
a(n) = 16 * A000332(n). - Alois P. Heinz, Oct 26 2020
E.g.f.: 2*exp(x)*x^4/3. - Stefano Spezia, Jul 17 2025
From Amiram Eldar, Jul 25 2025: (Start)
Sum_{n>=4} 1/a(n) = 1/12.
Sum_{n>=4} (-1)^n/a(n) = 2*log(2) - 4/3. (End)

A174794 a(0) = 0 and a(n) = (4*n^3 - 12*n^2 + 20*n - 9)/3 for n >= 1.

Original entry on oeis.org

0, 1, 5, 17, 45, 97, 181, 305, 477, 705, 997, 1361, 1805, 2337, 2965, 3697, 4541, 5505, 6597, 7825, 9197, 10721, 12405, 14257, 16285, 18497, 20901, 23505, 26317, 29345, 32597, 36081, 39805, 43777, 48005, 52497, 57261, 62305, 67637, 73265, 79197, 85441, 92005, 98897
Offset: 0

Views

Author

Roger L. Bagula, Mar 29 2010

Keywords

Comments

For n >= 1, a(n+1) = (4*n^3 + 8*n + 3)/3 is the number of evaluation points on the n-dimensional cube in Stenger's degree 7 cubature rule. - Franck Maminirina Ramaharo, Dec 18 2018

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[x*(1 + x)*(1 + 3*x^2)/(1 - x)^4, {x, 0, 50}], x]
  • Maxima
    a[0] : 0$ a[n] := (4*n^3 - 12*n^2 + 20*n - 9)/3$ makelist(a[n], n, 0, 50); /* Martin Ettl, Oct 21 2012 */

Formula

G.f.: x*(1 + x)*(1 + 3*x^2)/(1 - x)^4.
From Franck Maminirina Ramaharo, Dec 17 2018: (Start)
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4), n >= 5.
a(n) = 8*binomial(n - 1, 3) + 8*binomial(n - 1, 2) + 4*binomial(n - 1, 1) + 1, n >= 1.
E.g.f.: (9 - (9 - 12*x - 4*x^3)*exp(x))/3. (End)

Extensions

Definition replaced by polynomial - The Assoc. Eds. of the OEIS, Aug 10 2010

A349454 Number T(n,k) of endofunctions on [n] with exactly k fixed points, all of which are isolated; triangle T(n,k), n>=0, 0<=k<=n, read by rows.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 8, 3, 0, 1, 81, 32, 6, 0, 1, 1024, 405, 80, 10, 0, 1, 15625, 6144, 1215, 160, 15, 0, 1, 279936, 109375, 21504, 2835, 280, 21, 0, 1, 5764801, 2239488, 437500, 57344, 5670, 448, 28, 0, 1, 134217728, 51883209, 10077696, 1312500, 129024, 10206, 672, 36, 0, 1
Offset: 0

Views

Author

Alois P. Heinz, Dec 30 2021

Keywords

Examples

			Triangle T(n,k) begins:
        1;
        0,       1;
        1,       0,      1;
        8,       3,      0,     1;
       81,      32,      6,     0,    1;
     1024,     405,     80,    10,    0,   1;
    15625,    6144,   1215,   160,   15,   0,  1;
   279936,  109375,  21504,  2835,  280,  21,  0, 1;
  5764801, 2239488, 437500, 57344, 5670, 448, 28, 0, 1;
  ...
		

Crossrefs

Column k=0 gives A065440.
Row sums give A204042.
Main diagonal and first lower diagonal give A000012, A000004.
T(n+1,n-1) gives A000217.
T(n+3,n) gives A130809.
T(n+3,n-1) gives A102741 for n>=1.

Programs

  • Maple
    T:= (n, k)-> binomial(n, k)*(n-k-1)^(n-k):
    seq(seq(T(n, k), k=0..n), n=0..10);

Formula

T(n,k) = binomial(n,k) * (n-k-1)^(n-k).
From Mélika Tebni, Apr 02 2023: (Start)
E.g.f. of column k: -x / (LambertW(-x)*(1+LambertW(-x)))*x^k / k!.
Sum_{k=0..n} k^k*T(n,k) = A217701(n). (End)

A130812 If X_1,...,X_n is a partition of a 2n-set X into 2-blocks then a(n) is equal to the number of 6-subsets of X containing none of X_i, (i=1,...n).

Original entry on oeis.org

64, 448, 1792, 5376, 13440, 29568, 59136, 109824, 192192, 320320, 512512, 792064, 1188096, 1736448, 2480640, 3472896, 4775232, 6460608, 8614144, 11334400, 14734720, 18944640, 24111360, 30401280, 38001600, 47121984, 57996288, 70884352, 86073856, 103882240
Offset: 6

Views

Author

Milan Janjic, Jul 16 2007

Keywords

Comments

Number of n permutations (n>=6) of 3 objects u,v,z, with repetition allowed, containing n-6 u's. Example: if n=6 then n-6 =(0) zero u, a(1)=64. - Zerinvary Lajos, Aug 05 2008
a(n) is the number of 5-dimensional elements in an n-cross polytope where n>=6. - Patrick J. McNab, Jul 06 2015

Crossrefs

Programs

  • Magma
    [Binomial(2*n,6)+Binomial(n,2)*Binomial(2*n-4,2)- n*Binomial(2*n-2,4)-Binomial(n,3): n in [6..40]]; // Vincenzo Librandi, Jul 09 2015
  • Maple
    a:=n->binomial(2*n,6)+binomial(n,2)*binomial(2*n-4,2)-n*binomial(2*n-2,4)-binomial(n,3);
    seq(binomial(n,n-6)*2^6,n=6..32); # Zerinvary Lajos, Dec 07 2007
    seq(binomial(n+5, 6)*2^6, n=1..22); # Zerinvary Lajos, Aug 05 2008
  • Mathematica
    CoefficientList[Series[64/(1-x)^7,{x,0,30}],x] (* Vincenzo Librandi, Mar 21 2012 *)

Formula

a(n) = binomial(2*n,6) + binomial(n,2)*binomial(2*n-4,2) - n*binomial(2*n-2,4) - binomial(n,3).
a(n) = C(n,n-6)*2^6, n>=6. - Zerinvary Lajos, Dec 07 2007
G.f.: 64*x^6/(1-x)^7. - Colin Barker, Mar 20 2012

A210440 a(n) = 2*n*(n+1)*(n+2)/3.

Original entry on oeis.org

0, 4, 16, 40, 80, 140, 224, 336, 480, 660, 880, 1144, 1456, 1820, 2240, 2720, 3264, 3876, 4560, 5320, 6160, 7084, 8096, 9200, 10400, 11700, 13104, 14616, 16240, 17980, 19840, 21824, 23936, 26180, 28560, 31080, 33744, 36556, 39520, 42640, 45920, 49364, 52976
Offset: 0

Views

Author

Michel Marcus, Jan 20 2013

Keywords

Comments

Number of tin boxes necessary to build a tetrahedron with side length n, as shown in the link.
If "0" is prepended, a(n) is the convolution of 2n with itself. - Wesley Ivan Hurt, Mar 14 2015

Crossrefs

Cf. A000292, A028552, A033488 (partial sums), A046092, A130809.

Programs

Formula

a(n) = 4*A000292(n).
a(n+1)-a(n) = A046092(n+1).
From Bruno Berselli, Jan 20 2013: (Start)
G.f.: 4*x/(1-x)^4.
a(n) = -a(-n-2) = 4*a(n-1)-6*a(n-2)+4*a(n-3)-a(n-4).
a(n)-a(-n) = A217873(n).
a(n)+a(-n) = A016742(n).
(n-1)*a(n)-n*a(n-1) = A130809(n+1) with n>1. (End)
From Bruno Berselli, Jan 21 2013: (Start)
a(n) = n*A028552(n) - Sum_{i=0..n-1} A028552(i) for n>0.
4*A001296(n) = n*a(n) - Sum_{i=0..n-1} a(i) for n>0. (End)
G.f.: 2*x*W(0), where W(k) = 1 + 1/(1 - x*(k+2)*(k+4)/(x*(k+2)*(k+4) + (k+1)*(k+2)/W(k+1))); (continued fraction). - Sergei N. Gladkovskii, Aug 24 2013
a(n) = Sum_{i=1..n} i*(2*n-i+3). - Wesley Ivan Hurt, Oct 03 2013
From Amiram Eldar, Apr 30 2023: (Start)
Sum_{n>=1} 1/a(n) = 3/8.
Sum_{n>=1} (-1)^(n+1)/a(n) = 3*log(2) - 15/8. (End)
E.g.f.: 2*exp(x)*x*(6 + 6*x + x^2)/3. - Stefano Spezia, Jul 11 2025

A319575 a(n) = (2/3)*n*(n^3 - 6*n^2 + 11*n - 3).

Original entry on oeis.org

0, 2, 4, 6, 24, 90, 252, 574, 1136, 2034, 3380, 5302, 7944, 11466, 16044, 21870, 29152, 38114, 48996, 62054, 77560, 95802, 117084, 141726, 170064, 202450, 239252, 280854, 327656, 380074, 438540, 503502, 575424, 654786, 742084, 837830, 942552, 1056794, 1181116
Offset: 0

Views

Author

Peter Luschny, Oct 01 2018

Keywords

Crossrefs

Cf. A000012 (m=0), A005843 (m=1), A046092 (m=2), A130809 (m=3), this sequence (m=4), A319576 (m=5), A319577 (m=6).
Column n=4 of A122141.
Cf. A319574.

Programs

  • Maple
    a := n -> (2/3)*n*(n^3 - 6*n^2 + 11*n - 3):
    seq(a(n), n=0..38);
  • Mathematica
    A319575[n_] := 2/3*n*(n^3-6*n^2+11*n-3); Array[A319575, 50, 0] (* or *)
    LinearRecurrence[{5, -10, 10, -5, 1}, {0, 2, 4, 6, 24}, 50] (* Paolo Xausa, Feb 20 2024 *)
  • PARI
    concat(0, Vec(2*x*(1 + x)*(1 - 4*x + 7*x^2) / (1 - x)^5 + O(x^40))) \\ Colin Barker, Oct 02 2018

Formula

a(n) = [x^4] JacobiTheta3(x)^n.
a(n) = A319574(n,4).
From Colin Barker, Oct 02 2018: (Start)
G.f.: 2*x*(1 + x)*(1 - 4*x + 7*x^2) / (1 - x)^5.
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5) for n>4. (End)

A319576 a(n) = (4/15)*n*(n - 1)*(n^3 - 9*n^2 + 26*n - 9).

Original entry on oeis.org

0, 0, 8, 24, 48, 112, 312, 840, 2016, 4320, 8424, 15224, 25872, 41808, 64792, 96936, 140736, 199104, 275400, 373464, 497648, 652848, 844536, 1078792, 1362336, 1702560, 2107560, 2586168, 3147984, 3803408, 4563672, 5440872, 6448000, 7598976, 8908680, 10392984
Offset: 0

Views

Author

Peter Luschny, Oct 01 2018

Keywords

Crossrefs

Cf. A000012 (m=0), A005843 (m=1), A046092 (m=2), A130809 (m=3), A319575 (m=4), this sequence (m=5), A319577 (m=6).
Column n=5 of A122141.
Cf. A319574.

Programs

  • Maple
    a := n -> (4/15)*n*(n - 1)*(n^3 - 9*n^2 + 26*n - 9):
    seq(a(n), n=0..41);
  • Mathematica
    A319576[n_] := 4/15*n*(n-1)*(n^3-9*n^2+26*n-9); Array[A319576, 50, 0] (* or *)
    LinearRecurrence[{6, -15, 20, -15, 6, -1}, {0, 0, 8, 24, 48, 112}, 50] (* Paolo Xausa, Feb 20 2024 *)
  • PARI
    concat([0,0], Vec(8*x^2*(1 - 3*x + 3*x^2 + 3*x^3) / (1 - x)^6 + O(x^40))) \\ Colin Barker, Oct 02 2018

Formula

a(n) = [x^5] JacobiTheta3(x)^n.
a(n) = A319574(n,5).
From Colin Barker, Oct 02 2018: (Start)
G.f.: 8*x^2*(1 - 3*x + 3*x^2 + 3*x^3) / (1 - x)^6.
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.
(End)
Showing 1-10 of 19 results. Next