OFFSET
1,1
COMMENTS
Note that the definition here is subtly different from that of A092065.
Terms a(k) < 1000 correspond to primes. Beyond, numerators are probable primes. Note that A120271(3100) has 2187 digits. - M. F. Hasler, Feb 06 2008
Intersection of A000040 (the primes) and A120271 (numerators of partial sums of 1/(prime(i)-1)). - M. F. Hasler, Feb 06 2008
a(60) > 10000. - Jason Yuen, Aug 26 2024
EXAMPLE
1/(2-1) + 1/(3-1) = 3/2 and 3 is prime so a(1)=2.
MATHEMATICA
Position[Accumulate[1/(Prime[Range[3100]]-1)], _?(PrimeQ[ Numerator[ #]]&)]//Flatten (* Harvey P. Dale, Oct 16 2016 *)
PROG
(PARI) A120271(n) = numerator(sum(k=1, n, 1/(prime(k)-1)));
for (i=1, 500, if(isprime(A120271(i)), print1(i, ", ")));
(PARI) print_A092063( i=0 /* start testing at i+1 */)={local(s=sum(j=1, i, 1/(prime(j)-1))); while(1, while(!ispseudoprime(numerator(s+=1/(prime(i++)-1))), ); print1(i", "))} \\ M. F. Hasler, Feb 06 2008
CROSSREFS
KEYWORD
hard,nonn
AUTHOR
Mohammed Bouayoun (mohammed.bouayoun(AT)sanef.com), Feb 20 2004
EXTENSIONS
More terms from M. F. Hasler, Feb 06 2008
Edited by T. D. Noe, Oct 30 2008
Corrected by Harvey P. Dale, Oct 16 2016
a(48)-a(59) from Jason Yuen, Aug 26 2024
STATUS
approved