|
%I #27 Jul 02 2020 15:46:50
%S 2,3,5,7,11,16,24,36,54,81,121,181,271,406,609,913,1369,2053,3079,
%T 4618,6927,10390,15585,23377,35065,52597,78895,118342,177513,266269,
%U 399403,599104,898656,1347984,2021976,3032964,4549446,6824169,10236253,15354379,23031569,34547353,51821029,77731543,116597315
%N a(n) = floor(a(n-1)*3/2) bitwise-OR LSB(a(n-1)*n), with a(1)=2 (LSB = Least Significant Bit).
%F a(n) = floor(a(n-1)*3/2) bitwise-OR LSB(a(n-1)*n), with a(1)=2 (LSB = Least Significant Bit).
%e even index = x + (x>>1)
%e odd index = ((x + (x>>1)) | least significant bit of x
%e 3 = 2 + 1
%e 5 = (3 + 1) | 1
%e 7 = 5 + 2
%e 11 = (7 + 3) | 1
%e 16 = 11 + 5
%e 24 = (16 + 8) | 0
%e 36 = 24 + 12
%e 54 = (36 + 18) | 0
%e 81 = 54 + 27
%e 121 = (81 + 40) | 1
%t a[1] = 2; a[n_] := a[n] = BitOr[Floor[a[n - 1]*3/2], Mod[a[n - 1]*n, 2]]; Array[a, 45] (* _Amiram Eldar_, May 01 2020 *)
%o (Ruby) values = [2]; 100.times { |index| x = values.last; values << ((x + (x>>1)) | (x&1&index)) }; p values
%o (PARI) lista(nn) = {my(a=2, va = List(a)); for (n=2, nn, my(x = a + a>>1); if (n%2, a = bitor(x, a%2), a = x); listput(va, a);); Vec(va);} \\ _Michel Marcus_, May 03 2020
%Y Cf. A112088, A333588.
%K nonn,base
%O 1,1
%A _Simon Strandgaard_, Mar 27 2020
|
