cp's OEIS Frontend

A187841 Rank transform of the sequence floor(n*sqrt(2)+1/2); complement of A187842.

%I A187841 #5 Dec 04 2016 19:46:25
%S A187841 1,3,5,7,8,10,12,14,16,17,19,21,23,25,26,28,30,31,34,35,37,39,41,42,
%T A187841 44,46,48,50,51,53,55,57,59,60,62,64,66,68,69,71,73,74,77,78,80,82,83,
%U A187841 85,87,89,91,93,95,96,98,100,102,103,105,107,109,111,112,114,116,117,119,121,123,125,126,128,130,132,134,135,137,139,141,143,145,146
%N A187841 Rank transform of the sequence floor(n*sqrt(2)+1/2); complement of A187842.
%C A187841 See A187224.
%t A187841 r=2^(1/2);
%t A187841 seqA = Table[Floor[r*n+1/2], {n, 1, 220}]
%t A187841 seqB = Table[n, {n, 1, 220}];(*A000027*)
%t A187841 jointRank[{seqA_,
%t A187841    seqB_}] := {Flatten@Position[#1, {_, 1}],
%t A187841     Flatten@Position[#1, {_, 2}]} &[
%t A187841   Sort@Flatten[{{#1, 1} & /@ seqA, {#1, 2} & /@ seqB}, 1]];
%t A187841 limseqU =
%t A187841  FixedPoint[jointRank[{seqA, #1[[1]]}] &,
%t A187841    jointRank[{seqA, seqB}]][[1]] (*A187841*)
%t A187841 Complement[Range[Length[seqA]], limseqU]  (*A187842*)
%t A187841 (*by _Peter J. C. Moses_, Mar 13 2011*)
%Y A187841 Cf. A187224, A187842, A187839.
%K A187841 nonn
%O A187841 1,2
%A A187841 _Clark Kimberling_, Mar 13 2011