A381133 If n = (p_1^e_1)*(p_2^e_2)*(p_3^e_3)*... and min(p_1^e_1,p_2^e_2,...) = p_k^e_k then a(n) = e_k, a(1) = 0.

0, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 3, 2, 1, 1, 1, 5, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 6, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 2

Offset

1,4

Links

Table of n, a(n) for n=1..100.

Mathematica

Table[PrimeOmega[Min @@ (#[[1]]^#[[2]] & /@ FactorInteger[n])], {n, 100}]

Prog

(PARI) a(n) = if (n==1, 0, my(f=factor(n), v=vector(#f~, k, f[k, 1]^f[k, 2]), m=vecmin(v), i=select(x->(x==m), v, 1)); f[i[1], 2]); \\ Michel Marcus, Feb 19 2025

Crossrefs

Cf. A034684, A051904, A088388, A381131, A381132.

Sequence in context: A070012 A375034 A373365 * A071178 A366895 A326515

Adjacent sequences: A381130 A381131 A381132 * A381134 A381135 A381136

Keyword

nonn

Author

Ilya Gutkovskiy, Feb 14 2025

Status

approved