OFFSET
0,2
LINKS
Peter Luschny, Table of n, a(n) for n = 0..10000
FORMULA
MATHEMATICA
A363734[n_]:=If[n==0, 0, n+1+(n-1)DivisorSigma[0, n]]; Array[A363734, 100, 0] (* Paolo Xausa, Aug 06 2023 *)
PROG
(SageMath)
print([sum(n^k.divides(n) for k in srange(n+1)) for n in srange(57)])
(Python)
from sympy import divisor_count
def A363734(n): return (n-1)*divisor_count(n)+n+1 if n else 0 # Chai Wah Wu, Jun 28 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Jun 27 2023
STATUS
approved