OFFSET
2,3
LINKS
Pontus von Brömssen, Table of n, a(n) for n = 2..10000
EXAMPLE
For n = 80 = 2^4 * 3^0 * 5^1, gpf(80)! = 5! = 2^3 * 3^1 * 5^1. The ratios of the prime exponents are 4/3, 0/1, and 1/1, the greatest of which is 4/3, so a(80) = 4.
PROG
(Python)
from sympy import factorint, Rational
def A371148(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, Rational(f[p], v))
return a.p
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Pontus von Brömssen, Mar 13 2024
STATUS
approved