login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A189084 Zero-one sequence based on floor(sqrt(3)): a(A022838(k))=a(k); a(A054406(k))=1-a(k); a(1)=0. 3

%I #4 Mar 30 2012 18:57:23

%S 0,1,1,0,1,0,0,1,1,0,0,0,1,1,1,1,0,0,0,0,0,1,1,1,1,1,1,1,0,0,0,0,0,0,

%T 0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1,1,

%U 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1

%N Zero-one sequence based on floor(sqrt(3)): a(A022838(k))=a(k); a(A054406(k))=1-a(k); a(1)=0.

%e Let u=A022838=(Beatty sequence for sqrt(3)) and v=A054406=(complement of u). Then A189084 is the sequence a given by a(1)=0 and a(u(k))=a(k); a(v(k))=1-a(k).

%t r=3^(1/2); u[n_] := Floor[n*r]; (*A022838*)

%t a[1] = 0; h = 128;

%t c = (u[#1] &) /@ Range[2h];

%t d = (Complement[Range[Max[#1]], #1] &)[c]; (*A054406*)

%t Table[a[d[[n]]] = 1 - a[n], {n, 1, h - 1}]; (*A189084*)

%t Table[a[c[[n]]] = a[n], {n, 1, h}] (*A189084*)

%t Flatten[Position[%, 0]] (*A189085*)

%t Flatten[Position[%%, 1]] (*A189086*)

%Y Cf. A188967, A189085, A189086.

%K nonn

%O 1

%A _Clark Kimberling_, Apr 16 2011

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 09:04 EDT 2024. Contains 371240 sequences. (Running on oeis4.)