A343982 - OEIS

%I #14 May 08 2021 01:53:47

%S 1,6,54,135,486,495,516,1134,1863,2295,3375,4374,4875,5535,10935,

%T 11875,15435,19695,22295,23625,24057,34853,39015,39366,42875,43875,

%U 59265,64881,77625,84375,89667,100875,102375,114582,122625,142155,144495,161325,165375,205979,251505,268569

%N Numbers k that divide Sum_{j|k} j^(k/j).

%H Seiichi Manyama, <a href="/A343982/b343982.txt">Table of n, a(n) for n = 1..200</a>

%e 1^6 + 2^3 + 3^2 + 6^1 = 24 = 4 * 6. So 6 is a term.

%t q[n_] := Divisible[DivisorSum[n, #^(n/#) &], n]; Select[Range[10^5], q] (* _Amiram Eldar_, May 06 2021 *)

%o (PARI) isok(n) = sumdiv(n, d, Mod(d, n)^(n/d))==0;

%Y Cf. A007691, A055225, A067313, A336892.

%K nonn

%O 1,2

%A _Seiichi Manyama_, May 06 2021