|
%I #81 May 02 2024 12:45:56
%S 0,1,2,5,10,4,3,12,17,6,9,18,8,7,24,33,14,32,11,36,19,40,16,13,48,15,
%T 80,34,20,35,28,65,22,41,66,21,42,68,25,38,72,23,64,26,69,50,73,52,67,
%U 44,81,46,129,30,97,130,29,98,132,27,100,131,56,70,49,74,37,82
%N a(n) = n if n <= 2, otherwise the smallest number not occurring earlier having in its binary representation at least one bit in common with a(n-2), but none with a(n-1).
%C Conjectured to be a permutation of the nonnegative integers. [Comment modified by _N. J. A. Sloane_, Jan 10 2015]
%C This is a purely set-based version of A098550, using the binary representation of numbers.
%C An equivalent definition in terms of sets: S(0) = {}, S(1) = {0}, S(2} = {1}; thereafter S(n) is the smallest set (different from the S{i} with i < n) of positive integers such that S(n) meets S(n-2) but is disjoint from S(n-1). - _N. J. A. Sloane_, Mar 27 2022; corrected Aug 01 2022.
%H Chai Wah Wu, <a href="/A252867/b252867.txt">Table of n, a(n) for n = 0..50002</a> (First 10000 terms from Reinhard Zumkeller)
%H David L. Applegate, Hans Havermann, Bob Selcoe, Vladimir Shevelev, N. J. A. Sloane, and Reinhard Zumkeller, <a href="http://arxiv.org/abs/1501.01669">The Yellowstone Permutation</a>, arXiv preprint arXiv:1501.01669 [math.NT], 2015 and <a href="https://cs.uwaterloo.ca/journals/JIS/VOL18/Sloane/sloane9.html">J. Int. Seq. 18 (2015) 15.6.7</a>.
%H Chai Wah Wu, <a href="/A252867/a252867.png">Scatterplot of first million terms</a>
%H Chai Wah Wu, <a href="/A252867/a252867_1.png">Scatterplot of first million terms</a>, with red lines powers of 2.
%H Chai Wah Wu, <a href="/A252867/a252867.zip.txt">Gzipped file with first million terms</a> [Save file, delete .txt suffix, then open]
%e The sequence of sets is {}, {0}, {1}, {0,2}, {1,3}, {2}, {0,1}, {2,3}. After the initial 3 terms, S(n) is the minimum set (as ordered by A048793) that has a nonempty intersection with S(n-2) but empty intersection with S(n-1). [Typos corrected by _N. J. A. Sloane_, Aug 01 2022 at the suggestion of _Michel Dekking_.]
%e Comment from _N. J. A. Sloane_, Dec 31 2014: The binary expansions of the first few terms are:
%e 0 = 000000
%e 1 = 000001
%e 2 = 000010
%e 5 = 000101
%e 10 = 001010
%e 4 = 000100
%e 3 = 000011
%e 12 = 001100
%e 17 = 010001
%e 6 = 000110
%e 9 = 001001
%e 18 = 010010
%e 8 = 001000
%e 7 = 000111
%e 24 = 011000
%e 33 = 100001
%e 14 = 001110
%e 32 = 100000
%e 11 = 001011
%e 36 = 100100
%e 19 = 010011
%e 40 = 101000
%e ...
%p read("transforms") : # define ANDnos
%p A252867 := proc(n)
%p local a,known,i ;
%p option remember;
%p if n <=2 then
%p n;
%p else
%p for a from 3 do
%p known := false ;
%p for i from 1 to n-1 do
%p if procname(i) = a then
%p known := true;
%p break;
%p end if;
%p end do:
%p if not known then
%p if ANDnos(a, procname(n-1)) = 0 and ANDnos(a,procname(n-2)) > 0 then
%p return a;
%p end if;
%p end if;
%p end do:
%p end if
%p end proc:
%p seq(A252867(n),n=0..200) ; # _R. J. Mathar_, May 02 2024
%t a[n_] := a[n] = If[n<3, n, For[k=3, True, k++, If[FreeQ[Array[a, n-1], k], If[BitAnd[k, a[n-2]] >= 1 && BitAnd[k, a[n-1]] == 0, Return[k]]]]];
%t Table[a[n], {n, 0, 100}] (* _Jean-François Alcover_, Oct 03 2018 *)
%o (PARI) invecn(v,k,x)=for(i=1,k,if(v[i]==x,return(i)));0
%o alist(n)=local(v=vector(n,i,i-1), x); for(k=4, n, x=3; while(invecn(v, k-1, x)||!bitand(v[k-2], x)||bitand(v[k-1],x), x++); v[k]=x); v
%o (Haskell)
%o import Data.Bits ((.&.)); import Data.List (delete)
%o a252867 n = a252867_list !! n
%o a252867_list = 0 : 1 : 2 : f 1 2 [3..] where
%o f :: Int -> Int -> [Int] -> [Int]
%o f u v ws = g ws where
%o g (x:xs) = if x .&. u > 0 && x .&. v == 0
%o then x : f v x (delete x ws) else g xs
%o -- _Reinhard Zumkeller_, Dec 24 2014
%o (Python)
%o A252867_list, l1, l2, s, b = [0,1,2], 2, 1, 3, set()
%o for _ in range(10**2):
%o i = s
%o while True:
%o if not (i in b or i & l1) and i & l2:
%o A252867_list.append(i)
%o l2, l1 = l1, i
%o b.add(i)
%o while s in b:
%o b.remove(s)
%o s += 1
%o break
%o i += 1 # _Chai Wah Wu_, Dec 27 2014
%Y Cf. A098550, A252865, A048793, A252868.
%Y Reading this sequence mod 2 gives A253050 and A253051.
%Y Cf. A253581, A253582, A253589 (binary weight), A253603.
%Y Analyzed further in A303596, A303597, A303598, A303599, A305368.
%Y The graphs of A109812, A252867, A305369, A305372 all have roughly the same, mysterious, fractal-like structure. - _N. J. A. Sloane_, Jun 03 2018
%K nonn,look,base
%O 0,3
%A _Franklin T. Adams-Watters_, Dec 23 2014
|
