|
%I #33 Jul 12 2018 03:45:31
%S -1,2,-6,36,-300,3270,-43680,691992,-12670560,263281050,-6119720640,
%T 157325242140,-4431909081600,135757694361198,-4492575720622080,
%U 159723265791222000,-6071451523596103680,245720759937001346610,-10548874580411105832960,478801529559871868317140,-22909292930454154076160000,1152457216135660417348971990,-60807227650606789798265487360,3357825559218695417748978547080
%N E.g.f. A(x) satisfies: Sum_{n>=0} 1/n! * Product_{k=0..n} (n+1-k)*x + (k+1)*A(x) = 0.
%C A306065(n-1) = (-1)^n * a(n)/n for n >= 1.
%H Vaclav Kotesovec, <a href="/A306066/b306066.txt">Table of n, a(n) for n = 1..133</a> (terms 1..100 from Paul D. Hanna)
%F E.g.f. A(x) satisfies:
%F (1) Sum_{n>=0} 1/n! * Product_{k=0..n} (n+1-k)*x + (k+1)*A(x) = 0.
%F (2) Sum_{n>=0} 1/n! * Product_{k=0..n} (n+1-k - m)*x + (k+1)*A(x) = -m/x * A(x)^(2*m-2) / (x + A(x))^(m-2).
%F (3) Sum_{n>=0} 1/n! * Product_{k=0..n} (n+1-k)*x + (k+1 - p)*A(x) = -p/A(x) * x^(2*p-2) / (x + A(x))^(p-2).
%F (4) Sum_{n>=0} 1/n! * Product_{k=0..n} (n+1-k - m)*x + (k+1 - p)*A(x) = -(m*x + p*A(x)) * A(x)^(2*m-2) * x^(2*p-2) / (x + A(x))^(m+p-2).
%F (5) A(A(x)) = x.
%F (6) Sum_{n>=0} 1/n! * Product_{k=1..n} (n-k)*x + k*A(x) = -A(x)/x.
%F (7) Sum_{n>=0} 1/n! * Product_{k=1..n} (n-k+1)*x + (k-1)*A(x) = -x/A(x).
%F (8) Sum_{n>=0} 1/(n!*2^n) * Product_{k=1..n} (2*(n-k)+1)*x + (2*k-1)*A(x) = 1.
%F (9) Sum_{n>=0} 1/(n!*2^n) * Product_{k=1..n} (2*(n-k)+1 - m)*x + (2*k-1)*A(x) = (-A(x)/x)^( m*x/(x - A(x)) ).
%F (10) Sum_{n>=0} 1/(n!*2^n) * Product_{k=1..n} (2*(n-k)+1)*x + (2*k-1 - p)*A(x) = (-A(x)/x)^( p*A(x)/(x - A(x)) ).
%F (11) Sum_{n>=0} 1/(n!*2^n) * Product_{k=1..n} (2*(n-k)+1 - m)*x + (2*k-1 - p)*A(x) = (-A(x)/x)^( (m*x + p*A(x))/(x - A(x)) ).
%F a(n)/n! ~ (-1)^n * c * d^n / n^(3/2), where d = 2.45598128882155545489... and c = 0.2658048623687886... - _Vaclav Kotesovec_, Jul 12 2018
%e E.g.f.: A(x) = -x + 2*x^2/2! - 6*x^3/3! + 36*x^4/4! - 300*x^5/5! + 3270*x^6/6! - 43680*x^7/7! + 691992*x^8/8! - 12670560*x^9/9! + 263281050*x^10/10! - 6119720640*x^11/11! + 157325242140*x^12/12! - 4431909081600*x^13/13! + 135757694361198*x^14/14! - 4492575720622080*x^15/15! + 159723265791222000*x^16/16! + ...
%e such that
%e 0 = (x + A(x)) + (2*x + A(x))*(x + 2*A(x))/1! + (3*x + A(x))*(2*x + 2*A(x))*(x + 3*A(x))/2! + (4*x + A(x))*(3*x + 2*A(x))*(2*x + 3*A(x))*(x + 4*A(x))/3! + (5*x + A(x))*(4*x + 2*A(x))*(3*x + 3*A(x))*(2*x + 4*A(x))*(x + 5*A(x))/4! + ...
%e Also,
%e 1 = 1 + (x + A(x))/(1!*2) + (3*x + A(x))*(x + 3*A(x))/(2!*2^2) + (5*x + A(x))*(3*x + 3*A(x))*(x + 5*A(x))/(3!*2^3) + (7*x + A(x))*(5*x + 3*A(x))*(3*x + 5*A(x))*(x + 7*A(x))/(4!*2^4) + ...
%e EXAMPLES OF SUMS.
%e The e.g.f. A(x) satisfies the following sums.
%e (E1) Define
%e S1(m,p) = Sum_{n>=0} 1/n! * Product_{k=0..n} (n+1-k - m)*x + (k+1 - p)*A(x),
%e then
%e S1(m,p) = -(m*x + p*A(x)) * A(x)^(2*m-2) * x^(2*p-2) / (x + A(x))^(m+p-2).
%e (E2) Define
%e S2(m,p) = Sum_{n>=0} 1/(n!*2^n) * Product_{k=1..n} (2*(n-k)+1 - m)*x + (2*k-1 - p)*A(x),
%e then
%e S2(m,p) = (-A(x)/x)^( (m*x + p*A(x))/(x - A(x)) ).
%e Examples of S2(m,p):
%e S2(0,0) = 1,
%e S2(1,-1) = -A(x)/x,
%e S2(-1,1) = -x/A(x),
%e S2(1,2) = (-A(x)/x)^( (x + 2*A(x))/(x - A(x)) ).
%e Examples of S1(m,p):
%e S1(0,0) = 0,
%e S1(1,0) = -(x + A(x)) / x,
%e S1(2,0) = -2*A(x)^2 / x,
%e S1(3,0) = -3*A(x)^4 / (x*(x + A(x))),
%e S1(0,1) = -(x + A(x)) / A(x),
%e S1(0,2) = -2*x^2 / A(x),
%e S1(0,3) = -3*x^4 / (A(x)*(x + A(x))),
%e S1(1,1) = -(x + A(x)),
%e S1(2,2) = -2*x^2*A(x)^2 / (x + A(x)),
%e S1(3,3) = -3*x^4*A(x)^4 / (x + A(x))^3,
%e S1(2,1) = -A(x)^2 * (2*x + A(x)) / (x + A(x)),
%e S1(1,2) = -x^2 * (x + 2*A(x)) / (x + A(x)),
%e S1(4,3) = -x^4*A(x)^6 * (4*x + 3*A(x)) / (x + A(x))^5.
%o (PARI) {a(n) = my(A=[-1]); for(i=0,n, A=concat(A,0); A[#A] = -polcoeff( sum(m=0, #A, 1/m!*prod(k=0, m, (m+1-k)*x + (k+1)*x*Ser(A) ) ), #A)); n!*A[n]}
%o for(n=1,30, print1(a(n),", "))
%Y Cf. A306090, A306067, A306065.
%K sign
%O 1,2
%A _Paul D. Hanna_, Jun 27 2018
| 