login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A359503 Partial sums of A066839. 1
1, 2, 3, 6, 7, 10, 11, 14, 18, 21, 22, 28, 29, 32, 36, 43, 44, 50, 51, 58, 62, 65, 66, 76, 82, 85, 89, 96, 97, 108, 109, 116, 120, 123, 129, 145, 146, 149, 153, 165, 166, 178, 179, 186, 195, 198, 199, 215, 223, 231, 235, 242, 243, 255, 261, 275, 279, 282, 283 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
a(n) is the sum of all divisors d of k such that d^2 <= k where k ranges from 1 to n.
LINKS
FORMULA
a(n) = m*(6*n+5-m*(2*m+3))/6 + Sum_{k=1..n, i=1..floor(sqrt(k))} [(k-1) mod i] - [k mod i] where m = floor(sqrt(n)).
a(n) = m*(6*n+5-m*(2*m+3))/6 + Sum_{k=1..n, i=1..floor(sqrt(k))} (k-1) mod i - Sum_{k=1..n} A176314(k) where m = floor(sqrt(n)).
MATHEMATICA
Table[Select[Divisors[n], # <= Sqrt[n]&]//Total, {n, 1, 60}]//Accumulate (* Jean-François Alcover, Jan 26 2024 *)
PROG
(Python)
from itertools import takewhile
from sympy import divisors
def A359503(n): return sum(sum(takewhile(lambda x:x**2<=i, divisors(i))) for i in range(1, n+1))
CROSSREFS
Sequence in context: A062837 A190670 A226115 * A073170 A014689 A117206
KEYWORD
nonn
AUTHOR
Chai Wah Wu, Jan 24 2024
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 2 08:00 EST 2024. Contains 370461 sequences. (Running on oeis4.)