51nod 1675 序列变换

莫比乌斯反演的题目。
设
$f(k)$为满足$a_{b_x}=b_{a_y}$的情况下，$gcd(x,y)=k$的数量。
$F(k)$为满足$a_{b_x}=b_{a_y}$的情况下，$gcd(x,y)=k的倍数$的数量。
则

$$F(k)=\sum_{k|d} f(d)$$
由莫比乌斯反演第二种形式，得到
$$f(k)=\sum_{k|d} u(\frac{d}{k})*F(d)$$
当 $k=1$时，
$$f(1)=\sum_{d=1}^n u(d)*F(d)$$
$u(d)$通过线性筛实现， $F(d)$可以通过平均复杂度 $O(\sqrt n)$的暴力实现。
总时间复杂度 $O(n \sqrt n)$

//#include<bits/stdc++.h>
#include<stdio.h>
#include<cstring>
#include<cmath>
#include<iostream>
using namespace std;

const int MAXN=100100;
int p[MAXN],mu[MAXN],sum[MAXN],f[MAXN];
int a[MAXN],b[MAXN];
int n;

long long cal(int x)
{
    long long i,ret=0;
    for(i=x;i<=n;i+=x)
        sum[a[b[i]]]++;
    for(i=x;i<=n;i+=x)
        ret+=sum[b[a[i]]];
    for(i=x;i<=n;i+=x)
        sum[a[b[i]]]--;
    return ret;
}

int main()
{
    int pcnt,i,j;
    long long ans;
    pcnt=0;
    memset(p,0,sizeof(p));
    memset(f,0,sizeof(f));
    memset(sum,0,sizeof(sum));
    mu[1]=1;
    for(i=2;i<MAXN;i++)
    {
        if(!f[i])
        {
            p[++pcnt]=i;
            mu[i]=-1;
        }
        for(j=1;p[j]*i<MAXN;j++)
        {
            f[p[j]*i]=1;
            if(i%p[j]==0)
            {
                mu[p[j]*i]=0;
                break;
            }
            mu[p[j]*i]=-mu[i];
        }
    }
    while(~scanf("%d",&n))
    {
        for(i=1;i<=n;i++)
            scanf("%d",&a[i]);
        for(i=1;i<=n;i++)
            scanf("%d",&b[i]);
        ans=0;
        for(i=1;i<=n;i++)
        {
            if(mu[i])
                ans+=mu[i]*cal(i);
        }
        printf("%lld\n",ans);
    }
}