login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A218801 G.f. A(x) satisfies: 1+x = Sum_{n>=0} n^n*x^n*A(-n*x)^n/n! * exp(-n*x*A(-n*x)). 0
1, 1, 5, 60, 1539, 82539, 9208086, 2123763546, 1005501106845, 971130952310487, 1903911554135421599, 7548458861108146087406, 60351476147380872012216644, 971005203586845989294297793744, 31389082794561063490845804374258994, 2036350822333032259401319084453988921002 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
COMMENTS
Compare to the identity: 1/(1-x) = Sum_{n>=0} n^n*x^n/n! * exp(-n*x).
LINKS
EXAMPLE
G.f.: A(x) = 1 + x + 5*x^2 + 60*x^3 + 1539*x^4 + 82539*x^5 + 9208086*x^6 +...
where
1+x = 1 + x*A(-x)*exp(-x*A(-x)) + 2^2*x^2*A(-2*x)^2/2!*exp(-2*x*A(-2*x)) + 3^3*x^3*A(-3*x)^3/3!*exp(-3*x*A(-3*x)) + 4^4*x^4*A(-4*x)^4/4!*exp(-4*x*A(-4*x)) +...
PROG
(PARI) a(n)=local(A=1+x+sum(k=2, n-1, a(k)*x^k)+x^2*O(x^n)); if(n<2, 1, -(-1)^n*polcoeff(sum(k=0, n+2, k^k*x^k*subst(A, x, -k*x)^k/k!*exp(-k*x*subst(A, x, -k*x)+x^2*O(x^n))), n+1))
for(n=0, 15, print1(a(n), ", "))
CROSSREFS
Sequence in context: A084939 A171205 A214380 * A303066 A350360 A293454
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Nov 06 2012
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 03:51 EDT 2024. Contains 371264 sequences. (Running on oeis4.)