51nod 1675 序列变换 (简单莫比乌斯）

lyk有两序列a和b。
lyk想知道存在多少对x,y，满足以下两个条件。
1：gcd(x,y)=1。
2： abxabx = baybay 。

例如若a={1,1,1}，b={1,1,1}。那么存在7对，因为除了x=2,y=2或x=3,y=3外都满足条件。
Input
第一行一个数n(1<=n<=100000)。
接下来一行n个数，表示ai(1<=ai<=n)。
接下来一行n个数，表示bi(1<=bi<=n)。
Output
一行表示答案
Sample Input
3
1 1 1
1 1 1
Sample Output
7

就是利用莫比乌斯求出gcd=1的种数，难点在于ologn的统计吧，直接利用莫比乌斯跳就好，因为值不超过1e5

#include <bits/stdc++.h>
using namespace std;
const int N = 1e5+10;
typedef long long ll;
int a[N],b[N];
int viss[N];
int n;
bool vis[N+10];  
int prime[N+10],mu[N+10];  
int cnt;
void Init(){    
    memset(prime,0,sizeof(prime));  
    memset(mu,0,sizeof(mu));  
    memset(vis,0,sizeof(vis));  
    mu[1] = 1;  
    cnt = 0;  
    for(int i=2; i<N; i++){  
        if(!vis[i]){  
            prime[cnt++] = i;  
            mu[i] = -1;  
        }  
        for(int j=0; j<cnt&&i*prime[j]<N; j++){  
            vis[i*prime[j]] = 1;  
            if(i%prime[j]) mu[i*prime[j]] = -mu[i];  
            else{  
                mu[i*prime[j]] = 0;  
                break;  
            }  
        }  
    }
}

ll cal(int x)
{
    ll res=0;
    for(int i=x;i<=n;i+=x)
        viss[a[b[i]]]++;
    for(int i=x;i<=n;i+=x)
        res+=viss[b[a[i]]];
    for(int i=x;i<=n;i+=x)
        viss[a[b[i]]]--;
    return res;
}

int main()
{
    Init();
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d",&a[i]);
    for(int i=1;i<=n;i++)
        scanf("%d",&b[i]);
    ll res=0;
    for(int i=1;i<=n;i++)
    {
        res+=1ll*mu[i]*cal(i);
    }
    printf("%lld\n",res );
}