|
|
A060872
|
|
Sum of d*d' over all unordered pairs (d,d') with d*d' = n.
|
|
11
|
|
|
1, 2, 3, 8, 5, 12, 7, 16, 18, 20, 11, 36, 13, 28, 30, 48, 17, 54, 19, 60, 42, 44, 23, 96, 50, 52, 54, 84, 29, 120, 31, 96, 66, 68, 70, 180, 37, 76, 78, 160, 41, 168, 43, 132, 135, 92, 47, 240, 98, 150, 102, 156, 53, 216, 110, 224, 114, 116, 59, 360, 61, 124, 189, 256
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
a(n) is also the sum of all parts of all partitions of n into consecutive parts that differ by 2. - Omar E. Pol, May 05 2020
|
|
LINKS
|
|
|
FORMULA
|
a(n) = n * ceiling( d(n)/2) where d is the number of divisors function.
G.f.: x*f'(x), where f(x) = Sum_{k>=1} x^k^2/(1 - x^k). - Ilya Gutkovskiy, Apr 10 2017
|
|
EXAMPLE
|
a(4)=8 because pairs of factors are 1*4 and 2*2 and 1*4 + 2*2 = 8.
For n = 16 there are three partitions of 16 into consecutive parts that differ by 2 (including 16 as a valid partition). They are [16], [9, 7] and [7, 5, 3, 1]. The sum of the parts is [16] + [9 + 7] + [7 + 5 + 3 + 1] = 48, so a(16) = 48. - Omar E. Pol, May 05 2020
|
|
MATHEMATICA
|
Table[ n * Ceiling[ DivisorSigma[0, n] /2 ], {n, 1, 73} ]
|
|
PROG
|
(Magma) [n*Ceiling(DivisorSigma(0, n)/2): n in [1..70]]; // Vincenzo Librandi, Apr 12 2017
(Python)
from sympy import divisor_count
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|