%I #18 Jul 28 2022 03:39:45
%S 1,2,4,3,6,7,4,8,12,14,5,10,15,20,13,6,12,14,24,30,28,7,14,21,28,35,
%T 42,19,8,16,24,26,40,48,56,42,9,18,19,36,45,38,63,72,37,10,20,30,40,
%U 26,60,70,80,90,52,11,22,33,44,55,66,77,88,99,110,31,12,24
%N Triangle read by rows, n>=1, 1<=k<=n, T(n,k) = Sum_{c|n,d|k} phi(lcm(c,d)).
%C This is the lower triangular array of A216622, which is the main entry for this sequence.
%C T(n,1) = A000027(n).
%C T(n,n) = A062380(n).
%H Alois P. Heinz, <a href="/A216623/b216623.txt">Rows n = 1..141, flattened</a>
%e The first rows of the triangle are:
%e 1,
%e 2, 4,
%e 3, 6, 7,
%e 4, 8, 12, 14,
%e 5, 10, 15, 20, 13,
%e 6, 12, 14, 24, 30, 28,
%e 7, 14, 21, 28, 35, 42, 19,
%e 8, 16, 24, 26, 40, 48, 56, 42,
%e 9, 18, 19, 36, 45, 38, 63, 72, 37,
%p with(numtheory):
%p T:= (n, k)-> add(add(phi(ilcm(c, d)), c=divisors(n)), d=divisors(k)):
%p seq (seq (T(n, k), k=1..n), n=1..14); # _Alois P. Heinz_, Sep 12 2012
%t t[n_, k_] := Sum[ EulerPhi[ LCM[c, d]], {c, Divisors[n]}, {d, Divisors[k]}]; Table[t[n, k], {n, 1, 12}, {k, 1, n}] // Flatten (* _Jean-François Alcover_, Jul 23 2013 *)
%o (Sage) # uses[A216622]
%o for n in (1..9): [A216622(n,k) for k in (1..n)]
%Y Cf. A216620, A216621, A216622, A216624, A216625, A216626, A216627.
%K nonn,tabl
%O 1,2
%A _Peter Luschny_, Sep 12 2012
|