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
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
J. S. McCranie, A study of hyperperfect numbers, J. Int. Seqs. Vol. 3 (2000) #P00.1.3.
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
Name edited by Michel Marcus and Charles R Greathouse IV, Mar 11 2016
STATUS
approved