OFFSET
4,1
COMMENTS
a(n) >= A007917(n); a(n) is a power of some prime p <= n.
PROG
(PARI) works(n, m) = local(f, s, l, p, x); f = factor(n!); s = 0; l = matsize(f)[1]; if (f[l, 1] > m, return(0)); for (i = 1, l, p = f[i, 1]; x = 1; while (p^x <= m, x++); x--; s += f[i, 2]\x; if (f[i, 2]%x, s++)); s <= n; a(n) = local(f, m); f = factor(n); m = f[matsize(f)[1], 1]; while (!works(n, m), m++); m
CROSSREFS
KEYWORD
nonn,less
AUTHOR
David Wasserman, Jan 31 2005
STATUS
approved