【置换群】【枚举约数】hdu6038 Function

#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
#define MOD 1000000007ll
int n,m,a[100010],b[100010],zu;
bool vis[100010];
//ll Quick_Pow(ll a,ll p){
//	if(!p){
//		return 1;
//	}
//	ll res=Quick_Pow(a,p>>1);
//	res=res*res%MOD;
//	if(p&1ll){
//		res=(a%MOD*res)%MOD;
//	}
//	return res;
//}
int nowsiz,cnts[100010];
void dfs(int U){
	vis[U]=1;
	++nowsiz;
	if(!vis[b[U]]){
		dfs(b[U]);
	}
}
void df2(int U){
	vis[U]=1;
	++nowsiz;
	if(!vis[a[U]]){
		df2(a[U]);
	}
}
int main(){
	while(scanf("%d%d",&n,&m)!=EOF){
		memset(cnts,0,sizeof(cnts));
		++zu;
		for(int i=0;i<n;++i){
			scanf("%d",&a[i]);
		}
		for(int i=0;i<m;++i){
			scanf("%d",&b[i]);
		}
		for(int i=0;i<m;++i){
			if(!vis[i]){
				nowsiz=0;
				dfs(i);
				++cnts[nowsiz];
			}
		}
		memset(vis,0,sizeof(bool)*m);
		ll ans=1;
		for(int i=0;i<n;++i){
			if(!vis[i]){
				nowsiz=0;
				df2(i);
				int allcanb=0;
				for(int i=1;i*i<=nowsiz;++i){
					if(nowsiz%i==0){
						if(i*i==nowsiz){
							allcanb+=cnts[i]*i;
						}
						else{
							allcanb+=cnts[i]*i;
							allcanb+=cnts[nowsiz/i]*(nowsiz/i);
						}
					}
				}
				ans=ans*(ll)allcanb%MOD;
			}
		}
		printf("Case #%d: %lld\n",zu,ans);
		memset(vis,0,sizeof(bool)*n);
	}
	return 0;
}