A351946 - OEIS

%I #9 Apr 03 2022 13:53:13

%S 0,1,1,1,1,2,1,1,1,2,1,1,1,2,2,1,1,1,1,1,2,2,1,1,1,2,1,1,1,3,1,1,2,2,

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

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

%N a(n) = A051903(A181819(n)).

%H Michael De Vlieger, <a href="/A351946/b351946.txt">Table of n, a(n) for n = 1..10000</a>

%H <a href="/index/Eu#epf">Index entries for sequences computed from exponents in factorization of n</a>

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

%t {0}~Join~Array[Max[FactorInteger[#][[All, -1]]] &@ Apply[Times, Prime[FactorInteger[#][[All, -1]]]] &, 104, 2] (* _Michael De Vlieger_, Apr 03 2022 *)

%o (PARI)

%o A051903(n) = if((1==n),0,vecmax(factor(n)[, 2]));

%o A181819(n) = factorback(apply(e->prime(e),(factor(n)[,2])));

%o A351946(n) = A051903(A181819(n));

%Y Cf. A051903, A181819.

%Y Differs from A327500 for the first time at n=450, where a(450) = 2, while A327500(450) = 3. See A351947 for all such positions.

%K nonn

%O 1,6

%A _Antti Karttunen_, Apr 03 2022