|
%I #5 Dec 04 2016 19:46:25
%S 1,2,3,5,7,8,10,11,12,14,15,16,18,19,21,22,23,25,26,28,29,31,32,33,35,
%T 37,38,40,41,42,44,45,47,48,49,51,52,54,56,57,58,59,61,63,64,65,67,68,
%U 70,71,72,74,75,77,78,80,81,82,84,85,86,88,89,91,92,93,95,97,98,100,101,102,104,105,107,108,109,111,112,114,116,117,118,119
%N Rank transform of the sequence floor((-1+sqrt(5))n/2); complement of A187910.
%C See A187224.
%t r=(-1+5^(1/2))/2;
%t seqA = Table[Floor[r*n], {n, 1, 220}] (*A005206*)
%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]] (*A187909*)
%t Complement[Range[Length[seqA]], limseqU] (*A187910*)
%t (*by _Peter J. C. Moses_, Mar 15 2011*)
%Y Cf. A187224, A187910.
%K nonn
%O 1,2
%A _Clark Kimberling_, Mar 15 2011
|
