OFFSET
0,5
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000
MAPLE
with(numtheory):
d:= n-> n*add(i[2]/i[1], i=ifactors(n)[2]):
a:= n-> (f-> `if`(f<2, 0, min(factorset(f)[])))(d(n)):
seq(a(n), n=0..105); # Alois P. Heinz, Jun 08 2015
MATHEMATICA
d[n_] := n*Sum[i[[2]]/i[[1]], {i, FactorInteger[n]}];
a[n_] := Function[f, If[f<2, 0, Min[FactorInteger[f][[All, 1]]]]][d[n]]; a[0] = 0;
Table[a[n], {n, 0, 105}] (* Jean-François Alcover, Mar 23 2017, after Alois P. Heinz *)
PROG
(Python)
from sympy import primefactors, factorint
def A086134(n): return 0 if n <= 1 else min(primefactors(m)) if (m:=sum((n*e//p for p, e in factorint(n).items()))) > 1 else 0 # Chai Wah Wu, Nov 04 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Jul 23 2003
STATUS
approved