login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A161893
Denominators of S(n) = Sum_{j=2..n} (pi((j+1)^2) - pi(j^2))/(pi((j+1)^2)*pi(j^2)) where pi(k) = A000720(k).
1
4, 3, 18, 22, 30, 9, 11, 50, 15, 17, 78, 44, 48, 27, 122, 33, 72, 39, 170, 92, 198, 210, 57, 61, 258, 274, 73, 77, 81, 172, 362, 382, 200, 105, 438, 228, 240, 502, 526, 137, 566, 590, 153, 638, 658, 171, 714, 734, 189, 786, 818, 842, 217, 890, 914, 237, 974, 1006, 1038, 532, 1098, 564, 289, 297
OFFSET
2,1
COMMENTS
The sum converges rapidly to 1/2; S(100) = 0.4992..., S(500) = 0.49995....
EXAMPLE
First few fractions are 1/4, 1/3, 7/18, 9/22, 13/30, 4/9, 5/11, 23/50, 7/15, ...
PROG
(PARI) a(n) = denominator(sum(k=2, n, (primepi((k+1)^2) - primepi(k^2))/(primepi((k+1)^2)*primepi(k^2)))); \\ Michel Marcus, Aug 15 2022
CROSSREFS
Cf. A000720 (pi), A161892 (numerators).
Cf. A161621.
Sequence in context: A302851 A373946 A276083 * A192773 A183231 A241358
KEYWORD
nonn,frac,less
AUTHOR
Daniel Tisdale, Jun 21 2009
EXTENSIONS
Offset 2 and more terms from Michel Marcus, Aug 15 2022
STATUS
approved