A049784 a(n) = Sum_{k=1..n} T(n,k), array T as in A049783.

0, 0, 1, 1, 6, 3, 14, 12, 21, 21, 51, 21, 62, 69, 79, 63, 136, 90, 181, 129, 160, 186, 324, 160, 264, 306, 371, 297, 512, 285, 528, 444, 547, 628, 768, 372, 715, 820, 980, 652, 1079, 702, 1173, 1053, 1002, 1155, 1728, 960, 1417, 1292

Offset

1,5

Links

G. C. Greubel, Table of n, a(n) for n = 1..1000

Maple

seq( add(add(`mod`(n, floor(k/j)), j=1..k), k=1..n), n=1..60); # G. C. Greubel, Dec 12 2019

Mathematica

Table[Sum[Sum[Mod[n, Floor[k/j]], {j, k}], {k, n}], {n, 60}] (* G. C. Greubel, Dec 12 2019 *)

Prog

(PARI) a(n) = sum(k=1, n, sum(j=1, k, lift(Mod(n, k\j)) ));
vector(60, n, a(n)) \\ G. C. Greubel, Dec 12 2019
(Magma) [ &+[(&+[(n mod Floor(k/j)): j in [1..k]]): k in [1..n]]: n in [1..60]]; // G. C. Greubel, Dec 12 2019
(Sage) [sum(sum( n%floor(k/j) for j in (1..k)) for k in (1..n)) for n in (1..60)] # G. C. Greubel, Dec 12 2019
(GAP) List([1..60], n-> Sum([1..n], k-> Sum([1..k], j-> n mod Int(k/j)) )); # G. C. Greubel, Dec 12 2019

Crossrefs

Row sums of A049783.

Cf. A049785, A049786, A049787, A049788, A049789.

Sequence in context: A214498 A128756 A345056 * A341746 A097917 A116570

Adjacent sequences: A049781 A049782 A049783 * A049785 A049786 A049787

Keyword

nonn

Author

Clark Kimberling

Status

approved