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
Keywords
Programs
-
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 *)
-
Maxima
makelist(floor((4-sqrt(5))*n),n,1,100); /* Martin Ettl, Oct 17 2012 */
Comments