OFFSET
1,1
PROG
(PARI) f(n)= my(s=0); for(i=1, n, s=s+1/(prime(i)-i)); return(s);
for (i=1, 500, if(isprime(i) && isprime(numerator(f(i))), print1(i, ", ")));
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Mohammed Bouayoun (mohammed.bouayoun(AT)sanef.com), Feb 20 2004
EXTENSIONS
a(9)-a(12) from Amiram Eldar, Jul 22 2025
STATUS
approved
