OFFSET
0,9
COMMENTS
Partial sums of A001899.
LINKS
Robert Israel, Table of n, a(n) for n = 0..10000
MAPLE
g:= n -> nops(select(t -> t mod 5 = 4, numtheory:-divisors(n))):
g(0):= 0:
ListTools:-PartialSums(map(g, [$0..100])); # Robert Israel, Apr 29 2021
PROG
(PARI) a(n)=sum(k=0, n, (n\(5*k+4)))
(Maxima) A218447[n]:=sum(floor(n/(5*k+4)), k, 0, n)$
makelist(A218447[n], n, 0, 80); /* Martin Ettl, Oct 20 2012 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Oct 28 2012
STATUS
approved