hdu6038[找规律+循环节] 2017多校1

/*hdu6038[找规律+循环节] 2017多校1*/
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const double eps=1e-;
const int inf=0x3f3f3f3f;
const int mod=;
int n,m,lop,kase=,maxlop;
int a[],b[],mark[];
LL Arep[],Brep[];
LL fastpow(LL a,LL b){
    LL ret=,x=a;
    while(b){
        if(b&){
            ret=(ret*x)%mod;
        }
        x=(x*x)%mod;
        b>>=;
    }
    return ret;
}
void dfs(int now,int dps){
    if(mark[now]){
        Arep[dps]++;
        return;
    }
    mark[now]=;
    dfs(a[now],dps+);
}
void dfs2(int now,int dps){
    if(mark[now]){
        Brep[dps]++;
        return;
    }
    mark[now]=;
    dfs2(b[now],dps+);
}
void solve(){
    maxlop=;
    for(int i=;i<n;i++){
        if(!mark[a[i]]){
            dfs(i,);
        }
    }
    memset(mark,,sizeof(mark));
    for(int i=;i<m;i++){
        if(!mark[b[i]]){
            dfs2(i,);
        }
    }
    memset(mark,,sizeof(mark));
    LL ans=,sum=;
    for(int i=;i<=n;i++){
        if(Arep[i]){
            ans=;
            for(int j=;j*j<=i;j++){
                if(i%j==){
                    ans=(ans+(Brep[j]%mod)*(j%mod))%mod;
                    if(j*j!=i){
                        ans=(ans+(Brep[i/j]%mod)*((i/j)%mod))%mod;
                    }
                }
            }
            sum=(sum*fastpow(ans,Arep[i]))%mod;
        }
    }
    printf("Case #%d: %lld\n",kase++,sum);
}
int main(){
    //freopen("1006.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    while(~scanf("%d%d",&n,&m)){
        memset(mark,,sizeof(mark));
        memset(Arep,,sizeof(Arep));
        memset(Brep,,sizeof(Brep));
        for(int i=;i<n;scanf("%d",&a[i]),i++);
        for(int i=;i<m;scanf("%d",&b[i]),i++);
        solve();
    }
    return ;
}