OFFSET
1,1
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
EXAMPLE
n=33=3*11, mean=(3+11)/2=6.
MATHEMATICA
ffi[x_] := Flatten[FactorInteger[x]] lf[x_] := Length[FactorInteger[x]] ba[x_] := Table[Part[ffi[x], 2*w-1], {w, 1, lf[x]}] ep[x_] := Table[Part[ffi[x], 2*w], {w, 1, lf[x]}] Do[s=Apply[Plus, ba[n]]/lf[n]; If[IntegerQ[s]&&Greater[lf[n], 1], Print[n]], {n, 2, 1000}]
PROG
(Haskell)
a070005 n = a070005_list !! (n-1)
a070005_list = filter ((== 0) . a010055) a078174_list
-- Reinhard Zumkeller, Jun 01 2013
(PARI) lista(nn) = {for (n=2, nn, f = factor(n); if ((#f~ != 1) && (sum(k=1, #f~, f[k, 1]) % #f~ == 0), print1(n, ", ")); ); } \\ Michel Marcus, Mar 28 2015
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Labos Elemer, Apr 11 2002
STATUS
approved