OFFSET
1,4
COMMENTS
First differs from A088388 at n = 40.
EXAMPLE
For n = 12, gpf(n)! = 3! = 6 is not divisible by 12, but (3!)^2 = 36 is divisible by 12, so a(12) = 2.
PROG
(Python)
from sympy import factorint
def A362333(n):
f = factorint(n)
gpf = max(f, default=None)
a = 0
for p in f:
m = gpf
v = 0
while m >= p:
m //= p
v += m
a = max(a, -(-f[p]//v))
return a
CROSSREFS
KEYWORD
nonn
AUTHOR
Pontus von Brömssen, Apr 16 2023
STATUS
approved