login

%I #12 Aug 21 2022 11:33:52

%S 1,4,6,12,11,23,17,32,29,41,30,66,38,61,61,82,53,104,61,115,92,107,77,

%T 170,98,132,124,170,104,216,114,201,158,183,158,287,143,210,193,293,

%U 161,318,171,291,266,266,189,418,218,335,269,357,220,426,271,429,309,354,250

%N A051731 * A006218.

%F Inverse Mobius transform (A051731) of A006218. Row sums of triangle A143355.

%F a(n) = Sum_{i=1..n} tau(i)*A135539(n,i). - _Ridouane Oudra_, Jul 26 2022

%F a(n) = Sum_{d|n} A006218(d). - _Ridouane Oudra_, Jul 27 2022

%e a(4) = 12 = sum of row 4 terms of triangle A143355: (7, + 3 + 1 + 1).

%e a(4) = 12 = (1, 1, 0, 1) dot (1, 3, 5, 8) = (1 + 3 + 0 + 8), where (1, 1, 0, 1) = row 4 of A051731 and A006218 = (1, 3, 5, 8, 10, 14,...).

%o (PARI) row(n) = my(d=divisors(n)); vector(n, k, #select(x->(x>=k), d)); \\ A135539

%o a(n) = my(v=row(n)); sum(i=1, n, numdiv(i)*v[i]); \\ _Michel Marcus_, Jul 26 2022

%Y Cf. A006218, A051731, A143355.

%Y Cf. A135539.

%K nonn

%O 1,2

%A _Gary W. Adamson_, Aug 10 2008

%E Corrected typo in A-number in formula; added more terms - _R. J. Mathar_, Jan 19 2009