A372939 a(n) = Sum_{1 <= x_1, x_2, ... , x_n <= n} gcd(x_1, x_2, ... , x_n, n)^n.
1, 7, 53, 736, 6249, 185039, 1647085, 66912256, 1162222101, 39980466703, 570623341221, 53487845236736, 605750213184505, 44446670856053599, 1751575500459583437, 92232594468640915456, 1654480523772673528353, 236077997761634282622495
Offset: 1
Keywords
Programs
-
PARI
a(n) = sumdiv(n, d, moebius(n/d)*d^n*numdiv(d));
Formula
a(n) = Sum_{d|n} mu(n/d) * d^n * tau(d), where mu is the Moebius function A008683.