OFFSET
1
MATHEMATICA
u[n_] := n(n+1)/(2n+1)/6; (*A000330*)
a[1] = 0; a[2]=1; h = 128;
c = (u[#1] &) /@ Range[2h];
d = (Complement[Range[Max[#1]], #1] &)[c];
Table[a[d[[n]]] = 1 - a[n], {n, 1, h - 1}]; (*A189203*)
Table[a[c[[n]]] = a[n], {n, 1, h}] (*A189203*)
Flatten[Position[%, 0]] (*A189204*)
Flatten[Position[%%, 1]] (*A189205*)
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Apr 18 2011
STATUS
approved