OFFSET
1,2
COMMENTS
It is obvious that if n=2*p where p is a Sophie Germain odd prime then n=phi(phi(n)+sigma(n)). This sequence is a subsequence of A097646. Except for the first term all terms of this sequence are even. There is no other term up to 3*10^7.
LINKS
C. K. Caldwell, The Prime Glossary, Sophie Germain prime.
EXAMPLE
14592000 is in the sequence because 14592000=2*7296000, 7296000 is not a Sophie Germain odd prime and phi(phi(14592000)+sigma(14592000)) =14592000.
MATHEMATICA
Do[If[(!PrimeQ[n/2]||!PrimeQ[n+1])&&n==EulerPhi[EulerPhi[n]+ DivisorSigma[1, n]], Print[n]], {n, 30000000}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Farideh Firoozbakht, Sep 09 2004
EXTENSIONS
a(24)-a(33) from Donovan Johnson, Jan 18 2012
STATUS
approved