|
%I #13 Nov 09 2020 09:15:13
%S 1,4,9,18,31,48,69,98,135,180,233,294,363,440,525,626,743,876,1025,
%T 1190,1371,1568,1781,2010,2255,2516,2793,3086,3395,3720,4061,4434,
%U 4839,5276,5745,6246,6779,7344,7941,8570,9231,9924,10649,11406,12195,13016,13869
%N G.f.: (1 + Sum_{ i >= 0 } 2^i*x^(2^(i+1)-1)) / (1-x)^3.
%H Alois P. Heinz, <a href="/A063916/b063916.txt">Table of n, a(n) for n = 0..10000</a>
%p b:= proc(n) option remember; `if`(n<0, 0, 1+
%p (t-> 2*(b(floor(t))+b(ceil(t))))(n/2-1))
%p end:
%p a:= proc(n) option remember; `if`(n<0, 0, b(n)+a(n-1)) end:
%p seq(a(n), n=0..55); # _Alois P. Heinz_, Jul 10 2019
%t b[n_] := b[n] = If[EvenQ[n], 2b[n/2] + 2b[n/2-1], 4b[(n-1)/2]]+1;
%t b[0] = 0; b[1] = 1; b[2] = 3;
%t b /@ Range[100] // Accumulate (* _Jean-François Alcover_, Nov 09 2020 *)
%Y Partial sums of A063915.
%K nonn
%O 0,2
%A _N. J. A. Sloane_, Sep 01 2001
%E More terms from _Ralf Stephan_, Sep 15 2003
|
