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-3 of 3 results.

A317855 Decimal expansion of a constant related to the asymptotics of A122400.

Original entry on oeis.org

3, 1, 6, 1, 0, 8, 8, 6, 5, 3, 8, 6, 5, 4, 2, 8, 8, 1, 3, 8, 3, 0, 1, 7, 2, 2, 0, 2, 5, 8, 8, 1, 3, 2, 4, 9, 1, 7, 2, 6, 3, 8, 2, 7, 7, 4, 1, 8, 8, 5, 5, 6, 3, 4, 1, 6, 2, 7, 2, 7, 8, 2, 0, 7, 5, 3, 7, 6, 9, 7, 0, 5, 9, 2, 1, 9, 3, 0, 4, 6, 1, 1, 2, 1, 9, 7, 5, 7, 4, 6, 8, 5, 4, 9, 7, 8, 4, 5, 9, 3, 2, 4, 2, 2, 7
Offset: 1

Views

Author

Vaclav Kotesovec, Aug 09 2018

Keywords

Examples

			3.161088653865428813830172202588132491726382774188556341627278...
		

Crossrefs

Programs

  • Mathematica
    r = r /. FindRoot[E^(1/r)/r + (1 + E^(1/r)) * ProductLog[-E^(-1/r)/r] == 0, {r, 3/4}, WorkingPrecision -> 120]; RealDigits[(1 + Exp[1/r])*r^2][[1]]
  • PARI
    r=solve(r=.8,1,exp(1/r)/r + (1+exp(1/r))*lambertw(-exp(-1/r)/r))
    (1+exp(1/r))*r^2 \\ Charles R Greathouse IV, Jun 15 2021

Formula

Equals (1+exp(1/r))*r^2, where r = 0.873702433239668330496568304720719298213992... is the root of the equation exp(1/r)/r + (1+exp(1/r))*LambertW(-exp(-1/r)/r) = 0.

A244585 E.g.f.: Sum_{n>=1} (exp(n*x) - 1)^n / n.

Original entry on oeis.org

1, 5, 79, 2621, 149071, 12954365, 1596620719, 264914218301, 56934521042191, 15385666763366525, 5106110041462786159, 2041611328770984737981, 967972254733121945653711, 536962084044317668770841085, 344546100916295014902350596399
Offset: 1

Views

Author

Paul D. Hanna, Aug 21 2014

Keywords

Comments

Compare to: Sum_{n>=1} (1 - exp(-n*x))^n / n, the e.g.f. of A092552.

Examples

			E.g.f.: A(x) = x + 5*x^2/2! + 79*x^3/3! + 2621*x^4/4! + 149071*x^5/5! +...
where
A(x) = (exp(x)-1) + (exp(2*x)-1)^2/2 + (exp(3*x)-1)^3/3 + (exp(4*x)-1)^4/4 + (exp(5*x)-1)^5/5 + (exp(6*x)-1)^6/6 + (exp(7*x)-1)^7/7 +...
Exponentiation yields:
exp(A(x)) = 1 + x + 6*x^2/2! + 95*x^3/3! + 3043*x^4/4! + 167342*x^5/5! +...+ A243802(n)*x^n/n! +...
The O.G.F. begins:
F(x) = x + 5*x^2 + 79*x^3 + 2621*x^4 + 149071*x^5 + 12954365*x^6 +...
where
F(x) = x/(1-x) + 2*2!*x^2/((1-2*x)*(1-4*x)) + 3^2*3!*x^3/((1-3*x)*(1-6*x)*(1-9*x)) + 4^3*4!*x^4/((1-4*x)*(1-8*x)*(1-12*x)*(1-16*x)) + 5^4*5!*x^5/((1-5*x)*(1-10*x)*(1-15*x)*(1-20*x)*(1-25*x)) +...
		

Crossrefs

Programs

  • PARI
    {a(n) = n!*polcoeff( sum(m=1,n+1, (exp(m*x +x*O(x^n)) - 1)^m / m), n)}
    for(n=0,20,print1(a(n),", "))
    
  • PARI
    {a(n)=if(n<1, 0, polcoeff(sum(m=1, n, m^(m-1) * m! * x^m / prod(k=1, m, 1-m*k*x +x*O(x^n))), n))}
    for(n=0, 20, print1(a(n), ", "))

Formula

O.g.f.: Sum_{n>=1} n^(n-1) * n! * x^n / Product_{k=1..n} (1 - n*k*x).
a(n) ~ c * d^n * (n!)^2 / n^(3/2), where d = A317855 = (1+exp(1/r))*r^2 = 3.161088653865428813830172202588132491..., r = 0.873702433239668330496568304720719298... is the root of the equation exp(1/r)/r + (1+exp(1/r)) * LambertW(-exp(-1/r)/r) = 0, and c = 0.37498840921734807101035131780130551... . - Vaclav Kotesovec, Aug 21 2014

A244437 E.g.f.: exp( Sum_{n>=1} (1 - exp(-n*x))^n / n ).

Original entry on oeis.org

1, 1, 4, 41, 845, 30012, 1650475, 130216865, 13944696526, 1945060435587, 342412144747677, 74216506678085290, 19414505134246518741, 6029823819095965829293, 2193174302711080501699684, 923346371767630311443639677, 445468655004100653462280596881, 244137607569262412209821327718964
Offset: 0

Views

Author

Paul D. Hanna, Aug 21 2014

Keywords

Comments

Compare to: exp( Sum_{n>=1} (1 - exp(-x))^n/n ) = 1/(2-exp(x)), the e.g.f. of Fubini numbers (A000670).

Examples

			E.g.f.: A(x) = 1 + x + 4*x^2/2! + 41*x^3/3! + 845*x^4/4! + 30012*x^5/5! +...
where
log(A(x)) = (1-exp(-x)) + (1-exp(-2*x))^2/2 + (1-exp(-3*x))^3/3 + (1-exp(-4*x))^4/4 + (1-exp(-5*x))^5/5 + (1-exp(-6*x))^6/6 +...
Explicitly,
log(A(x)) = x + 3*x^2/2! + 31*x^3/3! + 675*x^4/4! + 25231*x^5/5! + 1441923*x^6/6! +...+ A092552(n)*x^n/n! +...
		

Crossrefs

Programs

  • Mathematica
    max = 20; s = Exp[Sum[(1 - Exp[-n x])^n/n, {n, 1, max}]] + O[x]^max; CoefficientList[s, x] Range[0, max-1]! (* Jean-François Alcover, Mar 31 2016 *)
  • PARI
    {a(n) = n!*polcoeff( exp( sum(m=1,n+1, (1 - exp(-m*x +x*O(x^n)))^m / m) ), n)}
    for(n=0,20,print1(a(n),", "))

Formula

E.g.f.: exp( Sum_{n>=1} A092552(n)*x^n/n! ), where A092552(n) = Sum_{k=1..n} k!*(k-1)! * Stirling2(n, k)^2.
a(n) ~ (n!)^2 / (2 * sqrt(Pi) * sqrt(1-log(2)) * n^(3/2) * log(2)^(2*n)). - Vaclav Kotesovec, Aug 21 2014
Showing 1-3 of 3 results.