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

A008546 Quintuple factorial numbers: Product_{k = 0..n-1} (5*k + 4).

Original entry on oeis.org

1, 4, 36, 504, 9576, 229824, 6664896, 226606464, 8837652096, 388856692224, 19053977918976, 1028914807624704, 60705973649857536, 3885182313590882304, 268077579637770878976, 19837740893195045044224, 1567181530562408558493696, 131643248567242318913470464
Offset: 0

Views

Author

Joe Keane (jgk(AT)jgk.org)

Keywords

Crossrefs

Programs

  • GAP
    List([0..20], n-> Product([0..n-1], k-> 5*k+4 )); # G. C. Greubel, Aug 20 2019
  • Magma
    [1] cat [(&*[5*k+4: k in [0..n-1]]): n in [1..20]]; // G. C. Greubel, Aug 20 2019
    
  • Maple
    f:= n-> product(5*k+4, k=0..n-1);
  • Mathematica
    FoldList[Times, 1, 5Range[0, 20] + 4] (* Vincenzo Librandi, Jun 10 2013 *)
    CoefficientList[Series[(1 - 5x)^(-4/5), {x, 0, 20}], x] Range[0, 20]! (* Vaclav Kotesovec, Jan 28 2015 *)
    Table[5^n Pochhammer[4/5, n], {n, 0, 20}] (* G. C. Greubel, Aug 20 2019 *)
  • PARI
    vector(20, n, n--; prod(j=0,n-1, 5*j+4) ) \\ G. C. Greubel, Aug 20 2019
    
  • Sage
    [5^n*rising_factorial(4/5, n) for n in (0..20)] # G. C. Greubel, Aug 20 2019
    

Formula

a(n) = 4*A034301(n) = (5*n - 1)(!^5), n >= 1, with a(0) = 1.
a(n) = A011801(n + 1, 1) (first column of triangle).
a(n) ~ (sqrt(2*Pi)/Gamma(4/5))*n^(n + 3/10)*(5/e)^n*(1 + 1/(300*n) + ...). - Joe Keane (jgk(AT)jgk.org), Nov 24 2001
G.f.: 1/(1 - 4*x/(1 - 5*x/(1 - 9*x/(1 - 10*x/(1 - 14*x/(1 - 15*x/(1 - 19*x/(1 - 20*x/(1 - 24*x/(1 - ... (continued fraction). - Philippe Deléham, Jan 08 2012
a(n) = (-1)^n*Sum_{k = 0..n} 5^k*s(n + 1, n + 1 - k), where s(n, k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
G.f.: ( 1 - 1/Q(0) )/x where Q(k) = 1 - x*(5*k - 1)/(1 - x*(5*k + 5)/Q(k + 1) ); (continued fraction); e.g.f. (1 - 5*x)^(-4/5). - Sergei N. Gladkovskii, Mar 20 2013
G.f.: 1/x - G(0)/(2*x), where G(k) = 1 + 1/(1 - x*(5*k - 1)/(x*(5*k - 1) + 1/G(k + 1))); (continued fraction). - Sergei N. Gladkovskii, May 27 2013
a(n) = 5^n * Gamma(n + 4/5) / Gamma(4/5). - Vaclav Kotesovec, Jan 28 2015
a(n) + (-5*n + 1)*a(n - 1) = 0. - R. J. Mathar, Sep 04 2016
G.f.: 1/(1 - 4*x - 20*x^2/(1 - 14*x - 90*x^2/(1 - 24*x - 210*x^2/(1 - 34*x - 380*x^2/(1 - 44*x - 600*x^2/(1 - 54*x - 870*x^2/(1 - ...))))))) (Jacobi continued fraction). - Nikolaos Pantelidis, Feb 29 2020
Sum_{n>=0} 1/a(n) = 1 + (e/5)^(1/5)*(Gamma(4/5) - Gamma(4/5, 1/5)). - Amiram Eldar, Dec 19 2022

A134278 A certain partition array in Abramowitz-Stegun order (A-St order), called M_3(6).

Original entry on oeis.org

1, 6, 1, 66, 18, 1, 1056, 264, 108, 36, 1, 22176, 5280, 3960, 660, 540, 60, 1, 576576, 133056, 95040, 43560, 15840, 23760, 3240, 1320, 1620, 90, 1, 17873856, 4036032, 2794176, 2439360, 465696, 665280, 304920, 249480, 36960, 83160, 22680, 2310
Offset: 1

Views

Author

Wolfdieter Lang, Nov 13 2007

Keywords

Comments

For the A-St order of partitions see the Abramowitz-Stegun reference given in A117506.
Partition number array M_3(6), the k=6 member in the family of a generalization of the multinomial number arrays M_3 = M_3(1) = A036040.
The sequence of row lengths is A000041 (partition numbers) [1, 2, 3, 5, 7, 11, 15, 22, 30, 42, ...].
The S2(6,n,m):=A049385(n,m) numbers (generalized Stirling2 numbers) are obtained by summing in row n all numbers with the same part number m. In the same manner the S2(n,m) (Stirling2) numbers A008277 are obtained from the partition array M_3= A036040.
a(n,k) enumerates unordered forests of increasing 6-ary trees related to the k-th partition of n in the A-St order. The m-forest is composed of m such trees, with m the number of parts of the partition.

Examples

			[1]; [6,1]; [66,18,1]; [1056,264,108,36,1]; [22176,5280,3960,660,540,60,1]; ...
There are a(4,3) = 108 = 3*6^2 unordered 2-forests with 4 vertices, composed of two 6-ary increasing trees, each with two vertices: there are 3 increasing labelings (1,2)(3,4); (1,3)(2,4); (1,4)(2,3) and each tree comes in six versions from the 6-ary structure.
		

Crossrefs

Cf. A049412 (row sums, also of triangle A049385).
Cf. A134273 (M_3(5) partition array).

Formula

a(n,k) = n!*Product_{j=1..n} (S2(6,j,1)/j!)^e(n,k,j)/e(n,k,j)! with S2(6,n,1) = A049385(n,1) = A008548(n) = (5*n-4)(!^5) (quintuple- or 5-factorials) and the exponent e(n,k,j) of j in the k-th partition of n in the A-St ordering of the partitions of n. Exponents 0 can be omitted due to 0!=1.

A045755 8-fold factorials: a(n) = Product_{k=0..n-1} (8*k+1).

Original entry on oeis.org

1, 1, 9, 153, 3825, 126225, 5175225, 253586025, 14454403425, 939536222625, 68586144251625, 5555477684381625, 494437513909964625, 47960438849266568625, 5035846079172989705625, 569050606946547836735625, 68855123440532288245010625, 8882310923828665183606370625
Offset: 0

Views

Author

Keywords

Crossrefs

Cf. k-fold factorials : A000142, A001147, A007559, A007696, A008548, A008542, A045754.

Programs

  • GAP
    List([0..20], n-> Product([0..n-1], j-> 8*j+1) ); # G. C. Greubel, Nov 11 2019
  • Magma
    [1] cat [(&*[8*j+1: j in [0..n-1]]): n in [1..20]]; // G. C. Greubel, Nov 11 2019
    
  • Maple
    a := n->product(8*k+1), k=0..(n-1));
  • Mathematica
    Table[8^n*Pochhammer[1/8, n], {n,0,20}] (* G. C. Greubel, Nov 11 2019 *)
  • PARI
    a(n)=prod(k=0, n, 8*k+1);
    
  • Sage
    [product( (8*j+1) for j in (0..n-1)) for n in (0..20)] # G. C. Greubel, Nov 11 2019
    

Formula

a(n+1) = (8*n+1)(!^8).
a(n) = Sum_{k=0..n} (-8)^(n-k)*A048994(n, k); A048994 = Stirling-1 numbers.
E.g.f.: (1-8*x)^(-1/8).
G.f.: 1+x/(1-9x/(1-8x/(1-17x/(1-16x/(1-25x/(1-24x/(1-33x/(1-32x/(1-... (continued fraction). - Philippe Deléham, Jan 07 2012
a(n) = (-7)^n*Sum_{k=0..n} (8/7)^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. [Mircea Merca, May 03 2012]
G.f.: 1/Q(0) where Q(k) = 1 - x*(8*k+1)/(1 - x*(8*k+8)/Q(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 20 2013
G.f.: 2/G(0), where G(k)= 1 + 1/(1 - 2*x*(8*k+1)/(2*x*(8*k+1) - 1 + 16*x*(k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 30 2013
G.f.: G(0)/2, where G(k)= 1 + 1/(1 - x*(8*k+1)/(x*(8*k+1) + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 05 2013
a(n) = 8^n * Gamma(n + 1/8) / Gamma(1/8). - Artur Jasinski,Aug 23 2016
a(n) ~ sqrt(2*Pi) * 8^n * n^(n - 3/8)/(Gamma(1/8)*exp(n)). - Ilya Gutkovskiy, Sep 10 2016
D-finite with recurrence: a(n) +(-8*n+7)*a(n-1)=0. - R. J. Mathar, Jan 17 2020
Sum_{n>=0} 1/a(n) = 1 + (e/8^7)^(1/8)*(Gamma(1/8) - Gamma(1/8, 1/8)). - Amiram Eldar, Dec 20 2022

Extensions

Additional comments from Philippe Deléham and Paul D. Hanna, Oct 29 2005
Edited by N. J. A. Sloane, Oct 14 2008 at the suggestion of Artur Jasinski.

A047055 Quintuple factorial numbers: a(n) = Product_{k=0..n-1} (5*k + 2).

Original entry on oeis.org

1, 2, 14, 168, 2856, 62832, 1696464, 54286848, 2008613376, 84361761792, 3965002804224, 206180145819648, 11752268311719936, 728640635326636032, 48818922566884614144, 3514962424815692218368, 270652106710808300814336, 22193472750286280666775552
Offset: 0

Views

Author

Joe Keane (jgk(AT)jgk.org)

Keywords

Comments

Hankel transform is A169621. - Paul Barry, Dec 03 2009

Crossrefs

Programs

  • GAP
    List([0..20], n-> Product([0..n-1], k-> (5*k+2) )); # G. C. Greubel, Aug 17 2019
  • Magma
    [1] cat [(&*[(5*k+2): k in [0..n-1]]): n in [1..20]]; // G. C. Greubel, Aug 17 2019
    
  • Maple
    a := n->product(5*i+2,i=0..n-1); [seq(a(j),j=0..30)];
  • Mathematica
    Table[5^n*Pochhammer[2/5, n], {n,0,20}] (* G. C. Greubel, Aug 17 2019 *)
    Join[{1},FoldList[Times,5*Range[0,20]+2]] (* Harvey P. Dale, Apr 03 2025 *)
  • PARI
    vector(20, n, n--; prod(k=0,n-1, 5*k+2)) \\ G. C. Greubel, Aug 17 2019
    
  • Sage
    [product((5*k+2) for k in (0..n-1)) for n in (0..20)] # G. C. Greubel, Aug 17 2019
    

Formula

E.g.f. (1-5*x)^(-2/5).
a(n) ~ sqrt(2*Pi)/Gamma(2/5)*n^(-1/10)*(5n/e)^n*(1 - (11/300)/n - ...). - Joe Keane (jgk(AT)jgk.org), Nov 24 2001
a(n) = A084940(n)/A000142(n)*A000079(n) = 5^n*Pochhammer(2/5, n) = 5^n*Gamma(n+2/5)*sin(2*Pi/5)*Gamma(3/5)/Pi. - Daniel Dockery (peritus(AT)gmail.com), Jun 13 2003
G.f.: 1/(1-2x/(1-5x/(1-7x/(1-10x/(1-12x/(1-15x/(1-17x/(1-20x/(1-22x/(1-25x/(1-.../(1-A047215(n+1)*x/(1-... (continued fraction). - Paul Barry, Dec 03 2009
a(n) = (-3)^n*Sum_{k=0..n} (5/3)^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
D-finite with recurrence: a(n) +(-5*n+3)*a(n-1) = 0. - R. J. Mathar, Dec 03 2012
G.f.: 1/G(0) where G(k) = 1 - x*(5*k+2)/( 1 - 5*x*(k+1)/G(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 23 2013
Sum_{n>=0} 1/a(n) = 1 + (e/5^3)^(1/5)*(Gamma(2/5) - Gamma(2/5, 1/5)). - Amiram Eldar, Dec 19 2022

A052562 a(n) = 5^n * n!.

Original entry on oeis.org

1, 5, 50, 750, 15000, 375000, 11250000, 393750000, 15750000000, 708750000000, 35437500000000, 1949062500000000, 116943750000000000, 7601343750000000000, 532094062500000000000, 39907054687500000000000
Offset: 0

Views

Author

Joe Keane (jgk(AT)jgk.org)

Keywords

Comments

A simple regular expression in a labeled universe.
For n >= 1 a(n) is the order of the wreath product of the symmetric group S_n and the Abelian group (C_5)^n. - Ahmed Fares (ahmedfares(AT)my-deja.com), May 07 2001

Crossrefs

Programs

  • Magma
    [5^n*Factorial(n): n in [0..20]]; // Vincenzo Librandi, Oct 05 2011
    
  • Maple
    spec := [S,{S=Sequence(Union(Z,Z,Z,Z,Z))},labeled]: seq(combstruct[count](spec,size=n), n=0..20);
    with(combstruct):A:=[N,{N=Cycle(Union(Z$5))},labeled]: seq(count(A,size=n)/5,n=1..16); # Zerinvary Lajos, Dec 05 2007
  • Mathematica
    Table[5^n*n!, {n, 0, 20}] (* Wesley Ivan Hurt, Sep 28 2013 *)
  • PARI
    {a(n) = 5^n*n!}; \\ G. C. Greubel, May 05 2019
    
  • Sage
    [5^n*factorial(n) for n in (0..20)] # G. C. Greubel, May 05 2019

Formula

a(n) = A051150(n+1, 0) (first column of triangle).
E.g.f.: 1/(1-5*x).
a(n) = 5*n*a(n-1) with a(0)=1.
G.f.: 1/(1-5*x/(1-5*x/(1-10*x/(1-10*x/(1-15*x/(1-15*x/(1-20*x/(1-... (continued fraction). - Philippe Deléham, Jan 08 2012
G.f.: 1/Q(0), where Q(k) = 1 - 5*x*(2*k+1) - 25*x^2*(k+1)^2/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Sep 28 2013
a(n) = n!*A000351(n). - R. J. Mathar, Aug 21 2014
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = e^(1/5) (A092514).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/5) (A092618). (End)

Extensions

Name changed by Arkadiusz Wesolowski, Oct 04 2011

A256890 Triangle T(n,k) = t(n-k, k); t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.

Original entry on oeis.org

1, 2, 2, 4, 12, 4, 8, 52, 52, 8, 16, 196, 416, 196, 16, 32, 684, 2644, 2644, 684, 32, 64, 2276, 14680, 26440, 14680, 2276, 64, 128, 7340, 74652, 220280, 220280, 74652, 7340, 128, 256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172, 256, 512, 72076, 1637860, 10978444, 27227908, 27227908, 10978444, 1637860, 72076, 512
Offset: 0

Views

Author

Dale Gerdemann, Apr 12 2015

Keywords

Comments

Related triangles may be found by varying the function f(x). If f(x) is a linear function, it can be parameterized as f(x) = a*x + b. With different values for a and b, the following triangles are obtained:
a\b 1.......2.......3.......4.......5.......6
The row sums of these, and similarly constructed number triangles, are shown in the following table:
a\b 1.......2.......3.......4.......5.......6.......7.......8.......9
The formula can be further generalized to: t(n,m) = f(m+s)*t(n-1,m) + f(n-s)*t(n,m-1), where f(x) = a*x + b. The following table specifies triangles with nonzero values for s (given after the slash).
a\b 0 1 2 3
-2 A130595/1
-1
0
With the absolute value, f(x) = |x|, one obtains A038221/3, A038234/4, A038247/5, A038260/6, A038273/7, A038286/8, A038299/9 (with value for s after the slash).
If f(x) = A000045(x) (Fibonacci) and s = 1, the result is A010048 (Fibonomial).
In the notation of Carlitz and Scoville, this is the triangle of generalized Eulerian numbers A(r, s | alpha, beta) with alpha = beta = 2. Also the array A(2,1,4) in the notation of Hwang et al. (see page 31). - Peter Bala, Dec 27 2019

Examples

			Array, t(n, k), begins as:
   1,    2,      4,        8,        16,         32,          64, ...;
   2,   12,     52,      196,       684,       2276,        7340, ...;
   4,   52,    416,     2644,     14680,      74652,      357328, ...;
   8,  196,   2644,    26440,    220280,    1623964,    10978444, ...;
  16,  684,  14680,   220280,   2643360,   27227908,   251195000, ...;
  32, 2276,  74652,  1623964,  27227908,  381190712,  4677894984, ...;
  64, 7340, 357328, 10978444, 251195000, 4677894984, 74846319744, ...;
Triangle, T(n, k), begins as:
    1;
    2,     2;
    4,    12,      4;
    8,    52,     52,       8;
   16,   196,    416,     196,      16;
   32,   684,   2644,    2644,     684,      32;
   64,  2276,  14680,   26440,   14680,    2276,     64;
  128,  7340,  74652,  220280,  220280,   74652,   7340,   128;
  256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172,   256;
		

Crossrefs

Programs

  • Magma
    A256890:= func< n,k | (&+[(-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n: j in [0..k]]) >;
    [A256890(n,k): k in [0..n], n in [0..10]]; // G. C. Greubel, Oct 18 2022
    
  • Mathematica
    Table[Sum[(-1)^(k-j)*Binomial[j+3, j] Binomial[n+4, k-j] (j+2)^n, {j,0,k}], {n,0, 9}, {k,0,n}]//Flatten (* Michael De Vlieger, Dec 27 2019 *)
  • PARI
    t(n,m) = if ((n<0) || (m<0), 0, if ((n==0) && (m==0), 1, (m+2)*t(n-1, m) + (n+2)*t(n, m-1)));
    tabl(nn) = {for (n=0, nn, for (k=0, n, print1(t(n-k, k), ", ");); print(););} \\ Michel Marcus, Apr 14 2015
    
  • SageMath
    def A256890(n,k): return sum((-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n for j in range(k+1))
    flatten([[A256890(n,k) for k in range(n+1)] for n in range(11)]) # G. C. Greubel, Oct 18 2022

Formula

T(n,k) = t(n-k, k); 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), where f(x) = x + 2.
Sum_{k=0..n} T(n, k) = A001715(n).
T(n,k) = Sum_{j = 0..k} (-1)^(k-j)*binomial(j+3,j)*binomial(n+4,k-j)*(j+2)^n. - Peter Bala, Dec 27 2019
Modified rule of Pascal: T(0,0) = 1, T(n,k) = 0 if k < 0 or k > n else T(n,k) = f(n-k) * T(n-1,k-1) + f(k) * T(n-1,k), where f(x) = x + 2. - Georg Fischer, Nov 11 2021
From G. C. Greubel, Oct 18 2022: (Start)
T(n, n-k) = T(n, k).
T(n, 0) = A000079(n). (End)

A085157 Quintuple factorials, 5-factorials, n!!!!!, n!5.

Original entry on oeis.org

1, 1, 2, 3, 4, 5, 6, 14, 24, 36, 50, 66, 168, 312, 504, 750, 1056, 2856, 5616, 9576, 15000, 22176, 62832, 129168, 229824, 375000, 576576, 1696464, 3616704, 6664896, 11250000, 17873856, 54286848, 119351232, 226606464, 393750000, 643458816
Offset: 0

Views

Author

Hugo Pfoertner, Jun 21 2003

Keywords

Comments

The term "Quintuple factorial numbers" is also used for the sequences A008546, A008548, A052562, A047055, A047056 which have a different definition. The definition given here is the one commonly used.

Examples

			a(12) = 168 because 12*a(12-5) = 12*a(7) = 12*14 = 168.
		

Crossrefs

Cf. n!:A000142, n!!:A006882, n!!!:A007661, n!!!!:A007662, n!!!!!!:A085158, 5-factorial primes: n!!!!!+1:A085148, n!!!!!-1:A085149.
Cf. A288092.

Programs

  • GAP
    a:= function(n)
        if n<1 then return 1;
        else return n*a(n-5);
        fi;
      end;
    List([0..40], n-> a(n) ); # G. C. Greubel, Aug 18 2019
    
  • Magma
    b:= func< n | (n lt 6) select n else n*Self(n-5) >;
    [1] cat [b(n): n in [1..40]]; // G. C. Greubel, Aug 18 2019
    
  • Maple
    a:= n-> `if`(n < 1, 1, n*a(n-5)) end proc; seq(a(n), n = 0..40); # G. C. Greubel, Aug 18 2019
  • Mathematica
    a[n_]:= If[n<1, 1, n*a[n-5]]; Table[a[n], {n,0,40}] (* G. C. Greubel, Aug 18 2019 *)
    Table[Times@@Range[n,1,-5],{n,0,40}] (* Harvey P. Dale, May 12 2020 *)
  • PARI
    a(n)=if(n<1, 1, n*a(n-5))
    for(n=0,50,print1(a(n),",")) \\ Herman Jamke (hermanjamke(AT)fastmail.fm), Oct 19 2006
    
  • Python
    def A085157(n):
        if n <= 0:
            return 1
        else:
            return n*A085157(n-5)
    n = 0
    while n <= 40:
        print(n,A085157(n))
        n = n+1 # A.H.M. Smeets, Aug 18 2019
  • Sage
    def a(n):
        if (n<1): return 1
        else: return n*a(n-5)
    [a(n) for n in (0..40)] # G. C. Greubel, Aug 18 2019
    

Formula

a(n) = 1 for n < 1, otherwise a(n) = n*a(n-5).
Sum_{n>=0} 1/a(n) = A288092. - Amiram Eldar, Nov 10 2020

A045756 Expansion of e.g.f. (1-9*x)^(-1/9), 9-factorial numbers.

Original entry on oeis.org

1, 1, 10, 190, 5320, 196840, 9054640, 498005200, 31872332800, 2326680294400, 190787784140800, 17361688356812800, 1736168835681280000, 189242403089259520000, 22330603564532623360000, 2835986652695643166720000, 385694184766607470673920000, 55925656791158083247718400000
Offset: 0

Views

Author

Keywords

Comments

Nine-fold factorials of numbers 9k+1, k = 0, 1, 2, ... - M. F. Hasler, Feb 14 2020

Crossrefs

Cf. A008542, A048994, A114806 (9-fold factorials), A132393.
Cf. k-fold factorials : A000142 ("1-fold"), A001147 (2-fold), A007559 (3), A007696 (4), A008548 (5), A008542 (6), A045754 (7), A045755 (8), A144773 (10), A256268 (combined table).

Programs

  • GAP
    List([0..20], n-> Product([0..n-1], j-> 9*j+1) ); # G. C. Greubel, Nov 11 2019
  • Magma
    [1] cat [(&*[9*j+1: j in [0..n-1]]): n in [1..20]]; // G. C. Greubel, Nov 11 2019
    
  • Maple
    seq( mul(9*j+1, j=0..n-1), n=0..20); # G. C. Greubel, Nov 11 2019
  • Mathematica
    Table[9^n*Pochhammer[1/9, n], {n,0,20}] (* G. C. Greubel, Nov 11 2019 *)
  • PARI
    vector(21, n, prod(j=0,n-2, 9*j+1) ) \\ G. C. Greubel, Nov 11 2019
    
  • Sage
    [product( (9*j+1) for j in (0..n-1)) for n in (0..20)] # G. C. Greubel, Nov 11 2019
    

Formula

a(n+1) = (9*n+1)(!^9) = Product_{k=0..n-1} (9*k+1), n >= 0.
E.g.f. (1-9*x)^(-1/9).
D-finite with recurrence: a(n) +(-9*n+8)*a(n-1)=0. - R. J. Mathar, Jan 17 2020
a(n) = A114806(9n-8). - M. F. Hasler, Feb 14 2020
a(n) = Sum_{k = 0..n} (-9)^(n - k) * A048994(n, k) = Sum_{k = 0..n} 9^(n - k) * A132393(n, k). Philippe Deléham, Sep 20 2008
a(n) = (-8)^n * sum_{k = 0..n} (9/8)^k * s(n + 1, n + 1 - k), where s(n, k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
a(n) = 9^n * Gamma(n + 1/9) / Gamma(1/9). - Artur Jasinski Aug 23 2016
a(n) ~ sqrt(2 * Pi) * 9^n * n^(n - 7/18)/(Gamma(1/9) * exp(n)). - Ilya Gutkovskiy, Sep 10 2016
Sum_{n>=0} 1/a(n) = 1 + (e/9^8)^(1/9)*(Gamma(1/9) - Gamma(1/9, 1/9)). - Amiram Eldar, Dec 21 2022

Extensions

a(0)=1 inserted; merged with A144772; formulas and programs changed accordingly by Georg Fischer, Feb 15 2020

A034301 a(n) = n-th quintic factorial number divided by 4.

Original entry on oeis.org

1, 9, 126, 2394, 57456, 1666224, 56651616, 2209413024, 97214173056, 4763494479744, 257228701906176, 15176493412464384, 971295578397720576, 67019394909442719744, 4959435223298761261056, 391795382640602139623424, 32910812141810579728367616, 2929062280621141595824717824
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • GAP
    List([1..20], n-> Product([1..n], k-> 5*k-1)/4 ); # G. C. Greubel, Aug 23 2019
  • Magma
    [&*[5*k-1: k in [1..n]]/4: n in [1..20]]; // G. C. Greubel, Aug 23 2019
    
  • Maple
    a:= n-> mul(5*k-1, k=1..n)/4: seq(a(n), n=1..20); # G. C. Greubel, Aug 23 2019
  • Mathematica
    Table[-5^(n+1)*Pochhammer[-1/5, n+1]/4, {n,20}] (* G. C. Greubel, Aug 23 2019 *)
  • PARI
    a(n) = prod(k=1,n, 5*k-1)/4;
    vector(20, n, a(n)) \\ G. C. Greubel, Aug 23 2019
    
  • Sage
    [-5^(n+1)*rising_factorial(-1/5, n+1)/4 for n in (1..20)] # G. C. Greubel, Aug 23 2019
    

Formula

a(n) = A008546(n)/4.
4*a(n) = (5*n-1)(!^5) = Product_{j=1..n} (5*j-1).
a(n) = (5*n)!/(5^n*n!*A008548(n)*2*A034323(n)*3*A034300(n)).
E.g.f.: (-1 + (1-5*x)^(-4/5))/4, a(0) = 0.
a(n) ~ sqrt(2*Pi) * 5/(4*Gamma(4/5)) * n^(13/10) * (5*n/e)^n * (1 + (241/300)/n + ...). - Joe Keane (jgk(AT)jgk.org), Nov 24 2001
D-finite with recurrence: a(n) +(-5*n+1)*a(n-1)=0. - R. J. Mathar, Feb 20 2020
Sum_{n>=1} 1/a(n) = 4*(e/5)^(1/5)*(Gamma(4/5) - Gamma(4/5, 1/5)). - Amiram Eldar, Dec 19 2022

Extensions

Terms a(17) onward added by G. C. Greubel, Aug 23 2019

A047056 Quintuple factorial numbers: Product_{k=0..n-1} (5*k+3).

Original entry on oeis.org

1, 3, 24, 312, 5616, 129168, 3616704, 119351232, 4535346816, 195019913088, 9360955828224, 496130658895872, 28775578215960576, 1812861427605516288, 123274577077175107584, 8999044126633782853632, 701925441877435062583296, 58259811675827110194413568
Offset: 0

Views

Author

Joe Keane (jgk(AT)jgk.org)

Keywords

Crossrefs

Programs

  • GAP
    List([0..20], n-> Product([0..n-1], k-> 5*k+3 )); # G. C. Greubel, Aug 20 2019
  • Magma
    [1] cat [(&*[5*k+3: k in [0..n-1]]): n in [1..20]]; // G. C. Greubel, Aug 20 2019
    
  • Maple
    seq(mul(5*k+3, k = 0..n-1), n = 0..20); # G. C. Greubel, Aug 20 2019
  • Mathematica
    Table[5^n*Pochhammer[3/5, n], {n,0,20}] (* G. C. Greubel, Aug 20 2019 *)
    Join[{1},FoldList[Times,5*Range[0,20]+3]] (* Harvey P. Dale, Oct 08 2020 *)
  • PARI
    vector(20, n, n--; prod(j=0,n-1, 5*j+3) ) \\ G. C. Greubel, Aug 20 2019
    
  • Sage
    [5^n*rising_factorial(3/5, n) for n in (0..20)] # G. C. Greubel, Aug 20 2019
    

Formula

E.g.f.: (1-5*x)^(-3/5).
a(n) ~ sqrt(2*Pi)/Gamma(3/5)*n^(1/10)*(5*n/e)^n*(1 - 11/300*n^ - 1 + ...). - Joe Keane (jgk(AT)jgk.org), Nov 24 2001
G.f.: 1/(1-3*x/(1-5*x/(1-8*x/(1-10*x/(1-13*x/(1-15*x/(1-18*x/(1-20*x/(1- ... (continued fraction). - Philippe Deléham, Jan 08 2012
a(n) = (-2)^n*Sum_{k=0..n} (5/2)^k*s(n+1,n+1-k), where s(n,k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
G.f.: 1/G(0) where G(k) = 1 - x*(5*k+3)/( 1 - 5*x*(k+1)/G(k+1) ); (continued fraction ). - Sergei N. Gladkovskii, Mar 23 2013
G.f.: G(0)/2, where G(k)= 1 + 1/(1 - (5*k+3)*x/((5*k+3)*x + 1/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 14 2013
D-finite with recurrence: a(n) +(-5*n+2)*a(n-1)=0. - R. J. Mathar, Jan 17 2020
Sum_{n>=0} 1/a(n) = 1 + (e/5^2)^(1/5)*(Gamma(3/5) - Gamma(3/5, 1/5)). - Amiram Eldar, Dec 19 2022
Previous Showing 11-20 of 60 results. Next