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