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

A002417 4-dimensional figurate numbers: a(n) = n*binomial(n+2, 3).

Original entry on oeis.org

1, 8, 30, 80, 175, 336, 588, 960, 1485, 2200, 3146, 4368, 5915, 7840, 10200, 13056, 16473, 20520, 25270, 30800, 37191, 44528, 52900, 62400, 73125, 85176, 98658, 113680, 130355, 148800, 169136, 191488, 215985, 242760, 271950, 303696, 338143, 375440, 415740
Offset: 1

Views

Author

Keywords

Comments

a(n) is 1/6 the number of colorings of a 2 X 2 hexagonal array with n+2 colors. - R. H. Hardin, Feb 23 2002
a(n) is the sum of all numbers that cannot be written as t*(n+1) + u*(n+2) for nonnegative integers t,u (see Schuh). - Floor van Lamoen, Oct 09 2002
a(n) is the total number of rectangles (including squares) contained in a stepped pyramid of n rows (or of base 2n-1) of squares. A stepped pyramid of squares of base 2*6 - 1 = 11, for instance, has the following vertices:
..........X.X
........X.X.X.X
......X.X.X.X.X.X
....X.X.X.X.X.X.X.X
..X.X.X.X.X.X.X.X.X.X
X.X.X.X.X.X.X.X.X.X.X.X
X.X.X.X.X.X.X.X.X.X.X.X - Lekraj Beedassy, Sep 02 2003
Partial sums of A002412. - Jonathan Vos Post, Mar 16 2006
a(n) equals -1 times the coefficient of x^3 of the characteristic polynomial of the (n + 2) X (n + 2) matrix with 2's along the main diagonal and 1's everywhere else (see Mathematica code below). - John M. Campbell, May 28 2011
a(n) is the n-th antidiagonal sum of the convolution array A213750. - Clark Kimberling, Jun 20 2012
Convolution of A000027 with A000384 (excluding 0). - Bruno Berselli, Dec 06 2012
The sequence is the binomial transform of (1, 7, 15, 13, 4, 0, 0, 0, ...). - Gary W. Adamson, Jul 31 2015
Also the number of 3-cycles in the (n+2)-triangular graph. - Eric W. Weisstein, Aug 14 2017

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 195.
  • K. -W. Lau, Solution to Problem 2495, Journal of Recreational Mathematics 2002-3 31(1) 79-80.
  • Fred. Schuh, Vragen betreffende een onbepaalde vergelijking, Nieuw Tijdschrift voor Wiskunde, 52 (1964-1965) 193-198.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Bisection of A002624.
a(n) = A093561(n+3, 4).
Cf. A220212 for a list of sequences produced by the convolution of the natural numbers with the k-gonal numbers.
Cf. A151974 (number of 4-cycles in the triangular graph), A290939 (5-cycles), A290940 (6-cycles).

Programs

  • GAP
    List([1..40], n-> n^2*(n+1)*(n+2)/6 ); # G. C. Greubel, Jul 03 2019
  • Magma
    /* A000027 convolved with A000384 (excluding 0): */ A000384:=func; [&+[(n-i+1)*A000384(i): i in [1..n]]: n in [1..40]]; // Bruno Berselli, Dec 06 2012
    
  • Magma
    [n*Binomial(n+2,3):n in [1..40]]; // Vincenzo Librandi, Aug 02 2015
    
  • Maple
    seq(n^2*(n+1)*(n+2)/6, n=1..50);
  • Mathematica
    Table[n Binomial[n+2, 3], {n, 40}]
    Table[-Coefficient[CharacteristicPolynomial[Array[KroneckerDelta[#1, #2] + 1 &, {n+2, n+2}], x], x^3], {n, 40}] (* John M. Campbell, May 28 2011 *)
    Nest[Accumulate, Range[1, 170, 4], 3] (* Vladimir Joseph Stephan Orlovsky, Jan 21 2012 *)
    LinearRecurrence[{5, -10, 10, -5, 1}, {1, 8, 30, 80, 175}, 40] (* Harvey P. Dale, Jan 11 2014 *)
    Table[n Pochhammer[n, 3]/6, {n, 40}] (* or *) CoefficientList[Series[ (1+3x)/(1-x)^5, {x,0,40}], x] (* Eric W. Weisstein, Aug 14 2017 *)
  • PARI
    a(n)=n^2*(n+1)*(n+2)/6 \\ Charles R Greathouse IV, Jun 10 2011
    
  • Sage
    [n*binomial(n+2,3) for n in (1..40)] # G. C. Greubel, Jul 03 2019
    

Formula

a(n) = n^2*(n+1)*(n+2)/6.
G.f.: x*(1+3*x)/(1-x)^5. - Simon Plouffe in his 1992 dissertation
a(n) = C(n+2, 2)*n^2/3. - Paul Barry, Jun 26 2003
a(n) = C(n+3, n)*C(n+1, 1). - Zerinvary Lajos, Apr 27 2005
a(n) = (binomial(n+3,n-1) - binomial(n+2,n-2))*(binomial(n+1,n-1) - binomial(n,n-2)). - Zerinvary Lajos, May 12 2006
a(n) = 5*a(n-1) -10*a(n-2) +10*a(n-3) -5*a(n-4) +a(n-5), n>5. - Wesley Ivan Hurt, Aug 01 2015
G.f.: x*hypergeometric2F1(2,4;1;x). - R. J. Mathar, Aug 09 2015
a(n) = A080852(4,n-1). - R. J. Mathar, Jul 28 2016
Sum_{n>=1} 1/a(n) = Pi^2/2 - 15/4. - Jaume Oliver Lafont, Jul 13 2017
E.g.f.: x*(6 + 18*x + 9*x^2 + x^3)*exp(x)/3!. - G. C. Greubel, Jul 03 2019
Sum_{n>=1} (-1)^(n+1)/a(n) = (Pi^2 + 27 - 48*log(2))/4. - Amiram Eldar, Jun 28 2020
a(n) = A000332(n+3) + 3*A000332(n+2). - Mircea Dan Rus, Jul 29 2020

Extensions

Edited and extended by Floor van Lamoen, Oct 09 2002

A139600 Square array T(n,k) = n*(k-1)*k/2+k, of nonnegative numbers together with polygonal numbers, read by antidiagonals upwards.

Original entry on oeis.org

0, 0, 1, 0, 1, 2, 0, 1, 3, 3, 0, 1, 4, 6, 4, 0, 1, 5, 9, 10, 5, 0, 1, 6, 12, 16, 15, 6, 0, 1, 7, 15, 22, 25, 21, 7, 0, 1, 8, 18, 28, 35, 36, 28, 8, 0, 1, 9, 21, 34, 45, 51, 49, 36, 9, 0, 1, 10, 24, 40, 55, 66, 70, 64, 45, 10, 0, 1, 11, 27, 46, 65, 81, 91, 92, 81, 55, 11
Offset: 0

Views

Author

Omar E. Pol, Apr 27 2008

Keywords

Comments

A general formula for polygonal numbers is P(n,k) = (n-2)*(k-1)*k/2 + k, where P(n,k) is the k-th n-gonal number.
The triangle sums, see A180662 for their definitions, link this square array read by antidiagonals with twelve different sequences, see the crossrefs. Most triangle sums are linear sums of shifted combinations of a sequence, see e.g. A189374. - Johannes W. Meijer, Apr 29 2011

Examples

			The square array of nonnegatives together with polygonal numbers begins:
=========================================================
....................... A   A   .   .   A    A    A    A
....................... 0   0   .   .   0    0    1    1
....................... 0   0   .   .   1    1    3    3
....................... 0   0   .   .   6    7    9    9
....................... 0   0   .   .   9    3    6    6
....................... 0   1   .   .   5    2    0    0
....................... 4   2   .   .   7    9    6    7
=========================================================
Nonnegatives . A001477: 0,  1,  2,  3,  4,   5,   6,   7, ...
Triangulars .. A000217: 0,  1,  3,  6, 10,  15,  21,  28, ...
Squares ...... A000290: 0,  1,  4,  9, 16,  25,  36,  49, ...
Pentagonals .. A000326: 0,  1,  5, 12, 22,  35,  51,  70, ...
Hexagonals ... A000384: 0,  1,  6, 15, 28,  45,  66,  91, ...
Heptagonals .. A000566: 0,  1,  7, 18, 34,  55,  81, 112, ...
Octagonals ... A000567: 0,  1,  8, 21, 40,  65,  96, 133, ...
9-gonals ..... A001106: 0,  1,  9, 24, 46,  75, 111, 154, ...
10-gonals .... A001107: 0,  1, 10, 27, 52,  85, 126, 175, ...
11-gonals .... A051682: 0,  1, 11, 30, 58,  95, 141, 196, ...
12-gonals .... A051624: 0,  1, 12, 33, 64, 105, 156, 217, ...
...
=========================================================
The column with the numbers 2, 3, 4, 5, 6, ... is formed by the numbers > 1 of A000027. The column with the numbers 3, 6, 9, 12, 15, ... is formed by the positive members of A008585.
		

Crossrefs

A formal extension negative n is in A326728.
Triangle sums (see the comments): A055795 (Row1), A080956 (Row2; terms doubled), A096338 (Kn11, Kn12, Kn13, Fi1, Ze1), A002624 (Kn21, Kn22, Kn23, Fi2, Ze2), A000332 (Kn3, Ca3, Gi3), A134393 (Kn4), A189374 (Ca1, Ze3), A011779 (Ca2, Ze4), A101357 (Ca4), A189375 (Gi1), A189376 (Gi2), A006484 (Gi4). - Johannes W. Meijer, Apr 29 2011
Sequences of m-gonal numbers: A000217 (m=3), A000290 (m=4), A000326 (m=5), A000384 (m=6), A000566 (m=7), A000567 (m=8), A001106 (m=9), A001107 (m=10), A051682 (m=11), A051624 (m=12), A051865 (m=13), A051866 (m=14), A051867 (m=15), A051868 (m=16), A051869 (m=17), A051870 (m=18), A051871 (m=19), A051872 (m=20), A051873 (m=21), A051874 (m=22), A051875 (m=23), A051876 (m=24), A255184 (m=25), A255185 (m=26), A255186 (m=27), A161935 (m=28), A255187 (m=29), A254474 (m=30).

Programs

  • Magma
    T:= func< n,k | k*(n*(k-1)+2)/2 >;
    A139600:= func< n,k | T(n-k, k) >;
    [A139600(n,k): k in  [0..n], n in [0..12]]; // G. C. Greubel, Jul 12 2024
    
  • Maple
    T:= (n, k)-> n*(k-1)*k/2+k:
    seq(seq(T(d-k, k), k=0..d), d=0..14);  # Alois P. Heinz, Oct 14 2018
  • Mathematica
    T[n_, k_] := (n + 1)*(k - 1)*k/2 + k; Table[T[n - k - 1, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Robert G. Wilson v, Jul 12 2009 *)
  • Python
    def A139600Row(n):
        x, y = 1, 1
        yield 0
        while True:
            yield x
            x, y = x + y + n, y + n
    for n in range(8):
        R = A139600Row(n)
        print([next(R) for  in range(11)]) # _Peter Luschny, Aug 04 2019
    
  • SageMath
    def T(n,k): return k*(n*(k-1)+2)/2
    def A139600(n,k): return T(n-k, k)
    flatten([[A139600(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Jul 12 2024

Formula

T(n,k) = n*(k-1)*k/2+k.
T(n,k) = A057145(n+2,k). - R. J. Mathar, Jul 28 2016
From Stefano Spezia, Apr 12 2024: (Start)
G.f.: y*(1 - x - y + 2*x*y)/((1 - x)^2*(1 - y)^3).
E.g.f.: exp(x+y)*y*(2 + x*y)/2. (End)

Extensions

Edited by Omar E. Pol, Jan 05 2009

A032091 Number of reversible strings with n-1 beads of 2 colors. 4 beads are black. String is not palindromic.

Original entry on oeis.org

2, 6, 16, 32, 60, 100, 160, 240, 350, 490, 672, 896, 1176, 1512, 1920, 2400, 2970, 3630, 4400, 5280, 6292, 7436, 8736, 10192, 11830, 13650, 15680, 17920, 20400, 23120, 26112, 29376, 32946, 36822, 41040, 45600, 50540, 55860, 61600, 67760, 74382, 81466, 89056
Offset: 6

Views

Author

Keywords

Comments

Also, number of 4-element subsets of the set {1,...,n-1} whose elements sum up to an odd integer, i.e., 4th column of the triangle A159916, cf. there. - M. F. Hasler, May 01 2009
Also, if the offset is changed to 3, so that a(3)=2, a(n) = number of non-equivalent (mod D_3) ways to place 2 indistinguishable points on a triangular grid of side n so that they are not adjacent. - Heinrich Ludwig, Mar 23 2014
Also, the number of binary strings of length n with exactly one pair of consecutive 0s and exactly three pairs of consecutive 1s. - Jeremy Dover, Jul 07 2016
From Petros Hadjicostas, May 19 2018: (Start)
Let k be an integer >= 2. The g.f. of the BHK[k] transform of the sequence (c(n): n >= 1), with g.f. C(x) = Sum_{n>=1} c(n)*x^n, is A_k(x) = (C(x)^k - C(x^2)^(k/2))/2 if k is even, and A_k(x) = (C(x)/2)*(C(x)^{k-1} - C(x^2)^{(k-1)/2}) if k is odd. This follows easily from the formulae in C. G. Bower's web link below about transforms.
When k is odd and c(n) = 1 for all n >= 1, we get C(x) = x/(1-x) and A_k(x) = (1/2)*(x/(1-x))*((x/(1-x))^{k-1} - (x^2/(1-x^2))^{(k-1)/2}). If (a_k(n): n >= 1) is the output sequence (with g.f. A_k(x)), then it can be proved (using Taylor expansions) that a_k(n) = (1/2)*(binomial(n-1, n-k) - binomial(floor((n-1)/2), floor((n-k)/2))) for n >= k+1. (Clearly, a_k(1) = ... = a_k(k) = 0.)
In this sequence, k = 5, and (according to C. G. Bower) a(n) = a_{k=5}(n) is the number of reversible non-palindromic compositions of n with 5 positive parts. If n = b_1 + b_2 + b_3 + b_4 + b_5 is such a composition of n (with b_i >= 1), then it is equivalent to the composition n = b_5 + b_4 + b_3 + b_2 + b_1, and each equivalent class has two elements because here linear palindromes are not allowed as compositions of n.
The fact that we are finding the BHK[5] transform of 1, 1, 1, ... means that each part of each composition of n can have exactly one color (see Bower's link below about transforms).
In each such composition replace each b_i with one black (B) ball followed by b_i - 1 white (W) balls. Then drop the first black (B) ball. We then get a reversible non-palindromic string of length n-1 that has 4 black balls and n-5 white balls. This process, applied to the equivalent compositions n = b_1 + b_2 + b_3 + b_4 + b_5 = b_5 + b_4 + b_3 + b_2 + b_1, gives two strings of length n-1 with 4 black balls and n-5 white balls that are mirror images of each other.
Hence, for n >= 2, a(n) = a_{k=5}(n) is also the number of reversible non-palindromic strings of length n-1 that have k-1 = 4 black balls and n-k = n-5 white balls. (Clearly, a(n) = a_{k=5}(n) > 0 only for n >= 6. For n=5, the composition 1+1+1+1+1, which corresponds to string BBBB, is discarded because it is palindromic.)
For k = 3 (an odd integer) we have a_k(n) = A002620(n-2) (for n >= 4), while for k = 7 (also an odd integer), we have a_k(n) = A032093(n) (for n >= 8).
For k = 4 (which is even), we have a_k(n) = A006584(n-2) (for n >= 5), while for k = 6 and k = 8 (which are also even integers), we get sequences A032092 and A032094, respectively. When k is even, the g.f. in these cases is A_k(x) = (C(x)^k - C(x^2)^(k/2))/2, where C(x) = x/(1-x). The formula for a_k(n) (given above) needs to be modified as well.
(End)
The formula for a(n) for this sequence was Ralf Stephan's conjecture 73. It was solved by Elizabeth Wilmer (see Proposition 2 in one of the links below). There is a minor typo in the original conjecture. - Petros Hadjicostas, Jul 04 2018

Examples

			From _Petros Hadjicostas_, May 19 2018: (Start)
For n=6, we have the following reversible non-palindromic compositions with 5 parts of n: 1+1+1+1+2 (= 2+1+1+1+1) and 1+1+1+2+1 (= 1+2+1+1+1). Using the process described in the comments, we get the following reversible non-palindromic strings with 4 black balls and n-5=1 white balls: BBBBW (= WBBBB) and BBBWB (= BWBBB).
For n=7, we get the following 6 compositions and 6 corresponding strings:
1+1+1+1+3 <-> BBBBWW
1+1+1+3+1 <-> BBBWWB
1+1+1+2+2 <-> BBBWBW
1+1+2+1+2 <-> BBWBBW
1+1+2+2+1 <-> BBWBWB
1+2+1+1+2 <-> BWBBBW
(End)
		

Crossrefs

a(n+6) = 2*A002624(n).
Fourth column of A274228. - Jeremy Dover, Jul 07 2016

Programs

  • Mathematica
    Table[If[EvenQ[n],(n^4-10n^3+32n^2-32n)/48,(n^4-10n^3+32n^2-38n+15)/48], {n,6,50}] (* or *)
    LinearRecurrence[{3,-1,-5,5,1,-3,1},{2,6,16,32,60,100,160},50] (* Harvey P. Dale, Apr 11 2016 *)
    CoefficientList[Series[-2/((x - 1)^5 (x + 1)^2), {x, 0, 42}], x] (* Robert G. Wilson v, Jun 20 2018 *)
  • PARI
    A032091(n)=polcoeff(2/(1-x)^5/(1+x)^2+O(x^(n-5)),n-6)
    A032091(n)=((n-5)*(n-3)*(n-1)^2+if(n%2==0,6*n-15))/48 \\ M. F. Hasler, May 01 2009

Formula

"BHK[ 5 ]" (reversible, identity, unlabeled, 5 parts) transform of 1, 1, 1, 1, ...
From M. F. Hasler, May 01 2009: (Start)
G.f.: -2*x^6 / ((x-1)^5*(x+1)^2). [corrected by Colin Barker, Mar 07 2015]
a(n) = [(n-5)(n-3)(n-1)^2 + (6n-15) X[2Z](n)]/48, where X[2Z] is the characteristic function of 2Z.
(End)
From Colin Barker, Mar 07 2015: (Start)
a(n) = (n^4-10*n^3+32*n^2-32*n)/48 if n is even.
a(n) = (n^4-10*n^3+32*n^2-38*n+15)/48 if n is odd.
(End)
a(n) = (2*n^4 - 20*n^3 + 64*n^2 + 6*(-1)^n*n - 70*n - 15*(-1)^n + 15)/96. - Ilya Gutkovskiy, Jul 08 2016
From Petros Hadjicostas, May 19 2018: (Start)
a(n) = (1/2)*(binomial(n-1, n-5) - binomial(floor((n-1)/2) - floor((n-5)/2))).
G.f.: (1/2)*(x/(1-x))*((x/(1-x))^4 - (x^2/(1-x^2))^2).
(End)
a(n) = 2*A002624(n-6) - Robert G. Wilson v, Jun 20 2018

A331461 Array read by antidiagonals: A(n,k) is the number of nonequivalent binary matrices with k columns and any number of nonzero rows with n ones in every column up to permutation of rows and columns.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 3, 1, 1, 1, 5, 8, 4, 1, 1, 1, 7, 23, 16, 5, 1, 1, 1, 11, 66, 93, 30, 6, 1, 1, 1, 15, 212, 652, 332, 50, 7, 1, 1, 1, 22, 686, 6369, 6414, 1062, 80, 8, 1, 1, 1, 30, 2389, 79568, 226041, 56712, 3117, 120, 9, 1, 1, 1, 42, 8682, 1256425, 12848128, 7295812, 441881, 8399, 175, 10, 1, 1
Offset: 0

Views

Author

Andrew Howroyd, Jan 18 2020

Keywords

Comments

A(n,k) is the number of non-isomorphic set multipartitions (multiset of sets) with k parts each part has size n.

Examples

			Array begins:
===========================================================
n\k | 0 1 2   3    4       5          6              7
----+-----------------------------------------------------
  0 | 1 1 1   1    1       1          1              1 ...
  1 | 1 1 2   3    5       7         11             15 ...
  2 | 1 1 3   8   23      66        212            686 ...
  3 | 1 1 4  16   93     652       6369          79568 ...
  4 | 1 1 5  30  332    6414     226041       12848128 ...
  5 | 1 1 6  50 1062   56712    7295812     1817321457 ...
  6 | 1 1 7  80 3117  441881  195486906   200065951078 ...
  7 | 1 1 8 120 8399 3006771 4298181107 17131523059493 ...
  ...
The A(2,3) = 8 matrices are:
  [1 0 0]  [1 1 0]  [1 1 1]  [1 1 0]  [1 1 0]  [1 1 1]  [1 1 0]  [1 1 1]
  [1 0 0]  [1 0 0]  [1 0 0]  [1 1 0]  [1 0 1]  [1 1 0]  [1 0 1]  [1 1 1]
  [0 1 0]  [0 1 0]  [0 1 0]  [0 0 1]  [0 1 0]  [0 0 1]  [0 1 1]
  [0 1 0]  [0 0 1]  [0 0 1]  [0 0 1]  [0 0 1]
  [0 0 1]  [0 0 1]
  [0 0 1]
		

Crossrefs

Programs

  • PARI
    \\ See A304942 for Blocks
    T(n,k)={Blocks(k, n*k, n)}
    { for(n=0, 7, for(k=0, 6, print1(T(n,k), ", ")); print) }

Formula

A306018(n) = Sum_{d|n} A(n/d, d).

A060099 G.f.: 1/((1-x^2)^3*(1-x)^4).

Original entry on oeis.org

1, 4, 13, 32, 71, 140, 259, 448, 742, 1176, 1806, 2688, 3906, 5544, 7722, 10560, 14223, 18876, 24739, 32032, 41041, 52052, 65429, 81536, 100828, 123760, 150892, 182784, 220116, 263568, 313956, 372096, 438957, 515508, 602889, 702240, 814891, 942172, 1085623
Offset: 0

Views

Author

Wolfdieter Lang, Apr 06 2001

Keywords

Comments

Fourth column (m=3) of triangle A060098.
Partial sums of A038163.
Equals the tetrahedral numbers, [1, 4, 10, 20, ...] convolved with the aerated triangular numbers, [1, 0, 3, 0, 6, 0, 10, ...]. [Gary W. Adamson, Jun 11 2009]

References

  • B. Broer, Hilbert series for modules of covariants, in Algebraic Groups and Their Generalizations..., Proc. Sympos. Pure Math., 56 (1994), Part I, 321-331. See p. 329.

Crossrefs

Cf. A001752 (for the similar series 1/((1-x)^4*(1-x^2))).
Cf. A028346 (for the similar series 1/((1-x)^4*(1-x^2)^2)).

Programs

  • Mathematica
    a[n_]:=If[OddQ[n],((1+n) (3+n) (5+n)^2 (7+n) (9+n))/5760,((2+n) (4+n) (6+n) (8+n) (15+10 n+n^2))/5760]; Map[a,Range[0,100]] (* Peter J. C. Moses, Mar 24 2013 *)
    CoefficientList[Series[1/((1-x^2)^3*(1-x)^4),{x,0,100}],x] (* Peter J. C. Moses, Mar 24 2013 *)
    LinearRecurrence[{4,-3,-8,14,0,-14,8,3,-4,1},{1,4,13,32,71,140,259,448,742,1176},40] (* Harvey P. Dale, Apr 06 2018 *)

Formula

a(n) = Sum_{} A060098(n+3, 3).
G.f.: 1/((1-x)^7*(1+x)^3).

A060098 Triangle of partial sums of column sequences of triangle A060086, read by rows.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 6, 8, 4, 1, 1, 9, 16, 13, 5, 1, 1, 12, 30, 32, 19, 6, 1, 1, 16, 50, 71, 55, 26, 7, 1, 1, 20, 80, 140, 140, 86, 34, 8, 1, 1, 25, 120, 259, 316, 246, 126, 43, 9, 1, 1, 30, 175, 448, 660, 622, 399, 176, 53, 10, 1
Offset: 0

Views

Author

Wolfdieter Lang, Apr 06 2001

Keywords

Comments

In the language of the Shapiro et al. reference (given in A053121) such a lower triangular (ordinary) convolution array, considered as a matrix, belongs to the Riordan-group. The g.f. for the row polynomials p(n,x) = Sum_{m=0..n} a(n,m)*x^m is (1/(1-x*z/((1-z^2)*(1-z))))/(1-z).
Row sums give A052534. Column sequences (without leading zeros) give A000012 (powers of 1), A002620(n+1), A002624, A060099-A060101 for m=0..5.
The bisections of the column sequences give triangles A060102 and A060556.
Riordan array (1/(1-x), x/((1-x)*(1-x^2))). - Paul Barry, Mar 28 2011

Examples

			p(3,x) = 1 + 4*x + 3*x^2 + x^3.
Triangle begins:
  1;
  1,  1;
  1,  2,  1;
  1,  4,  3,  1;
  1,  6,  8,  4,  1;
  1,  9, 16, 13,  5,  1;
  1, 12, 30, 32, 19,  6,  1;
  1, 16, 50, 71, 55, 26,  7,  1;
  ...
		

Crossrefs

Programs

  • Maple
    A060098 := proc(n,k) add( binomial(n-2*i,n-2*i-k)*binomial(k+i-1,i), i=0..floor(n/2)) ; end proc:
    seq(seq(A060098(n,k), k=0..n), n=0..12); # R. J. Mathar, Mar 29 2011
    # Recurrence after Philippe Deléham:
    T := proc(n, k) option remember;
    if k < 0 or k > n then 0 elif k = 0 or n = k then 1 else
    T(n-1, k) + T(n-1, k-1) + T(n-2, k) - T(n-3, k) fi end:
    for n from 0 to 9 do seq(T(n, k), k = 0..n) od;  # Peter Luschny, May 07 2023
  • Mathematica
    t[n_, k_] := Sum[ Binomial[n-2*j, n-2*j-k]*Binomial[k+j-1, j], {j, 0, n/2}]; Table[t[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 21 2013 *)

Formula

G.f. for column m >= 0: ((x/((1-x^2)*(1-x)))^m)/(1-x) = x^m/((1+x)^m*(1-x)^(2*m+1)).
Number triangle T(n,k) = Sum_{i=0..floor(n/2)} C(n-2*i,n-2*i-k)*C(k+i-1,i). - Paul Barry, Mar 28 2011
From Philippe Deléham, Apr 20 2023: (Start)
T(n, k) = 0 if k < 0 or if k > n; T(n, k) = 1 if k = 0 or k = n; otherwise:
T(n, k) = T(n-1, k) + T(n-1, k-1) + T(n-2, k) - T(n-3, k).
T(n, k) = A188316(n, k) + A188316(n-1, k). (End)

A164680 Expansion of x/((1-x)^3*(1-x^2)^3*(1-x^3)).

Original entry on oeis.org

1, 3, 9, 20, 42, 78, 139, 231, 372, 573, 861, 1254, 1791, 2499, 3432, 4629, 6162, 8085, 10492, 13455, 17094, 21503, 26832, 33201, 40795, 49764, 60333, 72687, 87096, 103785, 123075, 145236, 170646, 199626, 232617, 269997, 312277, 359898, 413448, 473438
Offset: 1

Views

Author

Alford Arnold, Aug 21 2009

Keywords

Comments

Convolution of A006918 with A001399, or of A002625 with A059841 (A000035 if offsets are respected),
or of A038163 with A022003 or of A057524 with A027656 or of A014125 with the aerated version of A000217,
or of A002624 with A103221, or of A002623 with A008731, or of other combinations of splitting the signature -/3,3,1 into two components.
If we apply the enumeration of Molien series as described in A139672,
this is row 45=9*5 of a table of values related to Molien series, i.e., the
product of the sequence on row 9 (A006918) with the sequence on row 5 (A001399).
This is associated with the root system E6, and can be described using the additive function on the affine E6 diagram:
1
|
2
|
1--2--3--2--1

Examples

			To calculate a(3), we consider the first three terms of A001399 = (1 1 2...)
and the first three terms of A006918 = (1 2 5 ...), to get the convolved a(3) = 1*5+1*2+2*1 = 9.
		

Crossrefs

Cf. A139672 (row 21).
For G2, the corresponding sequence is A001399.
For F4, the corresponding sequence is A115264.
For E7, the corresponding sequence is A210068.
For E8, the corresponding sequence is A045513.
See A210634 for a closely related sequence.

Programs

  • Magma
    R:=PowerSeriesRing(Integers(), 40); Coefficients(R!( 1/((1-x)^3*(1-x^2)^3*(1-x^3)) )); // G. C. Greubel, Jan 13 2020
  • Maple
    seq(coeff(series(x/((1-x)^3*(1-x^2)^3*(1-x^3)), x, n+1), x, n), n = 1..40); # G. C. Greubel, Jan 13 2020
  • Mathematica
    Rest@CoefficientList[Series[x/((1-x)^3*(1-x^2)^3*(1-x^3)), {x,0,40}], x] (* G. C. Greubel, Jan 13 2020 *)
    LinearRecurrence[{3,0,-7,3,6,0,-6,-3,7,0,-3,1},{1,3,9,20,42,78,139,231,372,573,861,1254},40] (* Harvey P. Dale, Aug 03 2025 *)
  • PARI
    Vec(1/(1-x)^3/(1-x^2)^3/(1-x^3)+O(x^40)) \\ Charles R Greathouse IV, Sep 23 2012
    
  • Sage
    x=PowerSeriesRing(QQ, 'x', 40).gen()
    1/((1-x)^3*(1-x^2)^3*(1-x^3))
    

Formula

a(n) = round( -(-1)^n*(n+3)*(n+7)/256 +(6*n^6 +180*n^5 +2070*n^4 +11400*n^3 +30429*n^2 +34290*n +9785)/103680 ) - R. J. Mathar, Mar 19 2012

Extensions

Edited and extended by R. J. Mathar, Aug 22 2009
Corrected link to index entries - R. J. Mathar, Aug 26 2009

A007009 Number of 3-voter voting schemes with n linearly ranked choices.

Original entry on oeis.org

1, 4, 12, 27, 54, 96, 160, 250, 375, 540, 756, 1029, 1372, 1792, 2304, 2916, 3645, 4500, 5500, 6655, 7986, 9504, 11232, 13182, 15379, 17836, 20580, 23625, 27000, 30720, 34816, 39304, 44217, 49572, 55404, 61731, 68590, 76000, 84000, 92610, 101871, 111804
Offset: 1

Views

Author

Keywords

Comments

With a(0) = 0 nontrivial integer solutions of (x + y)^3 = (x - y)^4. If x = a(n) then y = a(n + (-1)^n). - Thomas Scheuerle, Mar 22 2023

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A034828 (first differences).

Programs

  • Magma
    I:=[1,4,12,27,54,96,160]; [n le 7 select I[n] else 3*Self(n-1)-Self(n-2)- 5*Self(n-3)+5*Self(n-4)+Self(n-5)-3*Self(n-6)+Self(n-7): n in [1..50]]; // Vincenzo Librandi, Sep 21 2015
    
  • Maple
    a:= n-> (Matrix([[0$4, 1, 4, 12, 27]]). Matrix(8, (i, j)-> `if`(i=j-1, 1, `if`(j=1, [4, -4, -4, 10, -4, -4, 4, -1][i], 0)))^n)[1, 1]:
    seq(a(n), n=1..40);  # Alois P. Heinz, Aug 13 2008
  • Mathematica
    LinearRecurrence[{3, -1, -5, 5, 1, -3, 1}, {1, 4, 12, 27, 54, 96, 160}, 50] (* Vincenzo Librandi, Sep 21 2015 *)
  • PARI
    Vec(x*(1-x^3)/((1-x)^4*(1-x^2)^2) + O(x^100)) \\ Colin Barker, Jan 07 2016

Formula

G.f.: x*(1-x^3)/((1-x)^4*(1-x^2)^2) = x*(1+x+x^2)/((1-x)^5*(1+x)^2).
a(n) = (1/2)*Sum_{k=1..n+1} k*floor(k/2)*ceiling(k/2). - Vladeta Jovovic, Apr 29 2006
a(n) = A006009(n)/4.
a(n) = A007590(n+2)*A007590(n+1)/8. - Richard R. Forberg, Dec 03 2013
For n > 1, a(n) = A000332(n+3) - A002624(n-2). - Antal Pinter, Sep 20 2015
a(n) = (n^4 + 6*n^3 + 12*n^2 + 8*n)/32 for n even; a(n) = (n^4 + 6*n^3 + 12*n^2 + 10*n + 3)/32 for n odd. - Colin Barker, Jan 07 2016

Extensions

More terms from James Sellers, Sep 08 2000

A028346 Expansion of 1/((1-x)^4*(1-x^2)^2).

Original entry on oeis.org

1, 4, 12, 28, 58, 108, 188, 308, 483, 728, 1064, 1512, 2100, 2856, 3816, 5016, 6501, 8316, 10516, 13156, 16302, 20020, 24388, 29484, 35399, 42224, 50064, 59024, 69224, 80784, 93840, 108528, 125001, 143412, 163932, 186732, 212002, 239932, 270732, 304612, 341803
Offset: 0

Views

Author

Keywords

Comments

Equals triangle A152205 as an infinite lower triangular matrix * the triangular numbers: [1, 3, 6, ...]. - Gary W. Adamson, Feb 14 2010
a(n) is the number of partitions of n into four kinds of parts 1 and two kinds of parts 2. - Joerg Arndt, Mar 09 2016

Crossrefs

Cf. A152205, A001752 (for the similar series 1/((1-x)^4*(1-x^2))).

Programs

  • Magma
    [(n+4)*(2*n^4+32*n^3+172*n^2+352*n+15*(-1)^n+225)/960: n in [0..40]]; // Vincenzo Librandi, Feb 14 2016
  • Mathematica
    CoefficientList[Series[1/((1 - x)^4 (1 - x^2)^2), {x, 0, 40}], x] (* Vincenzo Librandi, Feb 14 2016 *)
    LinearRecurrence[{4, -4, -4, 10, -4, -4, 4, -1}, {1, 4, 12, 28, 58, 108, 188, 308}, 100] (* G. C. Greubel, Nov 25 2016 *)
  • PARI
    Vec(1/((1-x)^4*(1-x^2)^2)+O(x^99)) \\ Charles R Greathouse IV, Sep 26 2012
    

Formula

a(n) = (n+4)*(2*n^4 + 32*n^3 + 172*n^2 + 352*n + 15*(-1)^n + 225)/960. - R. J. Mathar, Apr 01 2010
From Antal Pinter, Jan 08 2016: (Start)
a(n) = C(n + 3, 3) + 2*C(n + 1, 3) + 3*C(n - 1, 3) + 4*C(n - 3, 3) + ...
a(n) = Sum_{i = 1..z} i*C(n + 5 - 2*i,3) where z = (2*n + 3 + (-1)^n)/4.
(End)
a(n) = Sum_{i = 0..n} A002624(i). - Antal Pinter, May 05 2016

A139672 Convolution of A008619 and A001400.

Original entry on oeis.org

1, 2, 5, 9, 17, 27, 44, 65, 97, 136, 191, 257, 346, 451, 587, 746, 946, 1177, 1461, 1786, 2178, 2623, 3151, 3746, 4443, 5223, 6126, 7131, 8283, 9558, 11007, 12603, 14403, 16377, 18588, 21003, 23692, 26618, 29858, 33372, 37244, 41430, 46022, 50972
Offset: 1

Views

Author

Alford Arnold, Apr 29 2008, May 01 2008

Keywords

Comments

This is row 21 of a table of values related to Molien series. It is the product of the sequence on row 3 (A008619) with the sequence on row 7 (A001400).
This table may be constructed by moving the rows of table A008284 to prime locations and generating the composite locations by multiplication in a manner similar to the calculation illustrated in the present sequence.
Rows 1 thru 20 and 22 thru 25 are as follows:

Programs

  • Maple
    a:= proc(n) local m, r; m:= iquo (n, 12, 'r'); r:= r+1; (19+ (145+ (260+ 15* (r+9)*r+ (405+ 90*r+ 216*m) *m) *m) *m) *m/5+ [0, 1, 2, 5, 9, 17, 27, 44, 65, 97, 136, 191][r]+ [0, 16, 37, 77, 128, 208, 307, 447, 616, 840, 1105, 1441][r]*m/2+ [0, 52, 119, 213, 328, 476, 651, 865, 1112, 1404, 1735, 2117][r]*m^2/2 end: seq (a(n), n=1..50); # Alois P. Heinz, Nov 10 2008
  • Mathematica
    CoefficientList[Series[x/((x^2+x+1)(x^2+1)(x+1)^3 (x-1)^6),{x,0,50}],x] (* or *) LinearRecurrence[{2,1,-3,0,-1,2,2,-1,0,-3,1,2,-1},{0,1,2,5,9,17,27,44,65,97,136,191,257},50] (* Harvey P. Dale, Feb 17 2016 *)

Formula

G.f.: x/((x^2+x+1)*(x^2+1)*(x+1)^3*(x-1)^6). - Alois P. Heinz, Nov 10 2008
a(n)= -A049347(n)/27 +(2*n+11)*(6*n^4+132*n^3+914*n^2+2068*n+1055)/69120 -(-1)^n*(51/512+n^2/256+11*n/256+A057077(n)/32 ). - R. J. Mathar, Nov 21 2008

Extensions

More terms from Alois P. Heinz, Nov 10 2008
Corrected A-number in definition. Added formula. - R. J. Mathar, Nov 21 2008
Showing 1-10 of 15 results. Next