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