This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A097631 #17 Feb 24 2019 17:16:26 %S A097631 0,1,5,58,907,19046,496869,15578130,570573623,23929861102, %T A097631 1131235173433,59529368839898,3451899685313523,218712237867226182, %U A097631 15034642075916533997,1114519318895861250082,88631119148029975177327,7526795487859400166772958,679859967684397018073935617 %N A097631 Number of unrooted directed trees on n nodes with a green root. %H A097631 C. Banderier, J.-M. Le Bars and V. Ravelomanana, <a href="https://arxiv.org/abs/math/0411138">Generating functions for kernels of digraphs</a>, arXiv:math/0411138 [math.CO], 2004. %F A097631 a(n) ~ 2^(n-1) * n^(n-2) * (1 - LambertW(1/2)) / (1 + LambertW(1/2)). - _Vaclav Kotesovec_, Feb 24 2019 %t A097631 (* Note: Mathematica's ProductLog is the Lambert W function. *) %t A097631 a[n_] := SeriesCoefficient[-ProductLog[-ProductLog[-2x]/2]/n - ProductLog[-2x] (ProductLog[-2x] + 2)/4, {x, 0, n}] n!; %t A097631 Array[a, 17] (* _Jean-François Alcover_, Feb 24 2019 *) %Y A097631 Equals A097629(n) - A097630(n). %K A097631 nonn %O A097631 1,3 %A A097631 _Ralf Stephan_, Aug 17 2004