|
MATHEMATICA
|
a[n_] := (8 n^2)^(-1) Sum[If[Mod[n, c] == 0, Sum[If[Mod[n, d] == 0, EulerPhi[c] EulerPhi[d] 2^(n^2/ LCM[c, d]), 0], {d, 1, n}], 0], {c, 1, n}] + (4 n)^(-1) Sum[If[Mod[n, d] == 0, EulerPhi[d] 2^(n^2/d), 0], {d, 1, n}] + If[Mod[n, 2] == 1, (4 n)^(-1) Sum[If[Mod[n, d] == 0 && Mod[d, 2] == 1, EulerPhi[d] (2^(n (n + 1)/(2 d)) - 2^(n^2/d)), 0], {d, 1, n}], (8 n)^(-1) Sum[If[Mod[n, d] == 0 && Mod[d, 2] == 1, EulerPhi[d] (2^(n^2/(2 d)) + 2^(n (n + 2)/(2 d)) - 2 2^(n^2/d)), 0], {d, 1, n}]] + (1/2) If[Mod[n, 2] == 1, 2^((n^2 - 3)/2), 7 2^(n^2/2 - 4)] + (4 n)^(-1) Sum[If[Mod[n, d] == 0, EulerPhi[d] 2^(n (n + d - 2 IntegerPart[d/2])/(2 d)), 0], {d, 1, n}] + If[Mod[n, 2] == 1, 2^((n^2 - 5)/4), 5 2^(n^2/4 - 3)];
|