|
|
A034913
|
|
Odd primes p such that q = (k*p+1)/(p-k) is prime for some p/2 < k < p.
|
|
2
|
|
|
3, 7, 13, 17, 23, 31, 47, 53, 67, 73, 79, 113, 137, 139, 151, 157, 163, 173, 193, 227, 257, 293, 307, 317, 337, 349, 353, 379, 401, 419, 457, 463, 467, 479, 487, 499, 541, 557, 577, 593, 599, 613, 617, 643, 653, 677, 691, 727, 733, 769, 787, 811, 829, 853
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Related to hyperperfect numbers of a certain form.
Old name was: Primes such that q=(k*a(n)+1)/(a(n)-k), is prime for some k, q>a(n).
|
|
LINKS
|
|
|
EXAMPLE
|
137 is a term since (132*137+1)/(137-132) = 3617 is prime.
|
|
PROG
|
(PARI) isok(p)=for (k=1, p-1, my(q = (k*p+1)/(p-k)); if ((q > p) && (denominator(q)==1) && isprime(q), return (1)); );
lista(nn) = forprime(p=3, nn, if (isok(p), print1(p, ", "))) \\ Michel Marcus, Mar 11 2016
(PARI) is(p)=fordiv(p^2+1, d, if(d>=p/2, return(0)); if(isprime((p*(p-d)+1)/d), return(isprime(p)))) \\ Charles R Greathouse IV, Mar 11 2016
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|