OFFSET
1,2
COMMENTS
See A118370 for the corresponding primes, 'divisorial primes'.
LINKS
Paolo P. Lava, Table of n, a(n) for n = 1..10000
EXAMPLE
The (positive) divisors of 6 are 1,2,3,6. As 6*3*2*1 + 1 = 37 is prime, 6 is a term.
MATHEMATICA
f[n_]:=PrimeQ[Times@@Divisors[n]+1]; lst={}; Do[If[f[n], AppendTo[lst, n]], {n, 7!}]; lst (* Vladimir Joseph Stephan Orlovsky, Dec 14 2009 *)
PROG
(PARI) for(n=1, 2500, s=1; fordiv(n, d, s=s*d); if(isprime(s+1), print1(n, ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Rick L. Shepherd, Apr 25 2006
STATUS
approved