OFFSET
1,1
COMMENTS
See A187224.
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..10000
MATHEMATICA
r=2^(1/2);
seqA = Table[Floor[r*n-1/2], {n, 1, 220}]
seqB = Table[n, {n, 1, 220}]; (* A000027 *)
jointRank[{seqA_,
seqB_}] := {Flatten@Position[#1, {_, 1}],
Flatten@Position[#1, {_, 2}]} &[
Sort@Flatten[{{#1, 1} & /@ seqA, {#1, 2} & /@ seqB}, 1]];
limseqU =
FixedPoint[jointRank[{seqA, #1[[1]]}] &,
jointRank[{seqA, seqB}]][[1]] (* A187839 *)
Complement[Range[Length[seqA]], limseqU] (* A187840 *)
(* Peter J. C. Moses, Mar 13 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Mar 13 2011
STATUS
approved