login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A065523
Numbers n such that sigma(n) is a prime power (A025475).
1
3, 7, 21, 31, 81, 93, 127, 217, 381, 400, 651, 889, 2667, 3937, 8191, 11811, 24573, 27559, 57337, 82677, 131071, 172011, 253921, 393213, 524287, 761763, 917497, 1040257, 1572861, 1777447, 2752491, 3120771, 3670009, 4063201, 5332341
OFFSET
1,1
LINKS
MATHEMATICA
s[n_] := DivisorSigma[1, n ]; Select[ Range[2, 10^6], !PrimeQ[ s[ # ]] && Mod[ s[ # ], s[ # ] - EulerPhi[ s[ # ]]] == 0 & ]
CROSSREFS
A proper subset of A065496.
Sequence in context: A003585 A108102 A357898 * A046528 A018572 A018641
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Nov 27 2001
EXTENSIONS
a(31)-a(35) from Donovan Johnson, Sep 05 2008
STATUS
approved