OFFSET
4,7
LINKS
G. C. Greubel, Table of n, a(n) for n = 4..1000
FORMULA
a(n) = Sum_{j=1..n-6} mod(n-3, floor((n-6)/j)). - G. C. Greubel, Dec 12 2019
MAPLE
seq( add(`mod`(n-3, floor((n-6)/j)), j=1..n-6), n=4..70); # G. C. Greubel, Dec 12 2019
MATHEMATICA
Table[Sum[Mod[n-3, Floor[(n-6)/j]], {j, n-6}], {n, 4, 70}] (* G. C. Greubel, Dec 12 2019 *)
PROG
(PARI) vector(70, n, sum(j=1, n-3, lift(Mod(n, (n-3)\j))) ) \\ G. C. Greubel, Dec 12 2019
(Magma) [n lt 7 select 0 else &+[((n-3) mod Floor((n-6)/j)): j in [1..n-6]]: n in [4..70]]; // G. C. Greubel, Dec 12 2019
(Sage) [sum( (n-3)%floor((n-6)/j) for j in (1..n-6)) for n in (4..70)] # G. C. Greubel, Dec 12 2019
(GAP) List([4..70], n-> Sum([1..n-6], j-> (n-3) mod Int((n-6)/j)) ); # G. C. Greubel, Dec 12 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved