The motivation is this n by n by n matrix/sum in Mathematica 14:
In[512]:= Clear[a, n, k, m]
nn = 40;
g[n_] := Total[Divisors[n]*MoebiusMu[Divisors[n]]]
s1 = 12;
s2 = 12;
s3 = 12;
Sum[Sum[Sum[g[GCD[a, b, c]]/(N[a, 20]^s1*b^s2*c^s3), {a, 1, nn}], {b,
1, nn}], {c, 1, nn}]
N[Zeta[s1]*Zeta[s2]*Zeta[s3]/Zeta[s1 + s2 + s3 - 1], 20]
Out[518]= 1.000738441321476615
Out[519]= 1.0007384413214766152