OFFSET
1,3
COMMENTS
Conjectured to consist entirely of integers.
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..300
FORMULA
The g.f. of A266489, G(x) = 1 + x*A'(x), satisfies: [x^n] G( x/G(x)^n ) = 0 for n>1.
EXAMPLE
G.f.: A(x) = x + x^2 + 4*x^3 + 33*x^4 + 436*x^5 + 8183*x^6 + 204086*x^7 + 6482641*x^8 + 254507098*x^9 + 12071123966*x^10 +...
such that G(x) = 1 + x*A'(x):
G(x) = 1 + x + 2*x^2 + 12*x^3 + 132*x^4 + 2180*x^5 + 49098*x^6 + 1428602*x^7 + 51861128*x^8 + 2290563882*x^9 + 120711239660*x^10 +...+ A266489(n)*x^n +...
satisfies: [x^n] G( x/G(x)^n ) = 0 for n>1.
PROG
(PARI) {a(n) = my(A=[1, 1]); for(i=2, n, A=concat(A, 0); A[#A] = -Vec(subst(Ser(A), x, x/Ser(A)^(#A-1)))[#A]); A[n+1]/n}
for(n=1, 30, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Feb 08 2016
STATUS
approved