51nod 1675 (莫比乌斯反演)

我们先不考虑 $a_{b_i}=b_{a_j}$ 这种情况，那么就很裸的莫比乌斯反演了。
设 $$f(x)=\sum _{i=1}^n\sum _{j=1}^n[\gcd (i,j)==x]$$
在设 $$g(x)=\sum _{x|d}f(d)$$
$$=\sum _{x|d}\sum _{i=1}^n\sum _{j=1}^n[\gcd (i,j)==d]$$
$$=\sum _{i=1}^n\sum _{j=1}^n[x|\gcd (i,j)]$$
$$=\sum _{x|i}^n\sum _{x|j}^{n}1$$
就得到如上式子，
那如果在考虑 $a_{b_i}=b_{a_j}$这个因素，只需加一个限制条件
$$f(x)=\sum _{i=1}^n\sum _{j=1}^n[\gcd (i,j)==x][a_{b_i}=b_{a_j}]$$
套用上面的方式，得到：$$g(x)=\sum _{x|i}^n\sum _{x|j}^n[a_{b_i}=b_{a_j}]$$
对于这个式子，可以 $log$得到

然后反演回去：$$f(x)=\sum _{x|d}u(\frac{d}{x})\ast g(d)$$
这里x 等于1，所以最终答案为：
$$f(1)=\sum _{d=1}^{n}u(d)\ast g(d)$$
$$=\sum _{d=1}^{n}u(d)\ast \sum _{d|i}^n\sum _{d|j}^n[a_{b_i}=b_{a_j}]$$

最终复杂度 $O(nlogn)$