hdu 6038 Function

转自：http://blog.csdn.net/cysjiang/article/details/76090902

题意：
给你一个a序列，代表0到n-1的排列；一个b序列代表0到m-1的排列。问你可以找出多少种函数关系，满足f(i)=b[f(a[i])];
分析：这个主要是找循环节
比如说：如果 a 序列是 2 0 1 那么我们可以发现

f(0) = b[f(a[0])] = b[f(2)]
f[1] = b[f(a[1])] = b[f(0)]
f[2] = b[f(a[2])] = b[f(1)]
那么f(0) f(1) f(2) 也是循环的 如果想找出这样的函数，必须值域里也存在同样长度的循环节或者存在其约数长度的循环节。
那么就是找两个序列的循环节，对于定义域里面的每一个循环节都找出来有多少种和他对应的。最后乘起来就是答案了

#include <bits/stdc++.h>
using namespace std;
const int N = 1e6+100;
int a[N],b[N],num1[N],num2[N];
typedef long long ll;
const int mod = 1e9+7;
int main()
{
    int n,m;
            int cas=1;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        memset(num2,0,sizeof(num2));
        for(int i=0;i<n;i++)
            scanf("%d",&a[i]);
        for(int i=0;i<m;i++)
            scanf("%d",&b[i]);
        int totfa=0;
        for(int i=0;i<n;i++)
        {
            int k=i;
            int tot=0;
            while(a[k]!=-1)
            {
                int t=k;
                k=a[k];
                a[t]=-1;
                tot++;
            }
            if(tot)num1[totfa++]=tot;
        }
        for(int i=0;i<m;i++)
        {
            int k=i;
            int tot=0;
            while(b[k]!=-1){
                int t=k;
                k=b[k];
                b[t]=-1;
                tot++;
            }
            if(tot)num2[tot]++;
        }
        long long total=1;
        for(int i=0;i<totfa;i++)
        {
        ll tol=0;
        for(int j=1;j*j<=num1[i];j++)
        {
            if(j*j==num1[i])
            {
                tol=(tol+num2[j]*j)%mod;
            }
            else if(num1[i]%j==0){
                tol=(tol+num2[j]*j%mod+num2[num1[i]/j]*(num1[i]/j)%mod)%mod;
            }
        }
        total*=tol;
        total%=mod;
        }
        printf("Case #%d: %lld\n",cas++,total);
    }
}