

A049789


a(n) = T(n,n4), array T as in A049783.


7



0, 0, 0, 0, 0, 0, 1, 0, 5, 5, 7, 4, 7, 10, 6, 9, 14, 10, 17, 10, 14, 16, 24, 12, 19, 20, 24, 18, 32, 16, 25, 23, 30, 31, 33, 19, 31, 34, 43, 28, 45, 26, 41, 37, 37, 41, 61, 27, 44, 39, 49, 41, 63, 47, 58, 35, 48, 54, 76, 37, 61, 68, 65, 46, 65
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

5,9


LINKS



FORMULA

a(n) = Sum_{j=1..n8} mod(n4, floor((n8)/j)).  G. C. Greubel, Dec 12 2019


MAPLE

seq( add(`mod`(n4, floor((n8)/j)), j=1..n8), n=5..70); # G. C. Greubel, Dec 12 2019


MATHEMATICA

Table[Sum[Mod[n4, Floor[(n8)/j]], {j, n8}], {n, 5, 70}] (* G. C. Greubel, Dec 12 2019 *)


PROG

(PARI) vector(70, n, sum(j=1, n4, lift(Mod(n, (n4)\j))) ) \\ G. C. Greubel, Dec 12 2019
(Magma) [ n lt 9 select 0 else &+[((n4) mod Floor((n8)/j)): j in [1..n8]]: n in [5..70]]; // G. C. Greubel, Dec 12 2019
(Sage) [sum( (n4)%floor((n8)/j) for j in (1..n8)) for n in (5..70)] # G. C. Greubel, Dec 12 2019
(GAP) List([5..70], n> Sum([1..n8], j> (n4) mod Int((n8)/j)) ); # G. C. Greubel, Dec 12 2019


CROSSREFS



KEYWORD

nonn


AUTHOR



STATUS

approved



