login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A067313 Numbers k such that sigma_k(k)/k is an integer, where sigma_k(k) is the sum of the k-th powers of the divisors of k (A023887). 2

%I #20 Aug 10 2020 13:19:14

%S 1,10,130,135,147,150,228,250,350,364,410,444,492,876,891,945,1014,

%T 1308,1372,1550,1690,1950,2050,2210,2373,2565,2850,3045,3050,3250,

%U 3375,3876,4108,4185,4905,4995,5050,5070,5145,5330,5439,5481,6150,6250,6321,6615,6890,7514

%N Numbers k such that sigma_k(k)/k is an integer, where sigma_k(k) is the sum of the k-th powers of the divisors of k (A023887).

%H David A. Corneth, <a href="/A067313/b067313.txt">Table of n, a(n) for n = 1..10000</a> (first 1000 terms from Amiram Eldar)

%e If k = 10, then sigma(10,10) = 1 + 1024 + 9765625 + 10000000000 = 1009766650 is divisible by k = 10.

%t Do[s=DivisorSigma[n, n]; If[IntegerQ[s/n], Print[n]], {n, 1, 10000}]

%t Select[Range[6500],Divisible[DivisorSigma[#,#],#]&] (* _Harvey P. Dale_, Feb 06 2019 *)

%o (PARI) isok(k) = !(sigma(k, k) % k); \\ _Michel Marcus_, Aug 10 2020

%o (PARI) is(n) = {my(d = divisors(n), v = vecsum(vector(#d - 1, i, Mod(d[i], n)^n))); lift(v)==0} \\ _David A. Corneth_, Aug 10 2020

%Y Cf. A023887.

%K nonn

%O 1,2

%A _Labos Elemer_, Jan 14 2002

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 05:43 EDT 2024. Contains 371264 sequences. (Running on oeis4.)