login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A137689 Indices m such that A128646(m)-1 is prime, where A128646 = denominator of partial sums of 1/(p(i)-1). 3
3, 4, 5, 7, 8, 9, 10, 11, 15, 16, 23, 24, 26, 47, 48, 54, 78, 79, 80, 243, 244, 245, 246, 247, 367, 368, 369, 370, 371, 372, 373, 374, 375, 447, 453, 635, 636, 1656, 1657, 1658, 1659, 1660 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Terms corresponding to indices m = a(k) > 1000 are not certified primes but at least probable primes. Is there a simple explanation for the large gaps between a(k)=80, a(k+1)=243 and a(k)=636, a(k+1)=1656?
LINKS
PROG
(PARI) print_A137689(i=0/*start checking at i+1*/)={local(s=sum(j=1, i, 1/(prime(j)-1))); while(1, while(!ispseudoprime(-1+denominator(s+=1/(prime(i++)-1))), ); print1(i", "))
CROSSREFS
Sequence in context: A278637 A191979 A072225 * A081690 A081688 A118170
KEYWORD
hard,more,nonn
AUTHOR
M. F. Hasler, Feb 07 2008
EXTENSIONS
Edited by T. D. Noe, Oct 30 2008
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 09:59 EDT 2024. Contains 371268 sequences. (Running on oeis4.)