OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..380 (first 100 terms from T. D. Noe)
FORMULA
EXAMPLE
MAPLE
a:= n-> n^(n+1)*add(i[2]/i[1], i=ifactors(n)[2]):
seq(a(n), n=0..30); # Alois P. Heinz, Jun 09 2015
# alternative
A068327 := proc(n)
A003415(n^n) ;
end proc:
seq( A068327(n), n=0..10) ; # R. J. Mathar, Oct 19 2021
MATHEMATICA
a312[n_] := Sum[ StirlingS2[n, k]*n!/(n - k)!, {k, 0, n}]; a3415[n_] := With[ {fi = FactorInteger[n]}, n*Total[ fi[[All, 2]] / fi[[All, 1]] ] ]; a3415[0] = a3415[1] = 0; a[n_] := a3415[ a312[n] ]; Table[ a[n], {n, 1, 16}] (* Jean-François Alcover, Mar 27 2013 *)
PROG
(Python)
from sympy import factorint
def A068327(n): return sum((n**(n+1)*e//p for p, e in factorint(n).items())) if n > 1 else 0 # Chai Wah Wu, Jun 12 2022
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
Reinhard Zumkeller, Feb 27 2002
STATUS
approved