A358301 Main diagonal of array in A358298.
2, 6, 20, 60, 124, 252, 388, 652, 924, 1332, 1748, 2428, 2988, 3948, 4788, 5908, 7028, 8692, 9964, 12052, 13748, 16004, 18124, 21204, 23476, 26996, 29972, 33788, 37196, 42124, 45548, 51188, 55732, 61412, 66532, 73348, 78484, 86548, 92956, 100924, 107772, 117692, 124556, 135476, 144036
Offset: 0
Keywords
Crossrefs
Programs
-
Mathematica
A005728[n_] := 1 + Sum[EulerPhi[i], {i, 1, n}]; Amn[m_, n_] := Sum[If[GCD[i, j] == 1, 1, 0], {i, 1, m}, {j, 1, n}]; Dmn[m_, n_] := A005728[m] + A005728[n] + 2 Sum[d = GCD[u, v]; If[d >= 1, (u+v)*EulerPhi[d]/d, 0], {u, 1, m}, {v, 1, n}] - 2*Amn[m, n]; Table[Dmn[n, n], {n, 0, 44}] (* Jean-François Alcover, Apr 18 2023, after Maple code in A358298 *)