OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..250
EXAMPLE
The value of 16 in this sequence corresponds to 1+2+3+4+6 = 16 with 1, 2, 3, 4 and 6 being the divisors of 36 <= sqrt(36).
MAPLE
A066839 := proc(n) a := 0 ; for k in numtheory[divisors](n) do if k^2 <= n then a := a+k ; fi; od: a ; end: A143837 := proc() rec := -1; for n from 1 do r := A066839(n) ; if r > rec then printf("%d, ", r) ; rec := r; fi; od: end: A143837() ; # R. J. Mathar, Nov 03 2008
PROG
(PARI) lista(nn) = {my(ms = 0); for (n=1, nn, sqn = sqrt(n); s = sumdiv(n, d, d*(d<=sqn)); if (s > ms, print1(s, ", "); ms = s); ); } \\ Michel Marcus, Oct 05 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
J. Lowell, Sep 02 2008
EXTENSIONS
More terms from R. J. Mathar, Nov 03 2008
Name edited by Jon E. Schoenfield at the suggestion of Joerg Arndt, Oct 12 2015
STATUS
approved