OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) = 4 because 4 = 2*2 and (2+1)/(2-1) = 3 is prime,
a(2) = 9 because 9 = 3*3 and (3+1)/(3-1) = 2 is prime,
a(3) = 15 because 15 = 5*3 and (5+1)/(3-1) = 3 is prime,
a(4) = 35 because 35 = 7*5 and (7+1)/(5-1) = 2 is prime,
a(5) = 39 because 39 = 13*3 and (13+1)/(3-1) = 7 is prime.
MAPLE
N:= 2000: # to get all terms <= N
P:= select(isprime, [2, seq(i, i=3..N/3, 2)]):
Res:= NULL:
for p in P do qmax:= min(p, N/p);
for q in P do
if q > qmax then break fi;
v:= (p+1)/(q-1);
if v::integer and isprime(v) then Res:= Res, q*p fiod od:
sort([Res]); # Robert Israel, Jul 28 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, May 26 2010
EXTENSIONS
More terms from R. J. Mathar, May 28 2010
STATUS
approved