OFFSET
1,2
MATHEMATICA
t[n_, k_] := (For[t1 = 0; d = 1, d <= n, d++, If[Mod[n, d] == 0, t1 = t1 + EulerPhi[d]*k^(n/d)]]; If[EvenQ[n], (t1 + (n/2)*(1+k)*k^(n/2))/(2*n), (t1 + n*k^((n+1)/2))/(2*n)]); a[n_] := Sum[t[n, k], {k, 1, n}]; Array[a, 20] (* Jean-François Alcover, Nov 02 2017, after Maple code for A081720 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, based on information supplied by Gary W. Adamson, Apr 05 2003
STATUS
approved