OFFSET
1,2
FORMULA
G.f.: (-1/2) * Sum_{k>0} phi(3*k) * log(1-3*x^k)/(3*k).
MATHEMATICA
a[n_] := DivisorSum[n, 3^(n/#-1)*EulerPhi[3*#]/(2*n) &]; Array[a, 30] (* Amiram Eldar, Jul 14 2023 *)
PROG
(PARI) a(n) = sumdiv(n, d, 3^(n/d-1)*eulerphi(3*d))/(2*n);
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jul 13 2023
STATUS
approved