OFFSET
1,1
COMMENTS
No additional terms up to 1 million. - Harvey P. Dale, Nov 13 2013
EXAMPLE
sigma(25) * p(25) = 3007 = 31*97.
MATHEMATICA
brilQ[n_]:=Module[{pf=Transpose[FactorInteger[n]][[1]]}, PrimeOmega[n] == 2&&IntegerLength[First[pf]]==IntegerLength[Last[pf]]]; Select[Range[ 100000], brilQ[DivisorSigma[1, #]Prime[#]]&] (* Harvey P. Dale, Nov 13 2013 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Giovanni Resta, Feb 13 2006
STATUS
approved