|
|
A049784
|
|
a(n) = Sum_{k=1..n} T(n,k), array T as in A049783.
|
|
7
|
|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,5
|
|
LINKS
|
|
|
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)) ));
(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
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|