A187895 Rank transform of the sequence A115384; complement of A187896.
1, 2, 4, 5, 7, 8, 9, 10, 12, 13, 14, 16, 17, 18, 20, 21, 23, 24, 25, 27, 28, 30, 31, 32, 33, 35, 37, 38, 40, 41, 42, 43, 45, 46, 47, 49, 50, 52, 53, 54, 55, 57, 59, 60, 61, 62, 63, 65, 66, 68, 70, 71, 72, 73, 74, 76, 78, 79, 80, 81, 82, 84, 86, 87, 89, 90, 91, 93, 94, 95, 97, 98, 99, 100, 102, 103, 105, 106, 107, 108, 109, 111, 113, 114, 116, 117, 118, 120, 121, 122, 123, 125
Offset: 1
Keywords
Programs
-
Mathematica
t = Nest[Flatten[# /. {0 -> {0, 1}, 1 -> {1, 0}}] &, {0}, 8]; t seqA = Table[Sum[t[[k]], {k, 1, n}], {n, 1, 256}] (* A115384 *) 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]] (* A187895 *) Complement[Range[Length[seqA]], limseqU] (* A187896 *) (* by Peter J. C. Moses, Mar 15 2011 *)
Comments