A187897 Rank transform of the sequence A159481; complement of A187898.
1, 2, 3, 5, 6, 8, 10, 11, 12, 13, 15, 16, 18, 19, 20, 21, 22, 24, 25, 26, 28, 29, 30, 32, 34, 35, 36, 38, 39, 40, 42, 43, 44, 46, 47, 48, 50, 51, 52, 54, 56, 57, 58, 60, 61, 63, 64, 65, 67, 68, 69, 71, 72, 74, 75, 76, 77, 79, 81, 82, 84, 85, 86, 87, 88, 90, 91, 92, 94, 95, 96, 98, 100, 101, 102, 103, 104, 106, 108, 109, 111, 112, 113, 114, 115, 117, 119, 120, 121, 123, 124, 125
Offset: 1
Keywords
Programs
-
Mathematica
t = Nest[Flatten[# /. {1 -> {1, 0}, 0 -> {0, 1}}] &, {1}, 8]; t seqA = Table[Sum[t[[k]], {k, 1, n}], {n, 1, 256}] (* A159481 *) 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]] (* A187897 *) Complement[Range[Length[seqA]], limseqU] (* A187898 *) (* by Peter J. C. Moses, Mar 15 2011 *)
Comments