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.
%I A266941 #27 Aug 27 2018 05:46:31 %S A266941 1,1,5,14,42,103,289,690,1771,4206,10142,23449,54786,123528,279480, %T A266941 619206,1366405,2969071,6425534,13727775,29187555,61439660,128620370, %U A266941 267044222,551527679,1130806020,2306746335,4676096006,9432394144,18920266428,37776372312 %N A266941 Expansion of Product_{k>=1} 1 / (1 - k*x^k)^k. %C A266941 This sequence is obtained from the generalized Euler transform in A266964 by taking f(n) = n, g(n) = n. - _Seiichi Manyama_, Nov 18 2017 %H A266941 Seiichi Manyama, <a href="/A266941/b266941.txt">Table of n, a(n) for n = 0..6086</a> (terms 0..1000 from Vaclav Kotesovec) %F A266941 From _Vaclav Kotesovec_, Jan 08 2016: (Start) %F A266941 a(n) ~ c * n^2 * 3^(n/3), where %F A266941 c = 3278974684037157122864203.021982619109776972432419491714093... if mod(n,3)=0 %F A266941 c = 3278974684037157122864202.999526122508793149896683112820555... if mod(n,3)=1 %F A266941 c = 3278974684037157122864203.001231135511323719311281438384212... if mod(n,3)=2 %F A266941 (End) %F A266941 In closed form, a(n) ~ (Product_{k>=4}((1 - k/3^(k/3))^(-k)) / ((1 - 2/3^(2/3))^2 * (1 - 1/3^(1/3))) + Product_{k>=4}((1 - (-1)^(2*k/3)*k/3^(k/3))^(-k)) / ((-1)^(2*n/3) * ((1 + 2*(-1)^(1/3)/3^(2/3))^2 * (1 - (-1)^(2/3)/3^(1/3)))) + Product_{k>=4}((1 - (-1)^(4*k/3)*k/3^(k/3))^(-k)) / ((-1)^(4*n/3) * ((1 + (-1/3)^(1/3)) * (1 - 2*(-1/3)^(2/3))^2))) * 3^(n/3) * n^2 / 54. - _Vaclav Kotesovec_, Apr 24 2017 %F A266941 a(0) = 1 and a(n) = (1/n) * Sum_{k=1..n} (Sum_{d|k} d^(2+k/d)) * a(n-k) for n > 0. - _Seiichi Manyama_, Nov 02 2017 %t A266941 nmax = 40; CoefficientList[Series[Product[1/(1-k*x^k)^k, {k, 1, nmax}], {x, 0, nmax}], x] %t A266941 nmax = 40; s = 1 - x; Do[s *= Sum[Binomial[k, j]*(-1)^j*k^j*x^(j*k), {j, 0, nmax/k}]; s = Expand[s]; s = Take[s, Min[nmax + 1, Exponent[s, x] + 1, Length[s]]];, {k, 2, nmax}]; CoefficientList[Series[1/s, {x, 0, nmax}], x] (* _Vaclav Kotesovec_, Aug 27 2018 *) %Y A266941 Cf. A006906, A265758, A266891, A266942, A266964, A266971. %K A266941 nonn %O A266941 0,3 %A A266941 _Vaclav Kotesovec_, Jan 06 2016