login
Rank transform of the sequence A115384; complement of A187896.
3

%I #10 Dec 04 2016 19:46:25

%S 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,

%T 35,37,38,40,41,42,43,45,46,47,49,50,52,53,54,55,57,59,60,61,62,63,65,

%U 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

%N Rank transform of the sequence A115384; complement of A187896.

%C A187895 is the rank transform (see A187224) of the sequence of partial sums of the Thue-Morse sequence A010060.

%t t = Nest[Flatten[# /. {0 -> {0, 1}, 1 -> {1, 0}}] &, {0}, 8]; t

%t seqA = Table[Sum[t[[k]], {k, 1, n}], {n, 1, 256}] (* A115384 *)

%t seqB = Table[n, {n, 1, 220}]; (* A000027 *)

%t jointRank[{seqA_,

%t seqB_}] := {Flatten@Position[#1, {_, 1}],

%t Flatten@Position[#1, {_, 2}]} &[

%t Sort@Flatten[{{#1, 1} & /@ seqA, {#1, 2} & /@ seqB}, 1]];

%t limseqU =

%t FixedPoint[jointRank[{seqA, #1[[1]]}] &,

%t jointRank[{seqA, seqB}]][[1]] (* A187895 *)

%t Complement[Range[Length[seqA]], limseqU] (* A187896 *)

%t (* by _Peter J. C. Moses_, Mar 15 2011 *)

%Y Cf. A187224, A187896, A187897.

%K nonn

%O 1,2

%A _Clark Kimberling_, Mar 15 2011