ICode9

考虑把矩阵消成上三角然后求对角线的值。 可以发现每一行只会消掉自己的倍数行，且系数为 \(1\)。 假设第 \(n\) 行 \(n\) 列的元素是 \(f[n]\)，有： \[f[n]=n^k-\sum_{d\mid n,d\ne n}f[d] \]\[f * 1=id^k \]\[f=id^k * \mu \]考虑每个质数幂处的这玩意儿是好算的，而且是考虑答案的乘积