OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..2000 from Robert Israel)
EXAMPLE
30 is a term, because 30^3 = 27000, sigma(30) = 72 and 27000 / 72 = 375.
MAPLE
select(t -> t^3 mod numtheory:-sigma(t) = 0, [$1..20000]); # Robert Israel, Sep 16 2022
MATHEMATICA
Select[Range[16000], Divisible[#^3, DivisorSigma[1, #]] &]
Select[Range[16000], PowerMod[#, 3, DivisorSigma[1, #]]==0&] (* Harvey P. Dale, Sep 04 2024 *)
PROG
(PARI) for(k=1, 10^6, if(k^3%sigma(k)==0, print1(k, ", "))) \\ Alexandru Petrescu, Sep 10 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Zdenek Cervenka, Aug 05 2022
STATUS
approved