%I #15 Nov 19 2024 03:16:24
%S 1,2,6,38,551,589,750,2510,3277,7801,9940,11359,16810,19875,21250,
%T 23700,30450,33000,34350,41005,117670,132300,144300,145150,161700,
%U 234635,242550,249000,308700,365750,407925,423801,437250,468150,524700,577850,585206,616585
%N Numbers k such that n | sigma_10(k) + phi(k)^10.
%C sigma_10(k) is the sum of the 10th powers of the divisors of k (A013958).
%H Amiram Eldar, <a href="/A055704/b055704.txt">Table of n, a(n) for n = 1..170</a>
%t Do[If[Mod[DivisorSigma[10, n]+EulerPhi[n]^10, n]==0, Print[n]], {n, 1, 10^5}]
%o (PARI) isok(n) = !((sigma(n, 10) + eulerphi(n)^10) % n); \\ _Michel Marcus_, Mar 02 2014
%Y Cf. A000010, A013958.
%K nonn,changed
%O 1,2
%A _Robert G. Wilson v_, Jun 09 2000
%E More terms from _Michel Marcus_, Mar 02 2014