login
a(1) = 1, a(n) = 2*a(n-1) + a([n/2]).
4

%I #10 Mar 30 2012 18:38:42

%S 1,3,7,17,37,81,169,355,727,1491,3019,6119,12319,24807,49783,99921,

%T 200197,401121,802969,1607429,3216349,6435717,12874453,25755025,

%U 51516169,103044657,206101633,412228073,824480953,1649011689

%N a(1) = 1, a(n) = 2*a(n-1) + a([n/2]).

%p A033489 := proc(n) option remember; if n = 1 then 1 else A033489(n-1)+A033489(round(2*(n-1)/2))+A033489(round((n-1)/2)); fi; end;

%o (PARI) a(n)=if(n<2,1,(2*a(n-1)+a(floor(n/2))))

%K nonn

%O 1,2

%A _N. J. A. Sloane_.

%E Better description and more terms from _Benoit Cloitre_, Jan 06 2004