OFFSET
1,1
COMMENTS
a(n) = gpf(n)^2 for n < gpf(n)^2; and for sqrt(n)-smooth numbers: a(n) <= gpf(n)*(least prime > n/gpf(n)).
MATHEMATICA
mfs={#, First[FactorInteger[#]][[1]]}&/@Range[3500];
Table[Select[Select[mfs, #[[2]]==Last[FactorInteger[n]][[1]]&], #[[1]]>n&][[1, 1]], {n, 2, 60}] (* Harvey P. Dale, Feb 18 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Sep 29 2002
EXTENSIONS
Corrected by Harvey P. Dale, Feb 18 2011
STATUS
approved