OFFSET
1,2
COMMENTS
EXAMPLE
6 is a term because 6^sigma(6) / sigma(6)^6 = 6^12 / 12^6 = 2176782336 / 2985984 = 729 (integer).
MAPLE
with(numtheory):
select(n->n &^ sigma(n) mod sigma(n)^n =0, [`$`(1..30000)]); # Muniru A Asiru, Mar 20 2018
PROG
(Magma) [n: n in[1..20000] | n^SumOfDivisors(n) mod SumOfDivisors(n)^n eq 0]
(GAP) Filtered([1..30000], n->PowerModInt(n, Sigma(n), Sigma(n)^n)=0); # Muniru A Asiru, Mar 20 2018
(PARI) isok(n) = my(s = sigma(n)); Mod(n, s^n)^s == 0; \\ Michel Marcus, Mar 23 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Mar 20 2018
STATUS
approved