login
A295521
a(n) = (1/n) * Sum_{d|n} mu(n/d)*(4^d - 3^d - 2^d + 1).
2
0, 2, 10, 39, 150, 545, 2010, 7320, 26880, 98775, 365010, 1353185, 5038950, 18830145, 70623958, 265737270, 1002976350, 3796197160, 14406059010, 54801140307, 208932573650, 798218035245, 3055417070010, 11716354754030, 45002103387120, 173117601112575
OFFSET
1,2
LINKS
FORMULA
a(n) = A027377(n) - A027376(n) - A001037(n) for n > 1.
MATHEMATICA
a[n_] := DivisorSum[n, MoebiusMu[n/#] * (4^# - 3^# - 2^# + 1) &] / n; Array[a, 26] (* Amiram Eldar, Oct 04 2023 *)
PROG
(PARI) {a(n) = sumdiv(n, d, moebius(n/d)*(4^d-3^d-2^d+1))/n}
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Nov 23 2017
STATUS
approved