OFFSET
1,1
LINKS
Paolo P. Lava, Table of n, a(n) for n = 1..1000
MAPLE
with(numtheory); P:=proc(q) local a, n, p; for n from 1 to q do
a:=n*add(op(2, p)/op(1, p), p=ifactors(n)[2]);
if isprime(n+a) and isprime(n-a) then print(n); fi;
od; end: P(10^5);
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Sep 18 2013
STATUS
approved