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.

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