login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

a(n) = Sum_{i|n, j|n, k|n} i*j*k/lcm(i,j,k).
3

%I #16 May 12 2021 10:16:51

%S 1,14,22,93,44,308,74,472,259,616,158,2046,212,1036,968,2123,344,3626,

%T 422,4092,1628,2212,602,10384,1227,2968,2548,6882,932,13552,1058,9006,

%U 3476,4816,3256,24087,1484,5908,4664,20768,1808,22792,1982,14694,11396,8428,2354,46706,3843,17178,7568

%N a(n) = Sum_{i|n, j|n, k|n} i*j*k/lcm(i,j,k).

%F If p is prime, a(p) = 4 + 3*p + p^2.

%t a[n_]:= Sum[i*j*k/LCM[i,j,k], {i, (d = Divisors[n])}, {j, d}, {k, d}]; Array[a, 50] (* _Amiram Eldar_, May 10 2021 *)

%o (PARI) a(n) = sumdiv(n, i, sumdiv(n, j, sumdiv(n, k, i*j*k/lcm([i, j, k]))));

%Y Cf. A060724, A344132, A344133, A344134.

%K nonn,mult

%O 1,2

%A _Seiichi Manyama_, May 10 2021