login
Nonprime numbers k that divide A124271(k) = Sum_{i=1..k} (prime(i)^k - 1) / (prime(i) - 1).
2

%I #12 Dec 10 2024 21:22:20

%S 1,9,15,121

%N Nonprime numbers k that divide A124271(k) = Sum_{i=1..k} (prime(i)^k - 1) / (prime(i) - 1).

%C The next term if it exists is greater than 1000.

%C Note that a(1) = 1, a(2) = 3^2 and a(4) = 11^2 are perfect squares.

%e 9 is a term because 9 divides A124271(9) = 108336139407.

%t s={};Do[If[!PrimeQ[k],If[Divisible[Sum[(Prime[i]^k-1)/(Prime[i]-1),{i,k}],k],AppendTo[s,k]]],{k,10^3}];s (* _James C. McMahon_, Dec 10 2024 *)

%Y Cf. A124271, A124273 (primes p that divide A124271(p)).

%K bref,hard,more,nonn

%O 1,2

%A _Alexander Adamchuk_, Oct 23 2006