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 A195983 #20 Nov 19 2014 05:48:10 %S A195983 1,0,1,2,7,32,189,1382,12090,123220,1433760,18750856,272243150, %T A195983 4344694984,75590125814,1423923579636,28872306351182,626986388773656, %U A195983 14518030696953650,357068944503961140,9296168815245928682,255409012110904594768,7385107814017356745382 %N A195983 G.f.: Product_{n>=1} [Product_{k=0..2*n} (1 - k*x)^(-(-1)^k*binomial(2*n,k))]^(1/n). %C A195983 Compare g.f. to the product: %C A195983 _ exp(x) = Product_{n>=1} [Product_{k=0..n} (1 + k*x)^(-(-1)^k*binomial(n,k))]^(1/n), %C A195983 more explicitly: %C A195983 _ exp(x) = [(1+x)^1/(1+0*x)^1] * [(1+x)^2/((1+0*x)^1*(1+2*x)^1)]^(1/2) * [(1+x)^3*(1+3*x)^1/((1+0*x)^1*(1+2*x)^3)]^(1/3) * [(1+x)^4*(1+3*x)^4/((1+0*x)^1*(1+2*x)^6*(1+4*x)^1)]^(1/4) * ... %H A195983 Vaclav Kotesovec, <a href="/A195983/b195983.txt">Table of n, a(n) for n = 0..320</a> %F A195983 G.f.: A(x) = exp( Sum_{n>=1} 2*A000670(n)*x^(n+1)/(n+1) ) where Sum_{n>=0} A000670(n)*x^n = Sum_{n>=0} n!*x^n/Product_{k=0..n} (1-k*x). %F A195983 a(n) ~ n! / (n^2 * (log(2))^n). - _Vaclav Kotesovec_, Nov 19 2014 %e A195983 G.f.: A(x) = 1 + x^2 + 2*x^3 + 7*x^4 + 32*x^5 + 189*x^6 + 1382*x^7 +... %e A195983 where %e A195983 A(x) = [(1-x)^2/((1-0*x)^1*(1-2*x)^1)] * [(1-x)^4*(1-3*x)^4/((1-0*x)^1*(1-2*x)^6*(1-4*x)^1)]^(1/2) * [(1-x)^6*(1-3*x)^20*(1-5*x)^6/((1-0*x)^1*(1-2*x)^15*(1-4*x)^15*(1-6*x)^1)]^(1/3) * ... %e A195983 Also, %e A195983 log(A(x)) = 2*x^2/2 + 6*x^3/3 + 26*x^4/4 + 150*x^5/5 + 1082*x^6/6 + 9366*x^7/7 +...+ 2*A000670(n)*x^(n+1)/(n+1) +... %e A195983 thus, the logarithmic derivative equals the series: %e A195983 A'(x)/A(x) = 2*x/(1-x) + 2*2!*x^2/((1-x)*(1-2*x)) + 2*3!*x^3/((1-x)*(1-2*x)*(1-3*x)) + 2*4!*x^4/((1-x)*(1-2*x)*(1-3*x)*(1-4*x)) +... %o A195983 (PARI) {a(n)=polcoeff(prod(m=1,n,prod(k=1,2*m,(1-k*x+x*O(x^n))^(-(-1)^k*binomial(2*m,k)/m))),n)} %o A195983 (PARI) /* Using o.g.f. exp( Sum_{n>=1} 2*A000670(n)*x^(n+1)/(n+1) ): */ %o A195983 {a(n)=polcoeff(exp(intformal(2*sum(m=1, n+1, m!*x^m/prod(k=1, m, 1-k*x+x*O(x^n))))), n)} %Y A195983 Cf. A000670, A084784. %K A195983 nonn %O A195983 0,4 %A A195983 _Paul D. Hanna_, Sep 25 2011