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 A004113 M1629 #38 Dec 26 2021 20:50:55 %S A004113 1,2,6,18,60,204,734,2694,10162,38982,151920,599244,2389028,9608668, %T A004113 38945230,158904230,652178206,2690598570,11151718166,46412717826, %U A004113 193891596436,812748036380,3417407089470,14410094628558,60920843101858,258169745573158,1096494947168142 %N A004113 Number of rooted trees with n nodes and 2-colored non-leaf nodes. %D A004113 N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence). %H A004113 Alois P. Heinz, <a href="/A004113/b004113.txt">Table of n, a(n) for n = 1..500</a> %H A004113 F. Harary, R. W. Robinson and A. J. Schwenk, <a href="https://doi.org/10.1017/S1446788700016190">Twenty-step algorithm for determining the asymptotic number of trees of various species</a>, J. Austral. Math. Soc., Series A, 20 (1975), 483-503. <a href="https://doi.org/10.1017/S1446788700033760">Errata</a>: Vol. A 41 (1986), p. 325. %H A004113 N. J. A. Sloane, <a href="/transforms.txt">Transforms</a> %H A004113 <a href="/index/Ro#rooted">Index entries for sequences related to rooted trees</a> %H A004113 <a href="/index/Tra#trees">Index entries for sequences related to trees</a> %F A004113 Shifts left and halves under EULER transform. %F A004113 a(n) ~ c * d^n / n^(3/2), where d = 4.49415643203339504537343052838796824... and c = 0.368722987377516657464802259... - _Vaclav Kotesovec_, Feb 28 2014 %p A004113 with(numtheory): etr:= proc(p) local b; b:= proc(n) option remember; `if`(n=0, 1, (add(add(d*p(d), d=divisors(j)) *b(n-j), j=1..n))/n) end end: b:= etr(a): a:= n-> `if`(n<=1, n, 2*b(n-1)): seq(a(n), n=1..30); # _Alois P. Heinz_, Sep 06 2008 %t A004113 etr[p_] := Module[{b}, b[n_] := b[n] = If[n == 0, 1, Sum[Sum[d*p[d], {d, Divisors[j]}]*b[n - j], {j, 1, n}]/n ]; b]; b = etr[a]; a[n_] := If[n <= 1, n, 2*b[n - 1]]; Table[a[n], {n, 1, 27}] (* _Jean-François Alcover_, Jan 29 2013, translated from _Alois P. Heinz_'s Maple program *) %Y A004113 Cf. A004114, A052316, A052317. %K A004113 nonn,eigen %O A004113 1,2 %A A004113 _N. J. A. Sloane_ %E A004113 Extended with better description from _Christian G. Bower_, Apr 15 1998