A187572 Rank transform of the sequence round(n/3); complement of A187473.
1, 2, 3, 4, 6, 7, 8, 10, 11, 12, 14, 15, 16, 18, 19, 20, 21, 22, 23, 25, 26, 27, 29, 30, 31, 33, 34, 35, 36, 37, 38, 40, 41, 42, 44, 45, 46, 48, 49, 50, 51, 52, 53, 55, 56, 57, 59, 60, 61, 63, 64, 65, 66, 67, 68, 70, 71, 72, 74, 75, 76, 78, 79, 80, 82, 83, 84, 86, 87, 88, 90, 91, 92, 93, 94, 95, 97, 98, 99, 101, 102, 103, 105, 106, 107, 108, 109, 110, 112, 113, 114, 116, 117, 118, 120, 121, 122
Offset: 1
Keywords
Programs
-
Mathematica
seqA = Table[Round[n/3], {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]] (*A187572*) Complement[Range[Length[seqA]], limseqU] (*A187573*) (*by Peter J. C. Moses, Mar 11 2011*)
Comments