OFFSET
1,2
MATHEMATICA
r=4-5^(1/2);
seqA = Table[Floor[r*n], {n, 1, 220}] (* A187330 *)
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]] (* A187907 *)
Complement[Range[Length[seqA]], limseqU] (* A187908 *)
(* Peter J. C. Moses, Mar 15 2011 *)
PROG
(Maxima) makelist(floor((4-sqrt(5))*n), n, 1, 100); /* Martin Ettl, Oct 17 2012 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Mar 15 2011
STATUS
approved