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.

A316143 Expansion of e.g.f. Product_{k>=1} 1 / (1 - (exp(x)-1)^k)^2.

Original entry on oeis.org

1, 2, 12, 92, 912, 10772, 148512, 2328692, 40842912, 791302772, 16767551712, 385382491892, 9542377300512, 253105962752372, 7156766466076512, 214814484529608692, 6819311473596695712, 228212485803422931572, 8028037725386962194912, 296094910181041530831092
Offset: 0

Views

Author

Vaclav Kotesovec, Jun 25 2018

Keywords

Comments

Self-convolution of A167137.
Conjecture: Let k be a positive integer. The sequence obtained by reducing a(n) modulo k is eventually periodic with the period dividing phi(k) = A000010(k). For example, modulo 7 we obtain the sequence [1, 2, 5, 1, 2, 6, 0, 2, 5, 1, 2, 6, 0, 2, 5, 1, 2, 6, 0, ...], with a preperiod of length 1 and an apparent period thereafter of 6 = phi(7). - Peter Bala, Mar 03 2023

Crossrefs

Programs

  • Mathematica
    nmax = 20; CoefficientList[Series[Product[1/(1-(Exp[x]-1)^k)^2, {k, 1, nmax}], {x, 0, nmax}], x] * Range[0, nmax]!

Formula

Sum_{k=0..n} binomial(n,k) * A167137(k) * A167137(n-k).
a(n) ~ n! * exp(Pi * sqrt(2*n/(3*log(2))) - Pi^2 * (1 - 1/log(2)) / 12) / (2^(7/4) * 3^(3/4) * n^(5/4) * (log(2))^(n - 1/4)).