login
Sum of primes dividing n^n+1 (with repetition).
1

%I #14 Jun 21 2018 02:01:57

%S 2,5,11,257,526,147,1030,1027,530793,31603,58685,2228292,113060,

%T 180326,163123,67280421584898,45957792327018709129,33414185,

%U 870543318650,4406613081350403,22864393425065,82812579069940,1576297793,27266016518

%N Sum of primes dividing n^n+1 (with repetition).

%H Harry J. Smith, <a href="/A064772/b064772.txt">Table of n, a(n) for n = 1..61</a>

%t Table[Total[Times@@@FactorInteger[n^n+1]],{n,25}] (* _Harvey P. Dale_, Sep 17 2011 *)

%o (PARI) sopfr(n)= { local(f,s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]*f[i, 2]); return(s) } { for (n=1, 61, write("b064772.txt", n, " ", sopfr(n^n + 1)) ) } \\ _Harry J. Smith_, Sep 24 2009

%Y Cf. A014566, A001414.

%K nonn

%O 1,1

%A _Jason Earls_, Oct 19 2001