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.

Previous Showing 11-20 of 21 results. Next

A051232 9-factorial numbers.

Original entry on oeis.org

1, 9, 162, 4374, 157464, 7085880, 382637520, 24106163760, 1735643790720, 140587147048320, 12652843234348800, 1252631480200531200, 135284199861657369600, 15828251383813912243200, 1994359674360552942643200, 269238556038674647256832000
Offset: 0

Views

Author

Keywords

Comments

For n >= 1, a(n) is the order of the wreath product of the symmetric group S_n and the Abelian group (C_9)^n. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 07 2001
a(n) = 9*A035023(n) = Product_{k=1..n} 9*k, n >= 1; a(0) := 1.
Pi^n/a(n) is the volume of a 2n-dimensional ball with radius 1/3. - Peter Luschny, Jul 24 2012

Crossrefs

Cf. A047058, A051188, A051189. a(n) = A051231(n-1, 0), A053116 (first column of triangle).

Programs

  • Magma
    [9^n*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Oct 05 2011
  • Maple
    with(combstruct):A:=[N,{N=Cycle(Union(Z$9))},labeled]: seq(count(A,size=n+1)/9, n=0..14); # Zerinvary Lajos, Dec 05 2007
  • Mathematica
    s=1;lst={s};Do[s+=n*s;AppendTo[lst, s], {n, 8, 2*5!, 9}];lst (* Vladimir Joseph Stephan Orlovsky, Nov 08 2008 *)

Formula

a(n) = n!*9^n =: (9*n)(!^9).
E.g.f.: 1/(1-9*x).
G.f.: 1/(1 - 9*x/(1 - 9*x/(1 - 18*x/(1 - 18*x/(1 - 27*x/(1 - 27*x/(1 - ...))))))), a continued fraction. - Ilya Gutkovskiy, Aug 09 2017
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = e^(1/9).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/9). (End)

A257625 Triangle read by rows: T(n, k) = t(n-k, k), where t(0,0) = 1, t(n,m) = 0 if n < 0 or m < 0, else t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), and f(n) = 6*n + 3.

Original entry on oeis.org

1, 3, 3, 9, 54, 9, 27, 621, 621, 27, 81, 6156, 18630, 6156, 81, 243, 57591, 408726, 408726, 57591, 243, 729, 526338, 7685847, 17166492, 7685847, 526338, 729, 2187, 4765473, 132656859, 568014201, 568014201, 132656859, 4765473, 2187
Offset: 0

Views

Author

Dale Gerdemann, May 10 2015

Keywords

Examples

			Array t(n,k) begins as:
    1,       3,          9,           27,             81, ...;
    3,      54,        621,         6156,          57591, ...;
    9,     621,      18630,       408726,        7685847, ...;
   27,    6156,     408726,     17166492,      568014201, ...;
   81,   57591,    7685847,    568014201,    30672766854, ...;
  243,  526338,  132656859,  16305974568,  1366261865802, ...;
  729, 4765473, 2175706332, 427278012876, 53552912878818, ...;
Triangle T(n,k) begins as:
     1;
     3,       3;
     9,      54,         9;
    27,     621,       621,        27;
    81,    6156,     18630,      6156,        81;
   243,   57591,    408726,    408726,     57591,       243;
   729,  526338,   7685847,  17166492,   7685847,    526338,     729;
  2187, 4765473, 132656859, 568014201, 568014201, 132656859, 4765473, 2187;
		

Crossrefs

Cf. A047058 (row sums), A142461, A257616.
See similar sequences listed in A256890.

Programs

  • Mathematica
    t[n_, k_, p_, q_]:= t[n, k, p, q] = If[n<0 || k<0, 0, If[n==0 && k==0, 1, (p*k+q)*t[n-1,k,p,q] + (p*n+q)*t[n,k-1,p,q]]];
    T[n_, k_, p_, q_]= t[n-k, k, p, q];
    Table[T[n,k,6,3], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Mar 01 2022 *)
  • Sage
    @CachedFunction
    def t(n,k,p,q):
        if (n<0 or k<0): return 0
        elif (n==0 and k==0): return 1
        else: return (p*k+q)*t(n-1,k,p,q) + (p*n+q)*t(n,k-1,p,q)
    def A257625(n,k): return t(n-k,k,6,3)
    flatten([[A257625(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Mar 01 2022

Formula

T(n, k) = t(n-k, k), where t(0,0) = 1, t(n,m) = 0 if n < 0 or m < 0, else t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), and f(n) = 6*n + 3.
Sum_{k=0..n} T(n, k) = A047058(n).
From G. C. Greubel, Mar 01 2022: (Start)
t(k, n) = t(n, k).
T(n, n-k) = T(n, k).
t(0, n) = T(n, 0) = A000244(n). (End)

A053101 a(n) = ((6*n+8)(!^6))/8(!^6), related to A034689 (((6*n+2)(!^6))/2 sextic, or 6-factorials).

Original entry on oeis.org

1, 14, 280, 7280, 232960, 8852480, 389509120, 19475456000, 1090625536000, 67618783232000, 4598077259776000, 340257717223424000, 27220617377873920000, 2340973094497157120000, 215369524693738455040000
Offset: 0

Views

Author

Keywords

Comments

Row m=7 of the array A(7; m,n) := ((6*n+m)(!^6))/m(!^6), m >= 0, n >= 0.

Crossrefs

Cf. A047058, A008542(n+1), A034689(n+1), A034723(n+1), A034724(n+1), A034787(n+1), A034788(n+1), A053100, this sequence, A053102, A053103 (rows m=0..10).

Programs

  • Magma
    m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(1/(1-6*x)^(7/3))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Aug 15 2018
  • Mathematica
    s=1;lst={s};Do[s+=n*s;AppendTo[lst, s], {n, 13, 5!, 6}];lst (* Vladimir Joseph Stephan Orlovsky, Nov 08 2008 *)
    With[{nn = 30}, CoefficientList[Series[1/(1 - 6*x)^(7/3), {x, 0, nn}], x]*Range[0, nn]!] (* G. C. Greubel, Aug 15 2018 *)
  • PARI
    x='x+O('x^30); Vec(serlaplace(1/(1-6*x)^(7/3))) \\ G. C. Greubel, Aug 15 2018
    

Formula

a(n) = ((6*n+8)(!^6))/8(!^6)= A034689(n+2)/8.
E.g.f.: 1/(1-6*x)^(7/3).

A051151 Generalized Stirling number triangle of first kind.

Original entry on oeis.org

1, -6, 1, 72, -18, 1, -1296, 396, -36, 1, 31104, -10800, 1260, -60, 1, -933120, 355104, -48600, 3060, -90, 1, 33592320, -13716864, 2104704, -158760, 6300, -126, 1, -1410877440, 609700608, -102114432, 8772624, -423360, 11592, -168
Offset: 1

Views

Author

Keywords

Comments

a(n,m) = R_n^m(a=0, b=6) in the notation of the given 1961 and 1962 references.
a(n,m) is a Jabotinsky matrix, i.e., the monic row polynomials E(n,x) := Sum_{m=1..n} a(n,m)*x^m = Product_{j=0..n-1} (x-6*j), n >= 1, and E(0,x) := 1 are exponential convolution polynomials (see A039692 for the definition and a Knuth reference).
This is the signed Stirling1 triangle A008275 with diagonal d >= 0 (main diagonal d = 0) scaled with 6^d.

Examples

			Triangle a(n,m) (with rows n >= 1 and columns m = 1..n) begins:
        1;
       -6,      1;
       72,    -18,      1;
    -1296,    396,    -36,    1;
    31104, -10800,   1260,  -60,   1;
  -933120, 355104, -48600, 3060, -90, 1;
   ...
3rd row o.g.f.: E(3,x) = 72*x - 18*x^2 + x^3.
		

Crossrefs

First (m=1) column sequence is: A047058(n-1).
Row sums (signed triangle): A008543(n-1)*(-1)^(n-1).
Row sums (unsigned triangle): A008542(n).
Cf. A008275 (Stirling1 triangle, b=1), A039683 (b=2), A051141 (b=3), A051142 (b=4), A051150 (b=5).

Formula

a(n, m) = a(n-1, m-1) - 6*(n-1)*a(n-1, m), n >= m >= 1; a(n, m) := 0, n < m; a(n, 0) := 0 for n >= 1; a(0, 0) = 1.
E.g.f. for the m-th column of the signed triangle: ((log(1 + 6*x)/6)^m)/m!.
a(n, m) = S1(n, m)*6^(n-m), with S1(n, m) := A008275(n, m) (signed Stirling1 triangle).

Extensions

Various sections edited by Petros Hadjicostas, Jun 08 2020

A053103 a(n) = ((6*n+10)(!^6))/10(!^6), related to A034724 (((6*n+4)(!^6))/4 sextic, or 6-factorials).

Original entry on oeis.org

1, 16, 352, 9856, 335104, 13404160, 616591360, 32062750720, 1859639541760, 119016930672640, 8331185147084800, 633170071178444800, 51919945836632473600, 4568955233623657676800, 429481791960623821619200
Offset: 0

Views

Author

Keywords

Comments

Row m=10 of the array A(7; m,n) := ((6*n+m)(!^6))/m(!^6), m >= 0, n >= 0.

Crossrefs

Cf. A047058, A008542(n+1), A034689(n+1), A034723(n+1), A034724(n+1), A034787(n+1), A034788(n+1), A053100, A053101, A053102, this sequence (rows m=0..10).

Programs

  • Magma
    m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(1/(1-6*x)^(8/3))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Aug 16 2018
  • Mathematica
    s=1;lst={s};Do[s+=n*s;AppendTo[lst, s], {n, 15, 5!, 6}];lst (* Vladimir Joseph Stephan Orlovsky, Nov 08 2008 *)
    With[{nn = 30}, CoefficientList[Series[1/(1 - 6*x)^(16/6), {x, 0, nn}], x]*Range[0, nn]!] (* G. C. Greubel, Aug 16 2018 *)
  • PARI
    x='x+O('x^30); Vec(serlaplace(1/(1-6*x)^(8/3))) \\ G. C. Greubel, Aug 16 2018
    

Formula

a(n) = ((6*n+10)(!^6))/10(!^6) = A034724(n+2)/10.
E.g.f.: 1/(1-6*x)^(8/3).

A051262 10-factorial numbers.

Original entry on oeis.org

1, 10, 200, 6000, 240000, 12000000, 720000000, 50400000000, 4032000000000, 362880000000000, 36288000000000000, 3991680000000000000, 479001600000000000000, 62270208000000000000000
Offset: 0

Views

Author

Keywords

Comments

For n >= 1 a(n) is the order of the wreath product of the symmetric group S_n and the Abelian group (C_10)^n. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 07 2001

Crossrefs

a(n) = A048176(n+1, 0)*(-1)^n (first column of unsigned triangle).

Programs

  • Magma
    [10^n*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Oct 05 2011
  • Maple
    with(combstruct):A:=[N,{N=Cycle(Union(Z$10))},labeled]: seq(count(A,size=n)/10,n=0..14); # Zerinvary Lajos, Dec 05 2007
  • Mathematica
    Array[#!*10^# &, 14, 0] (* Michael De Vlieger, Sep 04 2017 *)

Formula

a(n) = 10*A035279(n) = Product_{k=1..n} 10*k, n >= 1; a(0) := 1.
a(n) = n!*10^n =: (10*n)(!^10);
E.g.f.: 1/(1-10*x).
G.f.: 1/(1 - 10*x/(1 - 10*x/(1 - 20*x/(1 - 20*x/(1 - 30*x/(1 - 30*x/(1 - ...))))))), a continued fraction. - Ilya Gutkovskiy, May 12 2017
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = e^(1/10).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/10). (End)

A053100 a(n) = ((6*n+7)(!^6))/7, related to A008542 ((6*n+1)(!^6) sextic, or 6-factorials).

Original entry on oeis.org

1, 13, 247, 6175, 191425, 7082725, 304557175, 14923301575, 820781586625, 50067676784125, 3354534344536375, 244881007151155375, 19345599564941274625, 1644375963020008343125, 149638212634820759224375
Offset: 0

Views

Author

Keywords

Comments

Row m=7 of the array A(7; m,n) := ((6*n+m)(!^6))/m(!^6), m >= 0, n >= 0.

Crossrefs

Cf. A047058, A008542(n+1), A034689(n+1), A034723(n+1), A034724(n+1), A034787(n+1), A034788(n+1), this sequence, A053101, A053102, A053103 (rows m=0..10).

Programs

  • Magma
    m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(1/(1-6*x)^(13/6))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Aug 15 2018
  • Mathematica
    s=1;lst={s};Do[s+=n*s;AppendTo[lst, s], {n, 12, 5!, 6}];lst (* Vladimir Joseph Stephan Orlovsky, Nov 08 2008 *)
    With[{nn=20},CoefficientList[Series[1/(1-6x)^(13/6),{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Apr 20 2015 *)
  • PARI
    x='x+O('x^30); Vec(serlaplace(1/(1-6*x)^(13/6))) \\ G. C. Greubel, Aug 15 2018
    

Formula

a(n) = ((6*n+7)(!^6))/7(!^6) = A008542(n+2)/7.
E.g.f.: 1/(1-6*x)^(13/6).

A053102 a(n) = ((6*n+9)(!^6))/9(!^6), related to A034723 (((6*n+3)(!^6))/3 sextic, or 6-factorials).

Original entry on oeis.org

1, 15, 315, 8505, 280665, 10945935, 492567075, 25120920825, 1431892487025, 90209226682575, 6224436641097675, 466832748082325625, 37813452594668375625, 3289770375736148679375, 305948644943461827181875
Offset: 0

Views

Author

Keywords

Comments

Row m=9 of the array A(7; m,n) := ((6*n+m)(!^6))/m(!^6), m >= 0, n >= 0.

Crossrefs

Cf. A047058, A008542(n+1), A034689(n+1), A034723(n+1), A034724(n+1), A034787(n+1), A034788(n+1), A053100, A053101, this sequence, A053103 (rows m=0..10).

Programs

  • Magma
    m:=30; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!(1/(1-6*x)^(15/6))); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Aug 15 2018
  • Mathematica
    s=1;lst={s};Do[s+=n*s;AppendTo[lst, s], {n, 14, 5!, 6}];lst (* Vladimir Joseph Stephan Orlovsky, Nov 08 2008 *)
    With[{nn = 30}, CoefficientList[Series[1/(1 - 6*x)^(15/6), {x, 0, nn}], x]*Range[0, nn]!] (* G. C. Greubel, Aug 15 2018 *)
  • PARI
    x='x+O('x^30); Vec(serlaplace(1/(1-6*x)^(15/6))) \\ G. C. Greubel, Aug 15 2018
    

Formula

a(n) = ((6*n+9)(!^6))/9(!^6) = A034723(n+2)/9.
E.g.f.: 1/(1-6*x)^(15/6).

A196347 Triangle T(n, k) read by rows, T(n, k) = n!*binomial(n, k).

Original entry on oeis.org

1, 1, 1, 2, 4, 2, 6, 18, 18, 6, 24, 96, 144, 96, 24, 120, 600, 1200, 1200, 600, 120, 720, 4320, 10800, 14400, 10800, 4320, 720, 5040, 35280, 105840, 176400, 176400, 105840, 35280, 5040, 40320, 322560, 1128960, 2257920, 2822400, 2257920, 1128960, 322560, 40320
Offset: 0

Views

Author

Philippe Deléham, Oct 28 2011

Keywords

Comments

Unsigned version of A021012.
Equal to A136572*A007318.

Examples

			Triangle begins:
    1;
    1,   1;
    2,   4,    2;
    6,  18,   18,    6;
   24,  96,  144,   96,  24;
  120, 600, 1200, 1200, 600, 120;
  ...
		

Crossrefs

Programs

  • Magma
    /* As triangle */ [[Factorial(n)*Binomial(n, k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Sep 28 2015
  • Mathematica
    Table[n!*Binomial[n, j], {n, 0, 30}, {j, 0, n}] (* G. C. Greubel, Sep 27 2015 *)
  • Sage
    factorial(n)*binomial(n,k) # Danny Rorabaugh, Sep 27 2015
    

Formula

T(n,k) is given by (1,1,2,2,3,3,4,4,5,5,6,6,...) DELTA (1,1,2,2,3,3,4,4,5,5,6,6, ...) where DELTA is the operator defined in A084938.
Sum_{k>=0} T(m,k)*T(n,k) = (m+n)!.
T(2n,n) = A122747(n).
Sum_{k>=0} T(n,k)^2 = A010050(n) = (2n)!.
Sum_{k>=0} T(n,k)*x^k = A000007(n), A000142(n), A000165(n), A032031(n), A047053(n), A052562(n), A047058(n), A051188(n), A051189(n), A051232(n), A051262(n), A196258(n), A145448(n) for x = -1,0,1,2,3,4,5,6,7,8,9,10,11 respectively.
The row polynomials have the form (x + 1) o (x + 2) o ... o (x + n), where o denotes the black diamond multiplication operator of Dukes and White. See example E10 in the Bala link. - Peter Bala, Jan 18 2018

Extensions

Name exchanged with a formula by Peter Luschny, Feb 01 2015

A131182 Table T(n,k) = n!*k^n, read by upwards antidiagonals.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 2, 1, 0, 6, 8, 3, 1, 0, 24, 48, 18, 4, 1, 0, 120, 384, 162, 32, 5, 1, 0, 720, 3840, 1944, 384, 50, 6, 1, 0, 5040, 46080, 29160, 6144, 750, 72, 7, 1, 0, 40320, 645120, 524880, 122880, 15000, 1296, 98, 8, 1, 0, 362880, 10321920, 11022480, 2949120, 375000, 31104, 2058, 128, 9, 1
Offset: 0

Views

Author

Philippe Deléham, Sep 25 2007

Keywords

Comments

For k>0, T(n,k) is the n-th moment of the exponential distribution with mean = k. - Geoffrey Critzer, Jan 06 2019
T(n,k) is the minimum value of Product_{i=1..n} Sum_{j=1..k} r_j[i] where each r_j is a permutation of {1..n}. For the maximum value, see A331988. - Chai Wah Wu, Sep 01 2022

Examples

			The (inverted) table begins:
k=0: 1, 0,   0,    0,      0,       0, ... (A000007)
k=1: 1, 1,   2,    6,     24,     120, ... (A000142)
k=2: 1, 2,   8,   48,    384,    3840, ... (A000165)
k=3: 1, 3,  18,  162,   1944,   29160, ... (A032031)
k=4: 1, 4,  32,  384,   6144,  122880, ... (A047053)
k=5: 1, 5,  50,  750,  15000,  375000, ... (A052562)
k=6: 1, 6,  72, 1296,  31104,  933120, ... (A047058)
k=7: 1, 7,  98, 2058,  57624, 2016840, ... (A051188)
k=8: 1, 8, 128, 3072,  98304, 3932160, ... (A051189)
k=9: 1, 9, 162, 4374, 157464, 7085880, ... (A051232)
Main diagonal is 1, 1, 8, 162, 6144, 375000, ... (A061711).
		

Crossrefs

Main diagonal gives A061711.

Programs

  • Maple
    T:= (n,k)-> n!*k^n:
    seq(seq(T(d-k, k), k=0..d), d=0..12);  # Alois P. Heinz, Jan 06 2019
  • Python
    from math import factorial
    def A131182_T(n, k): # compute T(n, k)
        return factorial(n)*k**n # Chai Wah Wu, Sep 01 2022

Formula

From Ilya Gutkovskiy, Aug 11 2017: (Start)
G.f. of column k: 1/(1 - k*x/(1 - k*x/(1 - 2*k*x/(1 - 2*k*x/(1 - 3*k*x/(1 - 3*k*x/(1 - ...))))))), a continued fraction.
E.g.f. of column k: 1/(1 - k*x). (End)
Previous Showing 11-20 of 21 results. Next