%I #7 Mar 30 2012 18:57:24
%S 0,1,0,1,1,0,0,1,1,0,0,1,1,0,0,0,1,1,1,0,0,0,1,1,1,0,0,0,1,1,1,1,0,0,
%T 0,0,1,1,1,1,0,0,0,0,1,1,1,1,0,0,0,0,1,1,1,1,1,0,0,0,0,0,1,1,1,1,1,0,
%U 0,0,0,0,1,1,1,1,1,0,0,0,0,0,1,1,1,1,1,0,0,0,0,0,0,1,1,1,1,1,1,0,0,0,0,0,0,1,1,1,1,1,1,0,0,0,0,0,0,1,1,1,1,1,1,0,0,0,0,0,0,1,1,1,1,1,1
%N Zero-one sequence based on the sequence n(n+1)(2n+1)/6: a(A000330(k))=a(k); a(v(k))=1-a(k), where v=complement(A000330) and a(1)=0.
%t u[n_] := n(n+1)/(2n+1)/6; (*A000330*)
%t a[1] = 0; a[2]=1; h = 128;
%t c = (u[#1] &) /@ Range[2h];
%t d = (Complement[Range[Max[#1]], #1] &)[c];
%t Table[a[d[[n]]] = 1 - a[n], {n, 1, h - 1}]; (*A189203*)
%t Table[a[c[[n]]] = a[n], {n, 1, h}] (*A189203*)
%t Flatten[Position[%, 0]] (*A189204*)
%t Flatten[Position[%%, 1]] (*A189205*)
%Y Cf. A188967, A189204, A189205.
%K nonn
%O 1
%A _Clark Kimberling_, Apr 18 2011
|