A351946 a(n) = A051903(A181819(n)).

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, 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, 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

Offset: 1

Author

Antti Karttunen, Apr 03 2022

Keywords

• nonn

Links

• Michael De Vlieger, Table of n, a(n) for n = 1..10000
• Index entries for sequences computed from exponents in factorization of n

Crossrefs

Cf. A051903, A181819.

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

Programs

Mathematica

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

PARI

A051903(n) = if((1==n),0,vecmax(factor(n)[, 2]));
A181819(n) = factorback(apply(e->prime(e),(factor(n)[,2])));
A351946(n) = A051903(A181819(n));

Formula

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