|
|
A178379
|
|
Semiprimes of the form m*k such that (m+1)/(k-1) is prime.
|
|
3
|
|
|
4, 9, 15, 35, 39, 55, 77, 95, 111, 119, 183, 203, 209, 215, 219, 287, 299, 319, 335, 471, 527, 579, 707, 767, 791, 799, 815, 831, 939, 959, 989, 1007, 1055, 1079, 1191, 1199, 1211, 1263, 1343, 1371, 1415, 1623, 1655, 1691, 1703, 1799, 1829, 1839, 1967, 1983
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
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:
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|