OFFSET
1,2
COMMENTS
FORMULA
a(n) = 2^A048947(n). - R. J. Mathar, Sep 22 2016
MATHEMATICA
{ta=Table[0, {100}], u=1}; Do[If[IntegerQ[Log[2, DivisorSigma[1, n]]], Print[n]; ta[[u]]=n; u=u+1], {n, 1, 100000000}] DivisorSigma[1, ta]
PROG
(PARI) isok(n) = (n==1) || (ispower(sigma(n), , &r) && (r==2));
for(n=1, 1e7, if(isok(n), print1(sigma(n)", "))) \\ Altug Alkan, Nov 01 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Jun 02 2004
STATUS
approved