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

A007405 Dowling numbers: e.g.f.: exp(x + (exp(b*x) - 1)/b) with b=2.

Original entry on oeis.org

1, 2, 6, 24, 116, 648, 4088, 28640, 219920, 1832224, 16430176, 157554048, 1606879040, 17350255744, 197553645440, 2363935624704, 29638547505408, 388328781668864, 5304452565517824, 75381218537805824, 1112348880749130752, 17014743624340539392, 269360902955086379008
Offset: 0

Views

Author

Keywords

Comments

Binomial transform of A004211.
Equals leftmost term in iterates of M^n * [1,1,1,...], where M = a bidiagonal matrix with (1,3,5,7,...) in the main diagonal and (1,1,1,...) in the superdiagonal. - Gary W. Adamson, Apr 13 2009
This is the number of type B set partitions, see R. Suter. - Per W. Alexandersson, Dec 19 2022

Examples

			a(4) = 116 = sum of top row terms of M^3 = (49 + 44 + 18 + 4 + 1).
		

References

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

Crossrefs

Cf. A000110 (b=1), this sequence (b=2), A003575 (b=3), A003576 (b=4), A003577 (b=5), A003578 (b=6), A003579 (b=7), A003580 (b=8), A003581 (b=9), A003582 (b=10).

Programs

  • Magma
    m:=20; c:=2; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x +(Exp(c*x)-1)/c) )); [Factorial(n-1)*b[n]: n in [1..m-1]]; // G. C. Greubel, Feb 24 2019
  • Mathematica
    max = 19; f[x_]:= Exp[x + Exp[2x]/2 -1/2]; CoefficientList[Series[f[x], {x,0,max}], x]*Range[0, max]! (* Jean-François Alcover, Nov 22 2011 *)
    Table[Sum[Binomial[n, k] * 2^k * BellB[k, 1/2], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Apr 17 2020 *)
  • PARI
    x='x+O('x^66); Vec(serlaplace(exp(x+1/2*exp(2*x)-1/2))) \\ Joerg Arndt, May 13 2013
    
  • Sage
    @CachedFunction
    def S(n, k, m):
        if k > n or k < 0 : return 0
        if n == 0 and k == 0: return 1
        return S(n-1, k-1, m) + (m*(k+1)-1)*S(n-1, k, m)
    def A007405(n): return add(S(n, k, 2) for k in (0..n)) # Peter Luschny, May 20 2013
    
  • Sage
    b=2;
    def A007405_list(prec):
        P. = PowerSeriesRing(QQ, prec)
        return P( exp(x +(exp(b*x)-1)/b) ).egf_to_ogf().list()
    A007405_list(30) # G. C. Greubel, Feb 24 2019
    

Formula

E.g.f.: exp(x + (exp(2*x) - 1)/2).
Row sums of triangles A039755, A039756. - Philippe Deléham, Feb 20 2005
a(n) = sum of top row terms of M^n, M = an infinite square production matrix in which a diagonal of 1's is appended to the right of Pascal's triangle squared; as follows:
1, 1, 0, 0, 0, 0, ...
2, 1, 1, 0, 0, 0, ...
4, 4, 1, 1, 0, 0, ...
8, 12, 6, 1, 1, 0, ...
16, 32, 24, 8, 1, 1, ...
... - Gary W. Adamson, Aug 01 2011
G.f.: (G(0) - 1)/(x-1) where G(k) = 1 - 1/(1-(2*k+1)*x)/(1-x/(x-1/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 16 2013
G.f.: -G(0) where G(k) = 1 - (x*(2*k+1) - 2)/(x*(2*k+1) - 1 - x*(x*(2*k+1) - 1)/(x + (x*(2*k+1) - 2)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 29 2013
G.f.: 1/Q(0), where Q(k) = 1 - 2*(k+1)*x - 2*(k+1)*x^2/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, May 03 2013
G.f.: 1/Q(0), where Q(k) = 1 - x - x/(1 - x*(2*k+2)/Q(k+1)); (continued fraction). - Sergei N. Gladkovskii, May 13 2013
G.f.: 1/(1-x*Q(0)), where Q(k) = 1 + x/(1 - x + 2*x*(k+1)/(x - 1/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 19 2013
Conjecture: Let M_n be an n X n matrix whose elements are m_ij = 1 for i < j - 1, m_ij = -1 for i = j - 1, and m_ij = binomial(n - i,j - i) otherwise. Then a(n - 1) = Det(M_n). - Benedict W. J. Irwin, Apr 19 2017
a(n) = exp(-1/2) * Sum_{k>=0} (2*k + 1)^n / (2^k * k!). - Ilya Gutkovskiy, Apr 16 2020
a(n) = Sum_{k=0..n} binomial(n,k) * A187251(k) * A187251(n-k). - Vaclav Kotesovec, Apr 17 2020
a(n) ~ 2^(n + 1/2) * n^(n + 1/2) * exp(n/LambertW(2*n) - n - 1/2) / (sqrt(1 + LambertW(2*n)) * LambertW(2*n)^(n + 1/2)). - Vaclav Kotesovec, Jun 26 2022

Extensions

Name edited by G. C. Greubel, Feb 24 2019

A003575 Dowling numbers: e.g.f.: exp(x + (exp(b*x) - 1)/b) with b=3.

Original entry on oeis.org

1, 2, 7, 35, 214, 1523, 12349, 112052, 1120849, 12219767, 143942992, 1819256321, 24526654381, 350974470746, 5308470041299, 84554039118383, 1413794176669942, 24745966692370607, 452277149756692105, 8612255652371171012, 170517319084490074405
Offset: 0

Views

Author

Keywords

Comments

Named after the American mathematician Thomas Allan Dowling (b. 1941). - Amiram Eldar, Jun 06 2021

Crossrefs

Cf. A000110 (b=1), A007405 (b=2), this sequence (b=3), A003576 (b=4), A003577 (b=5), A003578 (b=6), A003579 (b=7), A003580 (b=8), A003581 (b=9), A003582 (b=10).

Programs

  • Magma
    m:=30; c:=3; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x+(Exp(c*x)-1)/c) )); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Feb 20 2019
    
  • Maple
    seq(coeff(series(n!*exp(z+(1/3)*exp(3*z)-(1/3)),z,n+1), z, n), n=0..30); # Muniru A Asiru, Feb 19 2019
  • Mathematica
    With[{nn=20},CoefficientList[Series[Exp[x+Exp[3x]/3-1/3],{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Jan 04 2019 *)
    Table[Sum[Binomial[n, k] * 3^k * BellB[k, 1/3], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Apr 17 2020 *)
  • PARI
    my(x = 'x + O('x^30)); Vec(serlaplace(exp(x + exp(3*x)/3 - 1/3))) \\ Michel Marcus, Feb 09 2018
    
  • Sage
    b=3;
    def A003575_list(prec):
        P. = PowerSeriesRing(QQ, prec)
        return P( exp(x +(exp(b*x)-1)/b) ).egf_to_ogf().list()
    A003575_list(30) # G. C. Greubel, Feb 20 2019

Formula

E.g.f.: exp(x + (exp(3*x) - 1)/3).
G.f.: 1/(1-x*Q(0)), where Q(k) = 1 + x/(1 - x + 3*x*(k+1)/(x - 1/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 19 2013
a(n) = exp(-1/3) * Sum_{k>=0} (3*k + 1)^n / (3^k * k!). - Ilya Gutkovskiy, Apr 16 2020
a(n) ~ 3^(n + 1/3) * n^(n + 1/3) * exp(n/LambertW(3*n) - n - 1/3) / (sqrt(1 + LambertW(3*n)) * LambertW(3*n)^(n + 1/3)). - Vaclav Kotesovec, Jun 26 2022

Extensions

Name clarified by G. C. Greubel, Feb 20 2019

A003576 Dowling numbers: e.g.f.: exp(x + (exp(b*x) - 1)/b) with b=4.

Original entry on oeis.org

1, 2, 8, 48, 352, 3008, 29440, 324096, 3947520, 52541440, 757260288, 11733385216, 194272854016, 3419584921600, 63707979972608, 1251489089060864, 25836869372608512, 558946705406427136, 12638569755079344128, 298003073694026432512, 7312035980392431353856
Offset: 0

Views

Author

Keywords

Crossrefs

Cf. A000110 (b=1), A007405 (b=2), A003575 (b=3), this sequence (b=4), A003577 (b=5), A003578 (b=6), A003579 (b=7), A003580 (b=8), A003581 (b=9), A003582 (b=10).

Programs

  • Magma
    m:=20; c:=4; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x+(Exp(c*x)-1)/c) )); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Feb 22 2019
    
  • Maple
    seq(coeff(series(factorial(n)*exp(z+(1/4)*exp(4*z)-(1/4)),z,n+1), z, n), n = 0 .. 20); # Muniru A Asiru, Feb 22 2019
  • Mathematica
    With[{m=20, b=4}, CoefficientList[Series[Exp[x+(Exp[b*x]-1)/b], {x,0,m}], x]*Range[0, m]!] (* G. C. Greubel, Feb 22 2019 *)
    Table[Sum[Binomial[n, k] * 4^k * BellB[k, 1/4], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Apr 17 2020 *)
  • PARI
    my(x='x+O('x^20)); b=4; Vec(serlaplace(exp(x+(exp(b*x)-1)/b))) \\ G. C. Greubel, Feb 22 2019
    
  • Sage
    m = 20; b=4; T = taylor(exp(x+(exp(b*x)-1)/b), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, Feb 22 2019

Formula

E.g.f.: exp(z + (exp(4*z) - 1)/4).
G.f.: 1/Q(0), where Q(k) = 1 - 2*x*(2*k+1) - 2*x^2*(2*k+2)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Sep 26 2013
a(n) = exp(-1/4) * Sum_{k>=0} (4*k + 1)^n / (4^k * k!). - Ilya Gutkovskiy, Apr 16 2020
a(n) ~ 4^(n + 1/4) * n^(n + 1/4) * exp(n/LambertW(4*n) - n - 1/4) / (sqrt(1 + LambertW(4*n)) * LambertW(4*n)^(n + 1/4)). - Vaclav Kotesovec, Jun 26 2022

A003577 Dowling numbers: e.g.f. exp(x + (exp(b*x)-1)/b) with b=5.

Original entry on oeis.org

1, 2, 9, 63, 536, 5307, 60389, 775988, 11062391, 172638727, 2921519374, 53221709973, 1037320865141, 21517178350762, 472862758184789, 10966587174511443, 267502464814857936, 6842498829509972687, 183057455239626138009, 5110016898453125496548
Offset: 0

Views

Author

Keywords

Crossrefs

Cf. A000110 (b=1), A007405 (b=2), A003575 (b=3), A003576 (b=4), this sequence (b=5), A003578 (b=6), A003579 (b=7), A003580 (b=8), A003581 (b=9), A003582 (b=10).

Programs

  • GAP
    b:=5;; a:=[1,2];; for n in [3..20] do a[n]:=2*a[n-1]+Sum([0..n-3],i->Binomial(n-2,i)*b^(n-2-i)*a[i+1]); od; Print(a); # Muniru A Asiru, Apr 10 2019
  • Magma
    m:=20; c:=5; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x +(Exp(c*x)-1)/c) )); [Factorial(n-1)*b[n]: n in [1..m-1]]; // G. C. Greubel, Feb 24 2019
    
  • Mathematica
    With[{m=20, b=5}, CoefficientList[Series[Exp[x +(Exp[b*x]-1)/b],{x,0,m}], x]*Range[0, m]!] (* G. C. Greubel, Feb 24 2019 *)
    Table[Sum[Binomial[n, k] * 5^k * BellB[k, 1/5], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Apr 17 2020 *)
  • PARI
    my(x='x+O('x^20)); b=5; Vec(serlaplace(exp(x +(exp(b*x)-1)/b))) \\ G. C. Greubel, Feb 24 2019
    
  • Sage
    m = 20; b=5; T = taylor(exp(x + (exp(b*x) -1)/b), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, Feb 24 2019
    

Formula

E.g.f.: exp(x + (exp(5*x) - 1)/5).
a(n) = exp(-1/5) * Sum_{k>=0} (5*k + 1)^n / (5^k * k!). - Ilya Gutkovskiy, Apr 16 2020
a(n) ~ 5^(n + 1/5) * n^(n + 1/5) * exp(n/LambertW(5*n) - n - 1/5) / (sqrt(1 + LambertW(5*n)) * LambertW(5*n)^(n + 1/5)). - Vaclav Kotesovec, Jun 26 2022

Extensions

Name clarified by Muniru A Asiru, Feb 24 2019

A003578 Dowling numbers: e.g.f. exp(x + (exp(b*x) - 1)/b), with b=6.

Original entry on oeis.org

1, 2, 10, 80, 772, 8648, 111592, 1631360, 26518672, 472528160, 9139219360, 190461416192, 4250569655872, 101040920561792, 2546488866632320, 67772341398044672, 1898177372174512384, 55780954727160472064, 1715291443214323558912, 55062161002484359565312
Offset: 0

Views

Author

Keywords

Crossrefs

Cf. A000110 (b=1), A007405 (b=2), A003575 (b=3), A003576 (b=4), A003577 (b=5), this sequence (b=6), A003579 (b=7), A003580 (b=8), A003581 (b=9), A003582 (b=10).

Programs

  • Magma
    m:=20; c:=6; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x +(Exp(c*x)-1)/c) )); [Factorial(n-1)*b[n]: n in [1..m-1]]; // G. C. Greubel, Feb 24 2019
    
  • Maple
    seq(coeff(series(factorial(n)*exp(z+(1/6)*exp(6*z)-(1/6)),z,n+1), z, n), n = 0 .. 20); # Muniru A Asiru, Feb 23 2019
  • Mathematica
    With[{nn=20},CoefficientList[Series[Exp[x+Exp[6x]/6-1/6],{x,0,nn}],x] Range[0,nn]!] (* Harvey P. Dale, Jun 12 2017 *)
    Table[Sum[Binomial[n, k] * 6^k * BellB[k, 1/6], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Apr 17 2020 *)
  • PARI
    my(x='x+O('x^20)); b=6; Vec(serlaplace(exp(x +(exp(b*x)-1)/b))) \\ G. C. Greubel, Feb 24 2019
    
  • Sage
    m = 20; b=6; T = taylor(exp(x + (exp(b*x) -1)/b), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, Feb 24 2019

Formula

E.g.f.: exp(x + (exp(6*x) - 1)/6).
a(n) = exp(-1/6) * Sum_{k>=0} (6*k + 1)^n / (6^k * k!). - Ilya Gutkovskiy, Apr 16 2020
a(n) ~ 6^(n + 1/6) * n^(n + 1/6) * exp(n/LambertW(6*n) - n - 1/6) / (sqrt(1 + LambertW(6*n)) * LambertW(6*n)^(n + 1/6)). - Vaclav Kotesovec, Jun 26 2022

A003580 Dowling numbers: e.g.f. exp(x + (exp(b*x)-1)/b) with b=8.

Original entry on oeis.org

1, 2, 12, 120, 1424, 19488, 307904, 5539712, 111259904, 2454487552, 58847153152, 1522019629056, 42209521995776, 1248370355347456, 39186678731423744, 1300179383923212288, 45436201241711542272, 1667242078056889843712, 64063345344029286727680
Offset: 0

Views

Author

Keywords

Crossrefs

Cf. A000110 (b=1), A007405 (b=2), A003575 (b=3), A003576 (b=4), A003577 (b=5), A003578 (b=6), A003579 (b=7), this sequence (b=8), A003581 (b=9), A003582 (b=10).

Programs

  • Magma
    m:=20; c:=8; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x +(Exp(c*x)-1)/c) )); [Factorial(n-1)*b[n]: n in [1..m-1]]; // G. C. Greubel, Feb 24 2019
    
  • Maple
    seq(coeff(series(factorial(n)*exp(z+(1/8)*exp(8*z)-(1/8)),z,n+1), z, n), n = 0 .. 20); # Muniru A Asiru, Feb 24 2019
  • Mathematica
    With[{m=20, b=8}, CoefficientList[Series[Exp[x +(Exp[b*x]-1)/b],{x,0,m}], x]*Range[0, m]!] (* G. C. Greubel, Feb 24 2019 *)
    Table[Sum[Binomial[n, k] * 8^k * BellB[k, 1/8], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Apr 17 2020 *)
  • PARI
    my(x='x+O('x^20)); b=8; Vec(serlaplace(exp(x +(exp(b*x)-1)/b))) \\ G. C. Greubel, Feb 24 2019
    
  • Sage
    m = 20; b=8; T = taylor(exp(x + (exp(b*x) -1)/b), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, Feb 24 2019

Formula

E.g.f.: exp(x + (exp(8*x) - 1)/8).
a(n) = exp(-1/8) * Sum_{k>=0} (8*k + 1)^n / (8^k * k!). - Ilya Gutkovskiy, Apr 16 2020
a(n) ~ 8^(n + 1/8) * n^(n + 1/8) * exp(n/LambertW(8*n) - n - 1/8) / (sqrt(1 + LambertW(8*n)) * LambertW(8*n)^(n + 1/8)). - Vaclav Kotesovec, Jun 26 2022

Extensions

Name clarified by Muniru A Asiru, Feb 24 2019

A003579 Dowling numbers: e.g.f. exp(x + (exp(b*x) - 1)/b), with b=7.

Original entry on oeis.org

1, 2, 11, 99, 1066, 13283, 190933, 3117900, 56729565, 1132679479, 24564972756, 574431351673, 14394977015245, 384489778509034, 10894501505088695, 326149933663962479, 10280153573323314858
Offset: 0

Views

Author

Keywords

Crossrefs

Cf. A000110 (b=1), A007405 (b=2), A003575 (b=3), A003576 (b=4), A003577 (b=5), A003578 (b=6), this sequence (b=7), A003580 (b=8), A003581 (b=9), A003582 (b=10).

Programs

  • GAP
    b:=7;; a:=[1,2];; for n in [3..20] do a[n]:=2*a[n-1]+Sum([0..n-3],i->Binomial(n-2,i)*b^(n-2-i)*a[i+1]); od; Print(a); # Muniru A Asiru, Apr 10 2019
  • Magma
    m:=20; c:=7; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x +(Exp(c*x)-1)/c) )); [Factorial(n-1)*b[n]: n in [1..m-1]]; // G. C. Greubel, Feb 24 2019
    
  • Mathematica
    With[{m=20, b=7}, CoefficientList[Series[Exp[x +(Exp[b*x]-1)/b],{x,0,m}], x]*Range[0, m]!] (* G. C. Greubel, Feb 24 2019 *)
    Table[Sum[Binomial[n, k] * 7^k * BellB[k, 1/7], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Apr 17 2020 *)
  • PARI
    Vec(serlaplace( exp(z + 1/7 * exp(7 * z) - 1/7) ) ) \\ Joerg Arndt, Feb 24 2019
    
  • Sage
    m = 20; b=7; T = taylor(exp(x + (exp(b*x) -1)/b), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, Feb 24 2019
    

Formula

E.g.f.: exp(x + (exp(7*x) - 1)/7).
a(n) = exp(-1/7) * Sum_{k>=0} (7*k + 1)^n / (7^k * k!). - Ilya Gutkovskiy, Apr 16 2020
a(n) ~ 7^(n + 1/7) * n^(n + 1/7) * exp(n/LambertW(7*n) - n - 1/7) / (sqrt(1 + LambertW(7*n)) * LambertW(7*n)^(n + 1/7)). - Vaclav Kotesovec, Jun 26 2022

Extensions

Name clarified by Muniru A Asiru, Feb 24 2019

A003582 Dowling numbers: e.g.f. exp(x + (exp(b*x)-1)/b) with b=10.

Original entry on oeis.org

1, 2, 14, 168, 2356, 37832, 701464, 14866848, 352943376, 9219925792, 261954304224, 8033968939648, 264411579439936, 9288709762556032, 346608927301622144, 13680000261825018368, 569006722158124974336, 24864267879086770135552, 1138321277772163220033024
Offset: 0

Views

Author

Keywords

Comments

In general, for b > 0, if e.g.f. = exp(x + (exp(b*x) - 1)/b), then a(n) ~ b^(n + 1/b) * n^(n + 1/b) * exp(n/LambertW(b*n) - n - 1/b) / (sqrt(1 + LambertW(b*n)) * LambertW(b*n)^(n + 1/b)). - Vaclav Kotesovec, Jun 26 2022

Crossrefs

Cf. A000110 (b=1), A007405 (b=2), A003575 (b=3), A003576 (b=4), A003577 (b=5), A003578 (b=6), A003579 (b=7), A003580 (b=8), A003581 (b=9), this sequence (b=10).

Programs

  • Magma
    m:=20; c:=10; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( Exp(x +(Exp(c*x)-1)/c) )); [Factorial(n-1)*b[n]: n in [1..m-1]]; // G. C. Greubel, Feb 24 2019
    
  • Maple
    seq(coeff(series(factorial(n)*exp(z+(1/10)*exp(10*z)-(1/10)),z,n+1), z, n), n = 0 .. 20); # Muniru A Asiru, Feb 24 2019
  • Mathematica
    With[{m=20, b=10}, CoefficientList[Series[Exp[x +(Exp[b*x]-1)/b],{x,0,m}], x]*Range[0, m]!] (* G. C. Greubel, Feb 24 2019 *)
    Table[Sum[Binomial[n, k] * 10^k * BellB[k, 1/10], {k, 0, n}], {n, 0, 20}] (* Vaclav Kotesovec, Apr 17 2020 *)
  • PARI
    my(x='x+O('x^20)); b=10; Vec(serlaplace(exp(x +(exp(b*x)-1)/b))) \\ G. C. Greubel, Feb 24 2019
    
  • Sage
    m = 20; b=10; T = taylor(exp(x + (exp(b*x) -1)/b), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, Feb 24 2019

Formula

E.g.f.: exp(x + (exp(10*x) - 1)/10).
a(n) = exp(-1/10) * Sum_{k>=0} (10*k + 1)^n / (10^k * k!). - Ilya Gutkovskiy, Apr 16 2020
a(n) ~ 10^(n + 1/10) * n^(n + 1/10) * exp(n/LambertW(10*n) - n - 1/10) / (sqrt(1 + LambertW(10*n)) * LambertW(10*n)^(n + 1/10)). - Vaclav Kotesovec, Jun 26 2022

Extensions

Name clarified by Muniru A Asiru, Feb 24 2019

A334162 a(0) = 1; thereafter a(n) = exp(-1/n) * Sum_{k>=0} (n*k + 1)^n / (n^k * k!).

Original entry on oeis.org

1, 2, 6, 35, 352, 5307, 111592, 3117900, 111259904, 4912490375, 261954304224, 16560019685937, 1222893826048000, 104189533522270666, 10132262911996769408, 1114216450970154278543, 137427598621356912082944, 18877351974681584403701519, 2869969478954093766868948480
Offset: 0

Views

Author

Ilya Gutkovskiy, Apr 16 2020

Keywords

Crossrefs

Programs

  • Mathematica
    Table[SeriesCoefficient[1/(1 - x) Sum[(x/(1 - x))^k/Product[(1 - n j x/(1 - x)), {j, 1, k}], {k, 0, n}], {x, 0, n}], {n, 0, 18}]
    Join[{1}, Table[n! SeriesCoefficient[Exp[x + (Exp[n x] - 1)/n], {x, 0, n}], {n, 1, 18}]]

Formula

a(n) = [x^n] (1/(1 - x)) * Sum_{k>=0} (x/(1 - x))^k / Product_{j=1..k} (1 - n*j*x/(1 - x)).
a(n) = n! * [x^n] exp(x + (exp(n*x) - 1) / n), for n > 0.
a(n) = A334165(n,n).

A364069 Dowling numbers: e.g.f. exp(x + (exp(b*x)-1)/b) with b = 63.

Original entry on oeis.org

1, 2, 67, 4355, 295234, 21036803, 1625419909, 140823067772, 13947448935109, 1570142163116087, 196457384808738412, 26717651072732512841, 3896182904620308595021, 605803757139146097600266, 100236348400243756326661039, 17619174544126256877550593743, 3280792242500933388439611444802
Offset: 0

Views

Author

Stefano Spezia, Jul 04 2023

Keywords

Comments

a(n) is the number of all 64-subgroups of R^n, where R^n is a near-vector space over a proper nearfield R.

Crossrefs

Cf. A000110 (b=1), A007405 (b=2), A003575 (b=3), A003576 (b=4), A003577 (b=5), A003578 (b=6), A003579 (b=7), A003580 (b=8), A003581 (b=9), A003582 (b=10), A364070 (b=624).
Row sums of the triangle A364072.
2nd row of the array A364074.

Programs

  • Mathematica
    With[{m=16, b=63}, CoefficientList[Series[Exp[x +(Exp[b*x]-1)/b], {x, 0, m}], x]*Range[0, m]!] (* or *)
    a[n_]:=Sum[Sum[Binomial[n,d]StirlingS2[n-d,k]63^(n-d-k),{d,0,n-k}],{k,0,n}]; Array[a,17,0]

Formula

E.g.f.: exp(x + (exp(63*x) - 1)/63).
a(n) = exp(-1/63) * Sum_{k>=0} (63*k + 1)^n / (63^k * k!).
a(n) ~ 63^(n + 1/63) * n^(n + 1/63) * exp(n/LambertW(63*n) - n - 1/63) / (sqrt(1 + LambertW(63*n)) * LambertW(63*n)^(n + 1/63)).
a(n) = Sum_{k=0..n} Sum_{d=0..n-k} binomial(n, d)*StirlingS2(n-d, k)*63^(n-d-k).
Showing 1-10 of 12 results. Next