login
Zero-one sequence based on the central polygonal numbers n^2-n+1: a(A002061(k))=a(k); a(A135668(k))=1-a(k), a(1)=0.
3

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

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

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

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

%N Zero-one sequence based on the central polygonal numbers n^2-n+1: a(A002061(k))=a(k); a(A135668(k))=1-a(k), a(1)=0.

%t u[n_] := n^2-n+1; (*A002061*)

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

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

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

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

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

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

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

%Y Cf. A188967, A189136, A189137, A189133, A002061, A135668.

%K nonn

%O 1

%A _Clark Kimberling_, Apr 17 2011