OFFSET
1,1
COMMENTS
From Robert Israel, Feb 26 2017: (Start)
Numbers n such that sigma(n) is in A033845.
Contains A139257. (End)
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= proc(n) local s; s:= numtheory:-sigma(n); s mod 6 = 0 and s = 2^padic:-ordp(s, 2)*3^padic:-ordp(s, 3) end proc:
select(filter, [$1..10000]); # Robert Israel, Feb 26 2017
MATHEMATICA
s3pQ[n_]:=Module[{s=DivisorSigma[1, n]}, s==3*EulerPhi[s]]; Select[ Range[ 400], s3pQ] (* Harvey P. Dale, May 28 2015 *)
PROG
(PARI) for(n=1, 500, if(sigma(n)==3*eulerphi(sigma(n)), print1(n, ", "))) \\ Derek Orr, Feb 26 2017
(PARI) is(n)=my(s=sigma(n)); s%3==0 && s==3*eulerphi(s) \\ Charles R Greathouse IV, Feb 27 2017
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Feb 22 2002
STATUS
approved