login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A227427 Numbers n such that sum_{i=1..n} i^sigma(i) == 0 (mod n). 8
1, 3, 17, 64, 223, 816, 3536, 22704, 27549 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

a(10) > 10^5. - Giovanni Resta, Jul 11 2013

LINKS

Table of n, a(n) for n=1..9.

EXAMPLE

1^sigma(1) + 2^sigma(2) + … + 16^sigma(16) + 17^sigma(17) =

1^1 + 2^3 + 3^4 + 4^7 + 5^6 + 6^12 + 7^8 + 8^15 + 9^13 + 10^18 + 11^12 + 12^28 + 13^14 + 14^24 + 15^24 + 16^31 + 17^18 = 21267649601053603536507860737702745369 and 21267649601053603536507860737702745369 / 17 = 1251038211826682560971050631629573257.

MAPLE

with(numtheory); ListA227427:=proc(q) local i, n;

for n from 1 to q do if add(i^sigma(i), i=1..n) mod n=0 then print(n);

fi; od; end: ListA227427(10^6);

MATHEMATICA

With[{nn=40000}, Transpose[Select[Thread[{Accumulate[Table[n^DivisorSigma[ 1, n], {n, nn}]], Range[nn]}], Divisible[#[[1]], #[[2]]]&]][[2]]] (* Harvey P. Dale, Jun 08 2016 *)

CROSSREFS

Cf. A000203, A227429

Sequence in context: A208996 A273949 A174285 * A120386 A061982 A101562

Adjacent sequences:  A227424 A227425 A227426 * A227428 A227429 A227430

KEYWORD

nonn

AUTHOR

Paolo P. Lava, Jul 11 2013

EXTENSIONS

a(7)-a(9) from Giovanni Resta, Jul 11 2013

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified November 19 05:29 EST 2018. Contains 317333 sequences. (Running on oeis4.)