login
A187907
Rank transform of the sequence floor((4 - sqrt(5))*n); complement of A187908.
3
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.
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
Sequence in context: A329938 A247429 A187232 * A024806 A084834 A118137
KEYWORD
nonn
AUTHOR
Clark Kimberling, Mar 15 2011
STATUS
approved