|
%I #12 Feb 12 2018 08:53:12
%S 1,3,4,7,8,9,12,14,15,17,19,21,22,24,26,28,30,31,33,35,37,38,41,42,44,
%T 46,47,50,51,52,55,56,58,60,62,63,65,68,69,71,73,74,76,78,80,82,83,85,
%U 86,89,90,92,94,96,98,99,101,103,104,107,108,111,112,114,116,117,119,121,123,124,126,128,130,132,134,135,137,139,140,142,144,146
%N Rank transform of the sequence floor(n*sqrt(2)-1/2); complement of A187840.
%C See A187224.
%H G. C. Greubel, <a href="/A187839/b187839.txt">Table of n, a(n) for n = 1..10000</a>
%t r=2^(1/2);
%t seqA = Table[Floor[r*n-1/2], {n, 1, 220}]
%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]] (* A187839 *)
%t Complement[Range[Length[seqA]], limseqU] (* A187840 *)
%t (* _Peter J. C. Moses_, Mar 13 2011 *)
%Y Cf. A187224, A187840.
%K nonn
%O 1,2
%A _Clark Kimberling_, Mar 13 2011
|
