%I #14 Oct 05 2017 09:54:24
%S 1,2,1,3,1,4,3,2,1,5,1,6,5,2,1,7,5,3,1,8,7,5,4,2,1,9,7,3,1,10,9,8,7,4,
%T 3,2,1,11,7,5,1,12,11,8,7,6,5,2,1,13,11,9,4,3,5,3,1,14,13,11,4,2,1,15,
%U 13,11,5,3,1,16,15,14,13,12,11,6,5,4,3,2,1
%N Let S be the sequence of rational numbers generated by these rules: 1 is in S, and if u/v is in S (with gcd(u, v) = 1), then (u+1)/v and u/(v+1) are in S, and duplicates are deleted as they occur; a(n) = the numerator of the n-th term of S.
%C See A293248 for the corresponding denominators.
%C The sequence S is a "rational" variant of A232559.
%C If r appears in S, then 1/r appears in S.
%C S is a permutation of the positive rational numbers:
%C - let f be the function u/v -> (u+1)/v
%C and g be the function u/v -> u/(v+1),
%C - let h^k be the k-th iterate of h,
%C - let r = u/v be a rational number in reduced form,
%C - without loss of generality, we can assume that u > v,
%C - according to Dirichlet's theorem on arithmetic progressions, we can choose a prime number p = k*u - 1 > u (where k > 2),
%C - we also have k*u - 1 > k*v,
%C - f^(p-1)(1) = p,
%C - g^(k*v-1)(f^(p-1)(1)) = p / (k*v) (and gcd(p, k*v)=1),
%C - f(g^(k*v-1)(f^(p-1)(1))) = (p+1) / (k*v) = (k*u) / (k*v) = u/v = r, QED.
%H Rémy Sigrist, <a href="/A293247/b293247.txt">Table of n, a(n) for n = 1..10000</a>
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Dirichlet%27s_theorem_on_arithmetic_progressions">Dirichlet's theorem on arithmetic progressions</a>
%H Rémy Sigrist, <a href="/A293247/a293247.png">Colorized scatterplot of A293247 vs. A293248 for n=1..100000</a>
%H Rémy Sigrist, <a href="/A293247/a293247.gp.txt">PARI program for A293247</a>
%e S(1) = 1 by definition; so a(1) = 1.
%e (1+1)/1 = 2 has not yet occurred; so S(2) = 2 and a(2) = 2.
%e 1/(1+1) = 1/2 has not yet occurred; so S(3) = 1/2 and a(3) = 1.
%e (2+1)/1 = 3 has not yet occurred; so S(4) = 3 and a(4) = 3.
%e 2/(1+1) = 1 has already occurred.
%o (PARI) See Links section.
%Y Cf. A232559, A293248.
%K nonn,frac
%O 1,2
%A _Rémy Sigrist_, Oct 03 2017