cp's OEIS Frontend

A187907 Rank transform of the sequence floor((4 - sqrt(5))*n); complement of A187908.

%I A187907 #14 Sep 09 2018 02:23:25
%S A187907 1,3,5,7,9,11,13,15,16,19,20,23,24,26,28,30,32,34,36,38,40,41,43,46,
%T A187907 48,49,51,53,55,57,59,61,63,64,66,68,71,73,74,76,78,80,82,84,86,88,89,
%U A187907 91,93,95,97,99,101,103,105,107,109,111,113,114,116,118,120,122,124,126,128,130,132,134,136,138,139,141,143,145,147,149,151,153,155
%N A187907 Rank transform of the sequence floor((4 - sqrt(5))*n); complement of A187908.
%C A187907 See A187224. A187232(n) = A187907(n) for n=1..20; A187232(21)=39 and A187907(21)=40.
%t A187907 r=4-5^(1/2);
%t A187907 seqA = Table[Floor[r*n], {n, 1, 220}] (* A187330 *)
%t A187907 seqB = Table[n, {n, 1, 220}];  (* A000027 *)
%t A187907 jointRank[{seqA_,
%t A187907    seqB_}] := {Flatten@Position[#1, {_, 1}],
%t A187907     Flatten@Position[#1, {_, 2}]} &[
%t A187907   Sort@Flatten[{{#1, 1} & /@ seqA, {#1, 2} & /@ seqB}, 1]];
%t A187907 limseqU =
%t A187907 FixedPoint[jointRank[{seqA, #1[[1]]}] &,
%t A187907    jointRank[{seqA, seqB}]][[1]] (* A187907 *)
%t A187907 Complement[Range[Length[seqA]], limseqU]  (* A187908 *)
%t A187907 (* _Peter J. C. Moses_, Mar 15 2011 *)
%o A187907 (Maxima) makelist(floor((4-sqrt(5))*n),n,1,100); /* _Martin Ettl_, Oct 17 2012 */
%Y A187907 Cf. A187224, A187908.
%K A187907 nonn
%O A187907 1,2
%A A187907 _Clark Kimberling_, Mar 15 2011