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

1, 3, 5, 7, 9, 11, 13, 15, 16, 19, 20, 23, 24, 26, 28, 30, 32, 34, 36, 38, 40, 41, 43, 46, 48, 49, 51, 53, 55, 57, 59, 61, 63, 64, 66, 68, 71, 73, 74, 76, 78, 80, 82, 84, 86, 88, 89, 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

Offset

1,2

Comments

See A187224. A187232(n) = A187907(n) for n=1..20; A187232(21)=39 and A187907(21)=40.

Links

Table of n, a(n) for n=1..81.

Mathematica

r=4-5^(1/2);
seqA = Table[Floor[r*n], {n, 1, 220}] (* A187330 *)
seqB = Table[n, {n, 1, 220}];  (* A000027 *)
jointRank[{seqA_,
   seqB_}] := {Flatten@Position[#1, {_, 1}],
    Flatten@Position[#1, {_, 2}]} &[
  Sort@Flatten[{{#1, 1} & /@ seqA, {#1, 2} & /@ seqB}, 1]];
limseqU =
FixedPoint[jointRank[{seqA, #1[[1]]}] &,
   jointRank[{seqA, seqB}]][[1]] (* A187907 *)
Complement[Range[Length[seqA]], limseqU]  (* A187908 *)
(* Peter J. C. Moses, Mar 15 2011 *)

Prog

(Maxima) makelist(floor((4-sqrt(5))*n), n, 1, 100); /* Martin Ettl, Oct 17 2012 */

Crossrefs

Cf. A187224, A187908.

Sequence in context: A329938 A247429 A187232 * A024806 A084834 A118137

Adjacent sequences: A187904 A187905 A187906 * A187908 A187909 A187910

Keyword

nonn

Author

Clark Kimberling, Mar 15 2011

Status

approved