OFFSET
1,1
MATHEMATICA
lst = Table[n^2 - EulerPhi[n], {n, 100000}]; Select[lst, PrimeQ[#] && PrimeQ[ ( # - 1)/DivisorSigma[0, # - 1] ] &]
PROG
(PARI) lista(nn) = {for (n=1, nn, if (isprime(p=n^2-eulerphi(n)) && (pp=p-1) && (type(r=pp/numdiv(pp))=="t_INT") && isprime(r), print1(p, ", ")); ); } \\ Michel Marcus, Jul 08 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Carlos Eduardo Olivieri, May 30 2015
STATUS
approved