OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
ff[x_] := Flatten[FactorInteger[x]] f1[x_] := Length[FactorInteger[x]] f2[x_] := Apply[Plus, Table[Part[ff[x], 2*w], {w, 1, f1[x]}]] Do[s=f2[n]/f1[n]; If[IntegerQ[s]&&Greater[s, 1]&&Greater[f1[n], 1], Print[n]], {n, 2, 10000}]
PROG
(PARI) is(n)=my(f=factor(n)[, 2]); #f>1 && vecsum(f)%#f==0 && vecmax(f)>1 \\ Charles R Greathouse IV, Oct 15 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Jan 16 2002
STATUS
approved