OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Least Prime Factor.
Eric Weisstein's World of Mathematics, Squareful.
FORMULA
MATHEMATICA
Select[ FactorInteger[#], #[[2]]>1&, 1][[1, 1]]& /@ Select[ Range[300], !SquareFreeQ[#]& ] (* Jean-François Alcover, Nov 06 2012 *)
PROG
(Python)
from math import isqrt
from sympy import mobius, factorint
def A046027(n):
def f(x): return n+sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1))
m, k = n, f(n)
while m != k:
m, k = k, f(k)
s = factorint(m)
return next(p for p in sorted(s) if s[p]>1) # Chai Wah Wu, Jul 22 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved