Search: id:A007889 Results 1-1 of 1 results found. %I A007889 %S A007889 1,1,2,7,36,246,2104,21652,260720,3598120,56010096,971055240,18558391936, %T A007889 387665694976,8787898861568,214868401724416,5636819806209792, %U A007889 157935254554567296,4707152127520549120,148704074888134683520 %N A007889 Number of intransitive (or alternating) trees: vertices are [0,n] and for no i home page %H A007889 Index entries for sequences related to trees %F A007889 E.g.f. satisfies: A(x) = exp( x*(1 + A(x))/2 ). E.g.f. A(x) equals the inverse function of 2*log(x)/(1+x). - Paul D. Hanna (pauldhanna(AT)juno.com), Mar 29 2008 %F A007889 E.g.f.: -2/x*LambertW(-1/2*x*exp(1/2*x)). - Vladeta Jovovic (vladeta(AT)eunet.rs), Mar 29 2008 %F A007889 Comments from Vladeta Jovovic and Paul D. Hanna, Apr 03 2008 (Start): Powers of e.g.f.: If A(x)^p = Sum_{n>=0} a(n,p)*x^n/n! then a(n,p) = (1/2^n)* Sum_{k=0..n} binomial(n,k)*p*(k+p)^(n-1). %F A007889 Let A(x) = e.g.f. of A007889, B(x) = e.g.f. of A138860 where B(x) = exp( x*[B(x) + B(x)^2]/2 ); then B(x) = A(x*B(x)) = (1/x)*Series_Reversion(x/ A(x)) and A(x) = B(x/A(x)) = x/Series_Reversion(x*B(x)). (End) %p A007889 f := (n)->1/(2^n*(n+1))*sum(binomial(n+1,k)*k^n,'k'=1..(n+1)); %o A007889 (PARI) {a(n)=local(A=1+x);for(i=0,n,A=exp(x*(1+A)/2 +x*O(x^n)));n!*polcoeff(A, n)} - Paul D. Hanna (pauldhanna(AT)juno.com), Mar 29 2008 %o A007889 (PARI) /* Coefficients of A(x)^p are given by: */ {a(n,p=1)=(1/2^n)*sum(k=0, n,binomial(n,k)*p*(k+p)^(n-1))} - Vladeta Jovovic and Paul D. Hanna, Apr 03 2008 %Y A007889 Cf. A038049. %Y A007889 Cf. A138860. %Y A007889 Sequence in context: A095793 A029768 A167199 this_sequence A125033 A034430 A143805 %Y A007889 Adjacent sequences: A007886 A007887 A007888 this_sequence A007890 A007891 A007892 %K A007889 nonn,easy,nice %O A007889 0,3 %A A007889 Alexander Postnikov [ apost(AT)math.mit.edu ] Search completed in 0.001 seconds