login
Numbers k such that sopfr((k-1)!) is divisible by k, where sopfr(k) = A001414(k) = sum of primes, with repetition, dividing k.
1

%I #28 May 17 2021 01:52:32

%S 1,2,45,53,177,436,1239,3651,6463,6869,10753,19450,29721,33289,88907,

%T 93682,1137232,1516121,4361271,9428534,43778664,74738670,271442366,

%U 775223371,835126289,1736463189,3088442241,5054888590,11184483614,16993011938,30788570768,33342871740

%N Numbers k such that sopfr((k-1)!) is divisible by k, where sopfr(k) = A001414(k) = sum of primes, with repetition, dividing k.

%C See A025281(k-1) for the values of sopfr((k-1)!).

%H Martin Ehrenstein, <a href="/A343424/b343424.txt">Table of n, a(n) for n = 1..35</a>

%e 45 is a term as sopfr(44!) = 585 which is divisible by 45.

%t sopfr[0] = sopfr[1] = 0; sopfr[n_] := Plus @@ Times @@@ FactorInteger[n]; sum = 0; s = {}; Do[sum += sopfr[n]; If[Divisible[sum, n + 1], AppendTo[s, n + 1]], {n, 0, 10^6}]; s (* _Amiram Eldar_, May 06 2021 *)

%o (PARI) sopfr(n) = (n=factor(n))[, 1]~*n[, 2]; \\ A001414

%o isok(k) = !(sopfr((k-1)!) % k); \\ _Michel Marcus_, May 06 2021

%Y Cf. A001414, A025281, A027746.

%K nonn

%O 1,2

%A _Scott R. Shannon_, Apr 15 2021

%E a(27) from _Amiram Eldar_, May 06 2021

%E a(28) and beyond from _Martin Ehrenstein_, May 16 2021