login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A300994
G.f. A(x) satisfies: [x^n] A(x)^n / (x*A(x))' = 0 for n > 1.
5
1, 1, 1, 2, 10, 78, 787, 9640, 138098, 2257718, 41409511, 841370530, 18753127340, 454916534472, 11932152340884, 336529124983248, 10156870523418506, 326668220673870022, 11154514687950236767, 403044704522955726502, 15364691402538733828594, 616307960867208404852222, 25948893032918425391596907
OFFSET
0,4
COMMENTS
Odd terms seem to occur only at positions 0, 1, and 2*A118113(k) for k >= 0.
Compare to identity: [x^n] (x*F(x))' / F(x)^(n+1) = 0 holds when F(0) = 1.
More generally, [x^n] G(x,k)^(k*(n+1)-1) / (x*G(x,k)^k)' = 0 is satisfied by an integer series G(x,k) when k is a fixed positive integer.
LINKS
FORMULA
G.f. A(x) satisfies: [x^n] A(x)^n / (A(x) + x*A'(x)) = 0 for n > 1.
a(n) ~ c * 2^n * (n-1)!, where c = 0.1261880758068409567445... - Vaclav Kotesovec, Oct 20 2020
EXAMPLE
G.f.: A(x) = 1 + x + x^2 + 2*x^3 + 10*x^4 + 78*x^5 + 787*x^6 + 9640*x^7 + 138098*x^8 + 2257718*x^9 + 41409511*x^10 + ...
such that [x^n] A(x)^n / (x*A(x))' = 0 for n > 1.
ILLUSTRATION OF DEFINITION.
The table of coefficients in A(x)^n / (x*A(x))' begins:
n=0: [1, -2, 1, -4, -29, -306, -3892, -57436, -961833, -17990518, ...];
n=1: [1, -1, 0, -3, -26, -279, -3594, -53588, -904770, -17034879, ...];
n=2: [1, 0, 0, -2, -21, -240, -3196, -48690, -834546, -15887984, ...];
n=3: [1, 1, 1, 0, -13, -185, -2674, -42548, -749180, -14525506, ...];
n=4: [1, 2, 3, 4, 0, -108, -1997, -34928, -646377, -12919990, ...];
n=5: [1, 3, 6, 11, 21, 0, -1124, -25545, -523467, -11040387, ...];
n=6: [1, 4, 10, 22, 54, 152, 0, -14048, -377328, -8851500, ...];
n=7: [1, 5, 15, 38, 104, 366, 1449, 0, -204288, -6313320, ...];
n=8: [1, 6, 21, 60, 177, 666, 3322, 17148, 0, -3380224, ...];
n=9: [1, 7, 28, 89, 280, 1083, 5750, 38090, 240717, 0, ...]; ...
in which the main diagonal consists of all zeros after the initial terms, illustrating that [x^n] A(x)^n / (x*A(x))' = 0 for n > 1.
RELATED SERIES.
(x*A(x))' = 1 + 2*x + 3*x^2 + 8*x^3 + 50*x^4 + 468*x^5 + 5509*x^6 + 77120*x^7 + 1242882*x^8 + 22577180*x^9 + 455504621*x^10 + ...
log(A(x)) = x + x^2/2 + 4*x^3/3 + 33*x^4/4 + 341*x^5/5 + 4252*x^6/6 + 61916*x^7/7 + 1026865*x^8/8 + 19065307*x^9/9 + 391349641*x^10/10 + ...
PROG
(PARI) {a(n) = my(A=[1, 1]); for(i=1, n, A = concat(A, 0); A[#A] = Vec( Ser(A)^(#A-1)/(x*Ser(A))' )[#A]); A[n+1]}
for(n=0, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Mar 29 2018
STATUS
approved