OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MAPLE
a:= n-> add(round(sqrt(d)), d=numtheory[divisors](n)):
seq(a(n), n=1..80); # Alois P. Heinz, Mar 02 2020
MATHEMATICA
Table[DivisorSum[n, Floor[1/2+Sqrt[#]]&], {n, 80}]
PROG
(PARI) a(n) = sumdiv(n, d, round(sqrt(d))); \\ Michel Marcus, Mar 03 2020
(Python)
from math import isqrt
from sympy import divisors
def A332931(n): return sum((m:=isqrt(d))+int(d-m*(m+1)>=1) for d in divisors(n, generator=True)) # Chai Wah Wu, Aug 03 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Harvey P. Dale, Mar 02 2020
STATUS
approved