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.

A202668 G.f. satisfies: A(x) = exp( Sum_{n>=1} (A(x) - (-1)^n)^n * x^n/n ).

This page as a plain text file.
%I A202668 #42 Oct 11 2024 05:32:22
%S A202668 1,2,4,12,42,158,618,2498,10360,43832,188420,820608,3613212,16057640,
%T A202668 71933768,324482500,1472604586,6719100254,30804229858,141829955338,
%U A202668 655541387406,3040527731790,14147444737654,66018910398574,308898542610666,1448867831911170
%N A202668 G.f. satisfies: A(x) = exp( Sum_{n>=1} (A(x) - (-1)^n)^n * x^n/n ).
%H A202668 Paul D. Hanna, <a href="/A202668/b202668.txt">Table of n, a(n) for n = 0..800</a>
%H A202668 Vaclav Kotesovec, <a href="/A202668/a202668.txt">Recurrence (of order 12)</a>
%F A202668 G.f. satisfies: A(x) = 1/(1-x*A(x)) * exp( Sum_{n>=1} 1/(1 - (-1)^n*x*A(x))^n * x^n/n ).
%F A202668 G.f. satisfies: A(x) = sqrt( (1 - (A(x)+1)^2*x^2)/(1 - (A(x)-1)^2*x^2) ) / (1 - (A(x)+1)*x).
%F A202668 G.f. satisfies: 0 = -(1+x) - x*A(x) + (1+x)*(1-x)^2*A(x)^2 - x*(1-x)^2*A(x)^3 - x^2*(1+x)*A(x)^4 + x^3*A(x)^5.
%F A202668 From _Vaclav Kotesovec_, Oct 11 2024: (Start)
%F A202668 a(n) ~ sqrt((-1 - s + (-1 - 2*r + 3*r^2)*s^2 + (-1 + 4*r - 3*r^2)*s^3 - r*(2 + 3*r)*s^4 + 3*r^2*s^5)/(1 - r*(1 + 3*s) + r^2*(-1 + 6*s - 6*s^2) + r^3*(1 - 3*s - 6*s^2 + 10*s^3))) / (2*sqrt(Pi) * n^(3/2) * r^(n - 1/2)), where r = 0.20089689587759865228481815120918189691453519374477284069915... and s = 2.3487742728380350386577466365052703249852809669846393564277... are positive real roots of the system of equations s^2*(1 + r^3*(-1 + s)^2*(1 + s)) = 1 + r^2*(-1 + s)^2*s^2 + r*(1 + s + s^2 + s^3) and 2*(-1 + r)^2*(1 + r)*s + 5*r^3*s^4 = r*(1 + 3*(-1 + r)^2*s^2 + 4*r*(1 + r)*s^3).
%F A202668 Numerically, a(n) ~ c * d^n / n^(3/2), where d = 1/r = 4.977677706923229216140896605827075562322447814212438341196056039... and c = 0.7100736662419384614471705442776864037581200760804364785319... (End)
%e A202668 G.f.: A(x) = 1 + 2*x + 4*x^2 + 12*x^3 + 42*x^4 + 158*x^5 + 618*x^6 + ...
%e A202668 where
%e A202668 log(A(x)) = (A(x) + 1)*x + (A(x) - 1)^2*x^2/2 + (A(x) + 1)^3*x^3/3 + (A(x) - 1)^4*x^4/4 + ...
%e A202668 log( A(x)*(1-x*A(x)) ) = 1/(1 + x*A(x))*x + 1/(1 - x*A(x))^2*x^2/2 + 1/(1 + x*A(x))^3*x^3/3 + 1/(1 - x*A(x))^4*x^4/4 + ...
%e A202668 From _Paul D. Hanna_, Oct 11 2024: (Start)
%e A202668 SPECIFIC VALUES.
%e A202668 A(t) = 2 at t = 0.195782060076367892865630673522992184838101...
%e A202668 where 12*t^3 - 4*t^2 - 15*t + 3 = 0.
%e A202668 A(t) = 3/2 at t = 0.1528468026979892250300352740045422934687...
%e A202668 where 45*t^3 - 18*t^2 - 260*t + 40 = 0.
%e A202668 A(1/6) = 1.5975588141693553913621853542774164447766461118908...
%e A202668 A(1/7) = 1.4422077780342017637064340698606478883307441400444...
%e A202668 A(1/8) = 1.3558965312086216338851741626422486193364696459775...
%e A202668 A(1/9) = 1.2992876417963412242026519185070094965390617289384...
%e A202668 A(1/10) = 1.258828814568496961617240364573696812116531654741...
%e A202668 (End)
%o A202668 (PARI) {a(n) = my(A=1+x); for(i=1, n, A=exp(sum(m=1, n, (A - (-1)^m +x*O(x^n))^m * x^m/m))); polcoeff(A, n)}
%o A202668 for(n=0,30,print1(a(n),", "))
%Y A202668 Cf. A202669, A185385, A202630, A202518, A155200.
%K A202668 nonn
%O A202668 0,2
%A A202668 _Paul D. Hanna_, Dec 22 2011