OFFSET
1,2
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
PROG
(Scheme, two different implementations)
(PARI) a(n) = my(f = factor(core(n))); if (!#f~, 1, vecmin(f[, 1])); \\ Michel Marcus, Oct 30 2016
(Python)
from sympy import primefactors
from sympy.ntheory.factor_ import core
def lpf(n): return 1 if n==1 else primefactors(n)[0]
def a(n): return lpf(core(n)) # Indranil Ghosh, May 17 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Oct 28 2016
STATUS
approved