|
%I #55 Dec 31 2023 11:29:48
%S 1,4,17,19,148,1577,3564,4388,5873,6639,8579,62500,376636,792949,
%T 996044,1174065,3333551,5179004,7516003
%N Numbers k such that k divides Sum_{j=1..k} j^j = A001923(k).
%C a(20) > 10^7. - _Hiroaki Yamanouchi_, Aug 25 2015
%p a:=0:
%p for n from 1 to 2000 do
%p a:=a+n^n:
%p if a mod n=0 then
%p print(n);
%p fi;
%p od: # Revised program from _R. J. Mathar_, Jun 18 2015
%t f=0; Do[ f=f+k^k; If[ IntegerQ[f/k], Print[k] ], {k,1,6639} ]
%o (PARI) for(n=1,10^4, s=sum(i=1,n,Mod(i,n)^i); if(!Mod(s,n), print1(n,", "))) \\ _Derek Orr_, Jun 18 2015
%o (Python)
%o from itertools import accumulate, count, islice
%o def A128981_gen(): # generator of terms
%o yield 1
%o for i, j in enumerate(accumulate(k**k for k in count(1)),start=2):
%o if j % i == 0:
%o yield i
%o A128981_list = list(islice(A128981_gen(),10)) # _Chai Wah Wu_, Jun 18 2022
%Y Cf. A062970, A001923.
%K hard,more,nonn
%O 1,2
%A _Alexander Adamchuk_, Apr 29 2007
%E a(11) and a(12) from _Jon E. Schoenfield_, May 09 2007
%E a(13) = 376636 from _Alexander Adamchuk_, May 03 2010
%E a(14)-a(16) from _Lars Blomberg_, May 10 2011
%E a(17) from _Giovanni Resta_, Jul 13 2015
%E a(18)-a(19) from _Hiroaki Yamanouchi_, Aug 25 2015
|
