login
a(0)=0; for n > 0, a(n) is the sum of the partial sums of binary digits of a(n-1) + 1, starting from the least significant digit, and given weights 1, 2, 2^2, etc.
0

%I #21 Aug 31 2017 14:20:22

%S 0,1,2,5,10,37,154,1125,9114,121957,1188762,19782757,449979290,

%T 7603084389,147015738266,4800786586725,197509352924058,

%U 6557890088131685,254650888299357082,7815799643948571749,315002221645968581530

%N a(0)=0; for n > 0, a(n) is the sum of the partial sums of binary digits of a(n-1) + 1, starting from the least significant digit, and given weights 1, 2, 2^2, etc.

%e a(4) is computed as follows: a(3) + 1 = 6 = 110_2; the partial sums of digits and their weights are (0,1), (0+1,2), and (0+1+1,4), so the sum of partial sums is 0*1 + 1*2 + 2*4 = 10.

%t PartialSums = Function[Accumulate[Reverse[IntegerDigits[#, 2]]]]

%t NestList[Total[# Power[2, Range[0, Length[#] - 1]]] &[PartialSums[1 + #]] &, 0, 20]

%K base,nonn,easy

%O 0,3

%A _Emanuel Landeholm_, Jul 18 2017