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

A014608 a(n) = (4n)!/(24^n).

Original entry on oeis.org

1, 1, 70, 34650, 63063000, 305540235000, 3246670537110000, 66475579247327250000, 2390461829733887910000000, 140810154080474667338550000000, 12868639981414579848070084500000000, 1746930746117010628955362040959500000000
Offset: 0

Views

Author

BjornE (mdeans(AT)algonet.se)

Keywords

Comments

a(n) is also the constant term in product 1 <= i,j <= n, i different from j (1 - x_i/x_j)^4. - Sharon Sela (sharonsela(AT)hotmail.com), Feb 16 2002

References

  • George E. Andrews, Richard Askey and Ranjan Roy, Special Functions, Cambridge University Press, 1998.

Crossrefs

Programs

  • Mathematica
    Table[(4n)!/24^n,{n,0,10}] (* Harvey P. Dale, Oct 15 2015 *)
  • PARI
    a(n)=(4*n)!/24^n;

Formula

From Amiram Eldar, Jan 26 2022: (Start)
Sum_{n>=0} 1/a(n) = (cos(2^(3/4)*3^(1/4)) + cosh(2^(3/4)*3^(1/4)))/2.
Sum_{n>=0} (-1)^n/a(n) = cos(6^(1/4))*cosh(6^(1/4)). (End)

A187783 De Bruijn's triangle, T(m,n) = (m*n)!/(n!^m) read by downward antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 6, 6, 1, 1, 1, 20, 90, 24, 1, 1, 1, 70, 1680, 2520, 120, 1, 1, 1, 252, 34650, 369600, 113400, 720, 1, 1, 1, 924, 756756, 63063000, 168168000, 7484400, 5040, 1
Offset: 0

Views

Author

Robert G. Wilson v, Jan 05 2013

Keywords

Comments

From Tilman Piesk, Oct 28 2014: (Start)
Number of permutations of a multiset that contains m different elements n times. These multisets have the signatures A249543(m,n-1) for m>=1 and n>=2.
In an m-dimensional Pascal tensor (the generalization of a symmetric Pascal matrix) P(x1,...,xn) = (x1+...+xn)!/(x1!*...*xn!), so the main diagonal of an m-dimensional Pascal tensor is D(n) = (m*n)!/(n!^m). These diagonals are the rows of this array (with m>0), which begins like this:
m\n:0 1 2 3 4 5
0: 1 1 1 1 1 1 ... A000012;
1: 1 1 1 1 1 1 ... A000012;
2: 1 2 6 20 70 252 ... A000984;
3: 1 6 90 1680 34650 756756 ... A006480;
4: 1 24 2520 369600 63063000 11732745024 ... A008977;
5: 1 120 113400 168168000 305540235000 623360743125120 ... A008978;
6: 1 720 7484400 137225088000 3246670537110000 88832646059788350720 ... A008979;
with columns: A000142 (n=1), A000680 (n=2), A014606 (n=3), A014608 (n=4), A014609 (n=5).
A089759 is the transpose of this matrix. A034841 is its diagonal. A141906 is its lower triangle. A120666 is the upper triangle of this matrix with indices starting from 1. A248827 are the diagonal sums (or the row sums of the triangle).
(End)

Examples

			T(3,5) = (3*5)!/(5!^3) = 756756 = A014609(3) = A006480(5) is the number of permutations of a multiset that contains 3 different elements 5 times, e.g., {1,1,1,1,1,2,2,2,2,2,3,3,3,3,3}.
		

Crossrefs

Cf. A089759 (transposed), A141906 (subtriangle), A120666 (subtriangle transposed), A060538 (1st row/column removed).
Main diagonal gives: A034841.
Row sums of the triangle: A248827.

Programs

  • Magma
    [Factorial(k*(n-k))/(Factorial(n-k))^k: k in [0..n], n in [0..10]]; // G. C. Greubel, Dec 26 2022
    
  • Mathematica
    T[n_, k_]:= (k*n)!/(n!)^k; Table[T[n, k-n], {k, 9}, {n, 0, k-1}]//Flatten
  • SageMath
    def A187783(n,k): return gamma(k*(n-k)+1)/(factorial(n-k))^k
    flatten([[A187783(n,k) for k in range(n+1)] for n in range(11)]) # G. C. Greubel, Dec 26 2022

Formula

T(m,n) = (m*n)!/(n!)^m.
A060540(m,n) = T(m,n)/m! . - R. J. Mathar, Jun 21 2023

Extensions

Row m=0 prepended by Tilman Piesk, Oct 28 2014

A089759 Table T(n,k), 0<=k, 0<=n, read by antidiagonals, defined by T(n,k) = (k*n)! / (n!)^k.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 6, 6, 1, 1, 1, 24, 90, 20, 1, 1, 1, 120, 2520, 1680, 70, 1, 1, 1, 720, 113400, 369600, 34650, 252, 1, 1, 1, 5040, 7484400, 168168000, 63063000, 756756, 924, 1, 1, 1, 40320, 681080400, 137225088000, 305540235000, 11732745024, 17153136, 3432, 1, 1
Offset: 0

Views

Author

Philippe Deléham, Jan 08 2004; revised Jun 08 2005

Keywords

Comments

T(n,k) is the number of lattice paths from {n}^k to {0}^k using steps that decrement one component by 1. - Alois P. Heinz, May 06 2013

Examples

			Row n=0: 1, 1,   1,      1,           1,               1, ... A000012
Row n=1: 1, 1,   2,      6,          24,             120, ... A000142
Row n=2: 1, 1,   6,     90,        2520,          113400, ... A000680
Row n=3: 1, 1,  20,   1680,      369600,       168168000, ... A014606
Row n=4: 1, 1,  70,  34650,    63063000,    305540235000, ... A014608
Row n=5: 1, 1, 252, 756756, 11732745024, 623360743125120, ... A014609
		

Crossrefs

Cf. A000680, A014606, A014608, A014609, A000984, A187783 (transposed version).
Main diagonal gives A034841.

Programs

  • Maple
    T:= (n, k)-> (k*n)!/(n!)^k:
    seq(seq(T(n, d-n), n=0..d), d=0..10);  # Alois P. Heinz, Aug 16 2012
  • Mathematica
    T[n_, k_] := (k*n)!/(n!)^k; Table[T[n-k, k], {n, 0, 10}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Dec 19 2015 *)

Extensions

Corrected by Alois P. Heinz, Aug 16 2012

A278073 Triangle read by rows, coefficients of the polynomials P(m, n) = Sum_{k=1..n} binomial(m*n, m*k)* P(m, n-k)*z with P(m, 0) = 1 and m = 3.

Original entry on oeis.org

1, 0, 1, 0, 1, 20, 0, 1, 168, 1680, 0, 1, 1364, 55440, 369600, 0, 1, 10920, 1561560, 33633600, 168168000, 0, 1, 87380, 42771456, 2385102720, 34306272000, 137225088000, 0, 1, 699048, 1160164320, 158411809920, 5105916816000, 54752810112000, 182509367040000
Offset: 0

Views

Author

Peter Luschny, Jan 22 2017

Keywords

Examples

			Triangle begins:
[1]
[0, 1]
[0, 1,    20]
[0, 1,   168,    1680]
[0, 1,  1364,   55440,   369600]
[0, 1, 10920, 1561560, 33633600, 168168000]
		

Crossrefs

Cf. A014606 (diagonal), A243664 (row sums), A002115 (alternating row sums), A281479 (central coefficients), A327023 (refinement).
Cf. A097805 (m=0), A131689 (m=1), A241171 (m=2), A278074 (m=4).

Programs

  • Maple
    P := proc(m, n) option remember; if n = 0 then 1 else
    add(binomial(m*n, m*k)*P(m, n-k)*x, k=1..n) fi end:
    for n from 0 to 6 do PolynomialTools:-CoefficientList(P(3,n), x) od;
    # Alternatively:
    A278073_row := proc(n)
    1/(1-t*((1/3)*exp(x)+(2/3)*exp(-(1/2)*x)*cos((1/2)*x*sqrt(3))-1));
    expand(series(%,x,3*n+1)); (3*n)!*coeff(%,x,3*n);
    PolynomialTools:-CoefficientList(%,t) end:
    for n from 0 to 6 do A278073_row(n) od;
  • Mathematica
    With[{m = 3}, Table[Expand[j!*SeriesCoefficient[1/(1 - t*(MittagLefflerE[m, x^m] - 1)), {x, 0, j}]], {j, 0, 21, m}]];
    Function[arg, CoefficientList[arg, t]] /@ % // Flatten
  • Sage
    R = PowerSeriesRing(ZZ, 'x')
    x = R.gen().O(30)
    @cached_function
    def P(m, n):
        if n == 0: return R(1)
        return expand(sum(binomial(m*n, m*k)*P(m, n-k)*x for k in (1..n)))
    def A278073_row(n): return list(P(3, n))
    for n in (0..6): print(A278073_row(n)) # Peter Luschny, Mar 24 2020

Formula

E.g.f.: 1/(1-t*((1/3)*exp(x)+(2/3)*exp(-(1/2)*x)*cos((1/2)*x*sqrt(3))-1)), nonzero terms.

A174266 Irregular triangle read by rows: T(n,k) = Sum_{i=0..k} (-1)^i * binomial(3*n+1,i) * binomial(k+3-i,3)^n, 0 <= k <= 3*(n-1).

Original entry on oeis.org

1, 1, 9, 9, 1, 1, 54, 405, 760, 405, 54, 1, 1, 243, 6750, 49682, 128124, 128124, 49682, 6750, 243, 1, 1, 1008, 83736, 1722320, 12750255, 40241088, 58571184, 40241088, 12750255, 1722320, 83736, 1008, 1, 1, 4077, 922347, 45699447, 789300477, 5904797049, 21475242671, 40396577931, 40396577931, 21475242671, 5904797049, 789300477, 45699447, 922347, 4077, 1
Offset: 1

Views

Author

Roger L. Bagula, Mar 14 2010

Keywords

Comments

From Yahia Kahloune, Jan 30 2014: (Start)
In general, let b(k,e,p) = Sum_{i=0..k} (-1)^i*binomial(e*p+1,i)*binomial(k+e-i,e)^p. Then T(n,k) = b(k,3,n).
With these coefficients we can calculate: Sum_{i=1..n} binomial(i+e-1,e)^p = Sum_{i=0..e*(p-1)} b(i,e,p)*binomial(n+e+i,e*p+1).
For example, A086020(n) = Sum_{i=1..n} binomial(2+i, 3)^2 = T(2,0)*binomial(n+3, 7) + T(2,1)*binomial(n+4,7) + T(2,2)*binomial(n+5,7) + T(2,3)*binomial(n+6,7) = (1/5040)*(20*n^7 + 210*n^6 + 854*n^5 + 1680*n^4 + 1610*n^3 + 630*n^2 + 36*n). (End)
T(n,k) is the number of permutations of 3 indistinguishable copies of 1..n with exactly k descents. A descent is a pair of adjacent elements with the second element less than the first. - Andrew Howroyd, May 06 2020

Examples

			Triangle begins:
  1;
  1,      9,         9,            1;
  1,     54,       405,          760,            405,       54,        1;
  1,    243,      6750,        49682,         128124,   128124,    49682, ... ;
  1,   1008,     83736,      1722320,       12750255, 40241088, 58571184, ... ;
  1,   4077,    922347,     45699447,      789300477, ... ;
  1,  16362,   9639783,   1063783164,    38464072830, ... ;
  1,  65511,  98361900,  23119658500,  1641724670475, ... ;
  1, 262116, 992660346, 484099087156, 64856779908606, ... ;
...
The T(2,1) = 9 permutations of 111222 with 1 descent are: 112221, 112212, 112122, 122211, 122112, 121122, 222111, 221112, 211122. - _Andrew Howroyd_, May 07 2020
		

Crossrefs

Row sums are A014606.
Similar triangles for e=1..6: A173018 (or A008292), A154283, this sequence, A236463, A237202, A237252.

Programs

  • Mathematica
    (* First program *)
    p[n_, x_]:= p[n,x]= (1-x)^(3*n+1)*Sum[(Binomial[k+1, 3])^n*x^k, {k, 0, Infinity}]/x^2;
    Table[CoefficientList[p[x, n], x], {n,10}]//Flatten (* corrected by G. C. Greubel, Mar 26 2022 *)
    (* Second program *)
    T[n_, k_]:= T[n, k]= Sum[(-1)^(k-j+1)*Binomial[3*n+1, k-j+1]*(j*(j^2-1)/2)^n, {j, 0, k+1}]/(3^n);
    Table[T[n, k], {n,10}, {k,3*n-2}]//Flatten (* G. C. Greubel, Mar 26 2022 *)
  • PARI
    T(n,k)={sum(i=0, k, (-1)^i*binomial(3*n+1, i)*binomial(k+3-i, 3)^n)} \\ Andrew Howroyd, May 06 2020
    
  • Sage
    @CachedFunction
    def T(n, k): return (1/3^n)*sum( (-1)^(k-j+1)*binomial(3*n+1, k-j+1)*(j*(j^2-1)/2)^n for j in (0..k+1) )
    flatten([[T(n, k) for k in (1..3*n-2)] for n in (1..10)]) # G. C. Greubel, Mar 26 2022

Formula

T(n,k) = [x^k] (1-x)^(3*n+1)*(Sum_{k>=0} (k*(k+1)*(k-1)/2)^n*x^k)/(3^n*x^2).
T(n,k) = T(n, 3*n-k).
From Yahia Kahloune, Jan 30 2014: (Start)
Sum_{i=1..n} binomial(2+i,3)^p = Sum_{i=0..3*p-3} T(p,i)*binomial(n+3+i,3*p+1).
binomial(n,3)^p = Sum_{i=0..3*p-3} T(p,i)*binomial(n+i,3*p). (End)
From Sergii Voloshyn, Dec 18 2024: (Start)
Let E be the operator (x^2)D*(1/x)*D*(x^2)*D, where D denotes the derivative operator d/dx. Then (1/6^n)* E^n(x^2/(1 - x)^4) = (row n generating polynomial)/(1 - x)^(3*n+4) = Sum_{i=0..k} (-1)^i * binomial(3*n+1,i) * binomial(k+3-i,3)^n.
For example, when n = 3 we have 1/216*E^3(x^2/(1 - x)^4) = x^2 (1 + 243x + 6750x^2 + 49682x^3 + 128124x^4 + 128124x^5 + 49682x^6 + 6750x^7 + 243x^8 + x^9)/(1 - x)^13. (End)

Extensions

Edited by Andrew Howroyd, May 06 2020

A014609 a(n) = (5n)!/(5!^n).

Original entry on oeis.org

1, 1, 252, 756756, 11732745024, 623360743125120, 88832646059788350720, 28837919555681211870935040, 18975581770994682860770223800320, 23183587808948692737291767860055162880, 49120458506088132224064306071170476903628800
Offset: 0

Views

Author

BjornE (mdeans(AT)algonet.se)

Keywords

Comments

a(n) is also the constant term in product 1 <= i,j <= n, i different from j (1 - x_i/x_j)^5. - Sharon Sela (sharonsela(AT)hotmail.com), Feb 18 2002

References

  • G. E. Andrews, R. Askey, and R. Roy, Special Functions, Cambridge University Press, 1998.

Crossrefs

Programs

  • PARI
    a(n)=if(n<0,0,(5*n)!/120^n)

A060538 Square array read by antidiagonals of number of ways of dividing n*k labeled items into n labeled boxes with k items in each box.

Original entry on oeis.org

1, 1, 2, 1, 6, 6, 1, 20, 90, 24, 1, 70, 1680, 2520, 120, 1, 252, 34650, 369600, 113400, 720, 1, 924, 756756, 63063000, 168168000, 7484400, 5040, 1, 3432, 17153136, 11732745024, 305540235000, 137225088000, 681080400, 40320, 1, 12870
Offset: 1

Views

Author

Henry Bottomley, Apr 02 2001

Keywords

Examples

			       1        1        1        1
       2        6       20       70
       6       90     1680    34650
      24     2520   369600 63063000
		

Crossrefs

Subtable of A187783.
Rows include A000012, A000984, A006480, A008977, A008978 etc.
Columns include A000142, A000680, A014606, A014608, A014609 etc.
Main diagonal is A034841.

Programs

  • PARI
    T(n,k)=(n*k)!/k!^n;
    for(n=1, 6, for(k=1, 6, print1(T(n,k), ", ")); print) \\ Harry J. Smith, Jul 06 2009

Formula

T(n, k) = (nk)!/k!^n = T(n-1, k)*binomial(nk, k) = T(n-1, k)*A060539(n, k) = A060540(n, k)*A000142(k).

A299041 Irregular table: T(n,k) equals the number of alignments of length k of n strings each of length 3.

Original entry on oeis.org

1, 1, 12, 30, 20, 1, 60, 690, 2940, 5670, 5040, 1680, 1, 252, 8730, 103820, 581700, 1767360, 3087000, 3099600, 1663200, 369600, 1, 1020, 94890, 2615340, 32186070, 214628400, 859992000, 2189325600, 3628409400, 3903900000, 2630628000, 1009008000, 168168000, 1, 4092, 979530, 58061420, 1411122300
Offset: 1

Views

Author

Peter Bala, Feb 02 2018

Keywords

Comments

An alignment of n strings of various lengths is a way of inserting blank characters into the n strings so that the resulting strings all have the same length. We don't allow insertion of a blank character into the same position in each of the n strings.
In this case, let s_1,...,s_n be n strings each of length 3 over an alphabet A. Let - be a gap symbol not in A and let A' = union of A and {-}. An alignment of the n strings is an n-tuple (s_1',...,s_n') of strings each of length >= 3 over the alphabet A' such that
(a) the strings s_i', 1 <= i <= n, have the same length. This common length is called the length of the alignment.
(b) deleting the gap symbols from s_i' yields the string s_i for 1 <= i <= n
(c) there is no value j such that all the strings s_i', 1 <= i <= n have a gap symbol at position j.
By writing the strings s_i' one under another we can consider an alignment of n strings as an n X L matrix, where L, the length of the alignment, ranges from a minimum value of 3 to a maximum value of 3*n. Each row of the matrix has 3 characters from the alphabet A and (L - 3) gap characters.
For example,
s_1' = ABC------
s_2' = ---DEF---
s_3' = ------GHI
is an alignment (of maximum length L = 9) of three strings s_1 = ABC, s_2 = DEF and s_3 = GHI each of length 3.
For the number of alignments of length k of n strings of length 1 (resp. 2) see A131689 (resp. A122193).

Examples

			Table begins
n\k| 3   4     5       6      7      8        9      10
- - - - - - - - - - - - - - - - - - - - - - - - - - - -
  1| 1
  2| 1  12    30      20
  3| 1  60   690    2940   5670    5040    1680
  4| 1 252  8730  103820 581700 1767360 3087000 3099600 ...
...
T(2,5) = 30: An alignment of length 5 will have two gap symbols on each line. There are C(5,2) = 10 ways of choosing the 2 positions to insert the gap symbols in the first string. The second string in the alignment must then have nongap symbols at these two positions leaving three positions in which to insert the remaining 1 nongap symbol, giving in total 10 x 3 = 30 possible alignments of 2 strings of 3 characters. Some examples are
  ABC--   ABC--   ABC--
  D--EF   -D-EF   --DEF
Row 2: Sum_{i = 3..n-1} C(i,3)^2 = C(n,4) + 12*C(n,5) + 30*C(n,6) + 20*C(n,7).
Row 3: Sum_{i = 3..n-1} C(i,3)^3 = C(n,4) + 60*C(n,5) + 690*C(n,6) + 2940*C(n,7) + 5670*C(n,8)+ 5040*C(n,9)+ 1680*C(n,10).
exp( Sum_{n >= 1} R(n,2)*x^n/n ) = (1 + x + 153*x^2 + 128793*x^3 + 319155321*x^4 + 1744213657689*x^5 + ....)^8
exp( Sum_{n >= 1} R(n,3)*x^n/n ) = (1 + x + 424*x^2 + 998584*x^3 + 6925040260*x^4 + 105920615923684*x^5 + ....)^27.
		

Crossrefs

Programs

  • Maple
    seq(seq(add( (-1)^(k-i) *binomial(k,i)*binomial(i,3)^n, i = 0..k ), k = 3..3*n), n = 1..6);
  • Mathematica
    nmax = 6; T[n_, k_] := Sum[(-1)^(k-i) Binomial[k, i] Binomial[i, 3]^n, {i, 0, k}]; Table[T[n, k], {n, 1, nmax}, {k, 3, 3n}] // Flatten (* Jean-François Alcover, Feb 20 2018 *)

Formula

T(n,k) = Sum_{i = 0..k} (-1)^(k-i)*binomial(k,i)*binomial(i,3)^n.
T(n,3) = 1; T(n,3*n) = (3*n)!/6^n = A014606(n)
T(n,k) = binomial(k,3)*( T(n-1,k) + 3*T(n-1,k-1) + 3*T(n-1,k-2) + T(n-1,k-3) ) for 3 <= k <= 3*n with boundary conditions T(n,3) = 1 for n >= 1 and T(n,k) = 0 if (k < 3) or (k > 3*n).
Double e.g.f.: exp(-x)*Sum_{n >= 0} exp(binomial(n,3)*y)*x^n/n! = 1 + (x^3/3!)*y + (x^3/3! + 12*x^4/4! + 30*x^5/5! + 20*x^6/6!)*y^2/2! + ....
n-th row polynomial R(n,x) = Sum_{i >= 3} binomial(i,3)^n*x^i/(1 + x)^(i+1) for n >= 1.
1/(1 - x)*R(n,x/(1 - x)) = Sum_{i >= 3} binomial(i,3)^n*x^i for n >= 1.
R(n,x) = x^3 o x^3 o ... o x^3 (n factors), where o is the black diamond product of power series defined in Dukes and White.
R(n,x) = coefficient of (z_1)^3*...*(z_n)^3 in the expansion of the rational function 1/(1 + x - x*(1 + z_1)*...*(1 + z_n)).
The polynomials Sum_{k = 3..3*n} T(n,k)*x^(k-3)*(1 - x)^(3*n-k) are the row polynomials of A174266.
Sum_{i = 3..n-1} binomial(i,3)^m = Sum_{k = 3..3*m} T(m,k)*binomial(n,k+1) for m >= 1. See Examples below.
x^3*R(n,-1 - x) = (-1)^n*(1 + x)^3*R(n,x).
R(n+1,x) = 1/3!*x^3*(d/dx)^3 ((1 + x)^3*R(n,x)) for n >= 1.
The zeros of R(n,x) belong to the interval [-1, 0].
Row sums R(n,1) = A062208(n); alternating row sums R(n,-1) = (-1)^n.
For k a nonzero integer, the power series A(k,x) := exp( Sum_{n >= 1} 1/k^3*R(n,k)*x^n/n ) appear to have integer coefficients. See the Example section.
Sum_{k = 3..3*n} T(n,k)*binomial(x,k) = ( binomial(x,3) )^n. Equivalently, Sum_{k = 3..3*n} (-1)^(n+k)*T(n,k)*binomial(x+k,k) = ( binomial(x+3,3) )^n. Cf. the Worpitzky-type identity Sum_{k = 1..n} A019538(n,k)* binomial(x,k) = x^n.
Sum_{k = 3..3*n} T(n,k)*binomial(x,k-3) = -binomial(x,3)^n + 3*binomial(x+1,3)^n - 3*binomial(x+2,3)^n + binomial(x+3,3)^n. These polynomials have their zeros on the vertical line Re x = -1/2 in the complex plane.

A266737 Number of words on {1,1,1,2,2,2,3,3,3,4,4,4,...,n,n,n} avoiding the pattern 1234.

Original entry on oeis.org

1, 1, 20, 1680, 173891, 21347262, 2977892253, 455912368540, 74876841353159, 12990339123973119, 2354973430941967605, 442587722191655715108, 85717352536181708342445, 17029266882947116165470103, 3457866959157770598680361537, 715559803849259851987691458500
Offset: 0

Views

Author

N. J. A. Sloane, Jan 06 2016

Keywords

Crossrefs

Extensions

More terms from Alois P. Heinz, Jan 14 2016

A292604 Triangle read by rows, coefficients of generalized Eulerian polynomials F_{2}(x).

Original entry on oeis.org

1, 1, 0, 5, 1, 0, 61, 28, 1, 0, 1385, 1011, 123, 1, 0, 50521, 50666, 11706, 506, 1, 0, 2702765, 3448901, 1212146, 118546, 2041, 1, 0, 199360981, 308869464, 147485535, 24226000, 1130235, 8184, 1, 0
Offset: 0

Views

Author

Peter Luschny, Sep 20 2017

Keywords

Comments

The generalized Eulerian polynomials F_{m}(x) are defined F_{m; 0}(x) = 1 for all m >= 0 and for n > 0:
F_{0; n}(x) = Sum_{k=0..n} A097805(n, k)*(x-1)^(n-k) with coeffs. in A129186.
F_{1; n}(x) = Sum_{k=0..n} A131689(n, k)*(x-1)^(n-k) with coeffs. in A173018.
F_{2; n}(x) = Sum_{k=0..n} A241171(n, k)*(x-1)^(n-k) with coeffs. in A292604.
F_{3; n}(x) = Sum_{k=0..n} A278073(n, k)*(x-1)^(n-k) with coeffs. in A292605.
F_{4; n}(x) = Sum_{k=0..n} A278074(n, k)*(x-1)^(n-k) with coeffs. in A292606.
The case m = 1 are the Eulerian polynomials whose coefficients are the Eulerian numbers which are displayed in Euler's triangle A173018.
Evaluated at x in {-1, 1, 0} these families of polynomials give for the first few m:
F_{m} : F_{0} F_{1} F_{2} F_{3} F_{4}
x = 1: A000012 A000142 A000680 A014606 A014608 ... (m*n)!/m!^n
x = 0: -- A000012 A000364 A002115 A211212 ... m-alternating permutations of length m*n.
Note that the constant terms of the polynomials are the generalized Euler numbers as defined in A181985. In this sense generalized Euler numbers are also generalized Eulerian numbers.

Examples

			Triangle starts:
[n\k][    0        1        2       3     4  5  6]
--------------------------------------------------
[0][      1]
[1][      1,       0]
[2][      5,       1,       0]
[3][     61,      28,       1,      0]
[4][   1385,    1011,     123,      1,    0]
[5][  50521,   50666,   11706,    506,    1, 0]
[6][2702765, 3448901, 1212146, 118546, 2041, 1, 0]
		

References

  • G. Frobenius. Über die Bernoullischen Zahlen und die Eulerschen Polynome. Sitzungsber. Preuss. Akad. Wiss. Berlin, pages 200-208, 1910.

Crossrefs

F_{0} = A129186, F_{1} = A173018, F_{2} is this triangle, F_{3} = A292605, F_{4} = A292606.
First column: A000364. Row sums: A000680. Alternating row sums: A002105.

Programs

  • Maple
    Coeffs := f -> PolynomialTools:-CoefficientList(expand(f), x):
    A292604_row := proc(n) if n = 0 then return [1] fi;
    add(A241171(n, k)*(x-1)^(n-k), k=0..n); [op(Coeffs(%)), 0] end:
    for n from 0 to 6 do A292604_row(n) od;
  • Mathematica
    T[n_, k_] /; 1 <= k <= n := T[n, k] = k (2 k - 1) T[n - 1, k - 1] + k^2 T[n - 1, k]; T[, 1] = 1; T[, _] = 0;
    F[2, 0][] = 1; F[2, n][x_] := Sum[T[n, k] (x - 1)^(n - k), {k, 0, n}];
    row[n_] := If[n == 0, {1}, Append[CoefficientList[ F[2, n][x], x], 0]];
    Table[row[n], {n, 0, 7}] (* Jean-François Alcover, Jul 06 2019 *)
  • Sage
    def A292604_row(n):
        if n == 0: return [1]
        S = sum(A241171(n, k)*(x-1)^(n-k) for k in (0..n))
        return expand(S).list() + [0]
    for n in (0..6): print(A292604_row(n))

Formula

F_{2; n}(x) = Sum_{k=0..n} A241171(n, k)*(x-1)^(n-k) for n>0 and F_{2; 0}(x) = 1.
Showing 1-10 of 40 results. Next