OFFSET
1,2
COMMENTS
In other words, a(n) is the sum of squares of terms in row n of A235791.
MAPLE
a:= n-> add(ceil((n+1)/k-(k+1)/2)^2, k=1..floor((sqrt(8*n+1)-1)/2)):
seq(a(n), n=1..60); # Alois P. Heinz, Dec 18 2020
PROG
(PARI) a(n) = sum(k=1, floor((sqrt(8*n+1)-1)/2), ceil((n+1)/k-(k+1)/2)^2) \\ Felix Fröhlich, Dec 19 2020; adapted from Maple code
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Dec 18 2020
STATUS
approved