|
%I #14 Oct 17 2016 16:40:43
%S 1,2,2,3,4,3,4,4,6,5,8,4,6,5,8,5,8,7,12,6,10,9,16,5,8,7,12,6,10,9,16,
%T 6,10,9,16,8,14,13,24,7,12,11,20,10,18,17,32,6,10,9,16,8,14,13,24,7,
%U 12,11,20,10,18,17,32,7,12,11,20,10,18,17,32,9,16,15,28,14,26,25,48,8,14,13,24
%N Guy Steele's sequence GS(4,5) (see A135416).
%H G. C. Greubel, <a href="/A135529/b135529.txt">Table of n, a(n) for n = 1..1000</a>
%F a(n) = A135533(n) + 1 - 2^(A000120(n)-1). - _Don Knuth_, Mar 01 2008
%p GS(4,5,200); # [see A135416].
%t i = 4; j = 5; Clear[a]; a[1] = 1; a[n_?EvenQ] := a[n] = {0, 1, a[n/2], a[n/2]+1, 2*a[n/2], 2*a[n/2]+1}[[i]]; a[n_?OddQ] := a[n] = {0, 1, a[(n-1)/2], a[(n-1)/2]+1, 2*a[(n-1)/2], 2*a[(n-1)/2]+1}[[j]]; Array[a, 83] (* _Jean-François Alcover_, Sep 12 2013 *)
%Y Cf. A135416.
%K nonn
%O 1,2
%A _N. J. A. Sloane_, based on a message from Guy Steele and _Don Knuth_, Mar 01 2008
|
