%I #23 Jan 04 2022 02:12:33
%S 0,1,2,4,6,9,15,25,35,49,70,100,160,256,416,676,936,1296,1800,2500,
%T 3550,5041,7171,10201,16261,25921,41377,66049,107169,173889,282309,
%U 458329,634349,877969,1215289,1682209,2335897,3243601,4504301,6255001,8881051,12609601
%N Number of subtrees of a complete binary tree.
%C Take the complete binary tree with n labeled nodes. Here is a poor picture of the tree with 6 nodes:
%C R
%C / \
%C / \
%C / \
%C o o
%C / \ /
%C o o o
%C Then the number of rooted subtrees of the graph is a(n).
%H A. V. Aho and N. J. A. Sloane, <a href="https://www.fq.math.ca/Scanned/11-4/aho-a.pdf">Some doubly exponential sequences</a>, Fibonacci Quarterly, Vol. 11, No. 4 (1973), pp. 429-437, <a href="http://neilsloane.com/doc/doubly.html">alternative link</a>.
%H <a href="/index/Ro#rooted">Index entries for sequences related to rooted trees</a>
%F a(0) = 0, a(1) = 1.
%F a(n) = 1 + a(floor((n-1)/2)) + a(ceiling((n-1)/2)) + a(floor((n-1)/2)) * a(ceiling((n-1)/2)) = (1+a(floor((n-1)/2))) * (1+a(ceiling((n-1)/2))).
%F If b(n) is sequence A005468, then a(n)=b(n+1)-1. From the definition of A005468, a(n) = b(floor((n+1)/2)*b(ceiling((n+1)/2). So for every odd n a(n) is a square: a(2n-1)=b(n)^2.
%F If c(n) is sequence A004019, then c(n)=a(2^n-1).
%F A004019 (and Aho and Sloane) give a closed formula for c(n) that translates to a(n) = nearest integer to b^((n+1)/2) - 1" where b = 2.25851...; the formula gives the asymptotic behavior of this sequence, however it does not compute the correct values for a(n) unless n+1 is a power of two.
%e For n = 3, the a(3) = 4 subtrees are:
%e R R R R
%e / \ / \
%e o o o o
%p a:= proc(n) option remember; `if`(n<2, n,
%p (h-> (1+a(h))*(1+a(n-1-h)))(iquo(n, 2)))
%p end:
%p seq(a(n), n=0..50); # _Alois P. Heinz_, Jan 02 2022
%t a[0] = 0; a[1] = 1; a[n_?EvenQ] := a[n] = (1 + a[n/2 - 1])*(1 + a[n/2]); a[n_?OddQ] := a[n] = (1 + a[(n-1)/2])^2; Table[a[n], {n, 0, 32}] (* _Jean-François Alcover_, Oct 19 2011 *)
%Y Cf. A004019, A005468.
%K easy,nice,nonn
%O 0,3
%A _Paolo Bonzini_, Mar 04 2009
|