OFFSET
1,6
LINKS
T. D. Noe, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = n/A053585(n).
EXAMPLE
a(36) = 4 because 36/3^2 = 4, 3^2 is highest power dividing 36 of largest prime dividing 36.
a(50) = 50/5^2 = 2.
MATHEMATICA
f[n_]:=Module[{c=Last[FactorInteger[n]]}, n/First[c]^Last[c]]; Array[ f, 110] (* Harvey P. Dale, Oct 14 2011 *)
PROG
(Python)
from sympy import factorint, primefactors
def a053585(n):
if n==1: return 1
p = primefactors(n)[-1]
return p**factorint(n)[p]
def a(n): return n/a053585(n) # Indranil Ghosh, May 19 2017
(PARI) a(n) = if(n>1, my(f=factor(n)); n/f[#f~, 1]^f[#f~, 2], 1); \\ Michel Marcus, Jan 10 2025
CROSSREFS
KEYWORD
nonn,easy,nice
AUTHOR
EXTENSIONS
More terms from James A. Sellers, Jan 21 2000
STATUS
approved