A353708 - OEIS

%I #23 May 06 2022 19:48:27

%S 0,1,2,4,5,3,8,12,6,16,9,7,18,24,13,32,34,10,17,20,14,11,33,36,22,19,

%T 40,44,21,64,42,15,65,48,26,66,37,25,72,38,23,73,96,50,27,68,100,35,

%U 128,28,29,67,98,52,129,74,30,49,97,70,130,41,45,80,82,39,132,88,43,131,84,56,136,69,51,58,76,133,144,90,46,160,81,31,134,192,57,47

%N a(0)=0, a(1)=1; thereafter a(n) = smallest nonnegative integer not among the earlier terms of the sequence such that a(n) and a(n-2) have no common 1-bits in their binary representations.

%C A set-theory analog of A121216.

%C This is a permutation of the nonnegative numbers.

%H Michael De Vlieger, <a href="/A353708/b353708.txt">Table of n, a(n) for n = 0..16384</a> (a(n) for n = 0..748 by N. J. A. Sloane)

%H <a href="/index/Per#IntegerPermutation">Index entries for sequences that are permutations of the nonnegative integers</a>

%p read(transforms) : # ANDnos def'd here

%p A353708 := proc(n)

%p     option remember;

%p     local c, i, known ;

%p     if n <= 2 then

%p         n;

%p     else

%p         for c from 1 do

%p             known := false ;

%p             for i from 1 to n-1 do

%p                 if procname(i) = c then

%p                     known := true;

%p                     break ;

%p                 end if;

%p             end do:

%p             if not known and ANDnos(c, procname(n-2)) =0 then

%p                 return c;

%p             end if;

%p         end do:

%p     end if;

%p end proc: # Following R. J. Mathar's program for A109812.

%p [seq(A353708(n), n=0..256)] ;

%p # second Maple program:

%p b:= proc() false end: t:= 2:

%p a:= proc(n) option remember; global t; local k; if n<2 then n

%p       else for k from t while b(k) or Bits[And](k, a(n-2))>0

%p       do od; b(k):=true; while b(t) do t:=t+1 od; k fi

%p     end:

%p seq(a(n), n=0..100);  # _Alois P. Heinz_, May 06 2022

%t nn = 87; c[_] = -1; a[0] = c[0] = 0; a[1] = c[1] = 1; u = 2; Do[k = u; While[Nand[c[k] == -1, BitAnd[a[n - 2], k] == 0], k++]; Set[{a[n], c[k]}, {k, n}]; If[k == u, While[c[u] > -1, u++]], {n, 2, nn}]; Array[a, nn + 1, 0] (* _Michael De Vlieger_, May 06 2022 *)

%Y Cf. A109812, A121216.

%K nonn

%O 0,3

%A _N. J. A. Sloane_, May 06 2022