A078316 - OEIS

%I #14 Sep 08 2024 22:33:43

%S 1,1,1,2,1,1,2,1,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1,

%T 1,1,1,2,1,1,2,1,2,1,2,1,1,2,3,1,1,1,1,2,1,1,3,1,1,1,1,1,2,1,1,1,1,2,

%U 1,2,1,1,1,1,1,3,1,1,1,2,2,2,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,2,1,1,1,1,1,1

%N Maximum exponent in the prime factorization of n*rad(n)+1, where rad = A007947 (squarefree kernel).

%H Reinhard Zumkeller, <a href="/A078316/b078316.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = A051903(A078310(n)).

%t a[n_] := Max[FactorInteger[1 + n * Times @@ FactorInteger[n][[;;, 1]]][[;;, 2]]]; Array[a, 100] (* _Amiram Eldar_, Sep 07 2024 *)

%o (Haskell)

%o a078316 = a051903 . a078310  -- _Reinhard Zumkeller_, Jul 23 2013

%o (PARI) a(n) = vecmax(factor(1 + n * vecprod(factorint(n)[, 1]))[, 2]); \\ _Amiram Eldar_, Sep 07 2024

%Y Cf. A007947, A051903, A078310, A078313, A078315.

%K nonn

%O 1,4

%A _Reinhard Zumkeller_, Nov 23 2002