OFFSET
1
MATHEMATICA
u[n_] := n(3n-2); (*A000567*)
a[1] = 0; h = 128;
c = (u[#1] &) /@ Range[2h];
d = (Complement[Range[Max[#1]], #1] &)[c]; (*A183220*)
Table[a[d[[n]]] = 1 - a[n], {n, 1, h - 1}]; (*A189219*)
Table[a[c[[n]]] = a[n], {n, 1, h}] (*A189219*)
Flatten[Position[%, 0]] (*A189220*)
Flatten[Position[%%, 1]] (*A189221*)
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Apr 18 2011
STATUS
approved