OFFSET
1,1
COMMENTS
All Mersenne primes are in the sequence. Because if p=2^q-1 is prime then 2*phi(p+1)=2*phi(2^q)=2^q=p+1=sigma(p). There are no other prime terms. - Farideh Firoozbakht, Aug 14 2014
No terms beyond a(80) up to n = 5 million. - Harvey P. Dale, Mar 19 2016
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..758 (terms below 10^10, terms 1..80 from Harvey P. Dale)
MAPLE
select(t -> numtheory:-sigma(t) = 2*numtheory:-phi(t+1), [$1..10^6]); # Robert Israel, Aug 14 2014
MATHEMATICA
With[{nn=250000}, Position[Thread[{DivisorSigma[1, Range[nn]], 2*EulerPhi[ Range[ 2, nn+1]]}], {x_, x_}]]//Flatten (* Harvey P. Dale, Mar 19 2016 *)
PROG
(PARI) isok(n) = sigma(n) == 2*eulerphi(n+1); \\ Michel Marcus, Nov 24 2013
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Mar 09 2002
STATUS
approved