POJ2369 Permutations 置换群循环节

题目链接：POJ2369

题目大意：求循环节长度的最小公倍数

如果对这个知识点不是很懂，或者对这个知识点不明白，可以看我之前写的，传送门

/*
2017年8月2日10:10:54
POJ2369 循环节LCM 
AC
*/
#include<stdio.h>
#include<algorithm>
#include<string.h>
using  namespace std; 
typedef long long ll;
const int maxn=1010;
int	a[maxn],b[maxn],c[maxm],vis[maxn],l[maxn];
int n,tot,ans;

/*gcd*/
ll gcd(ll a,ll b){  
    return b>0?gcd(b,a%b):a;  
}

/*lcm*/
/*先除再乘，否则先乘可能会爆掉*/ 
ll lcm(ll x,ll y){
	return x/gcd(x,y)*y;
} 

void dfs(int u){
	for(int i=0;i<n;i++){
		if(b[i]==u&&!vis[i]){
			c[tot]++;
			vis[i]=1;
			dfs(a[i]);
		}
	} 
}

int main(){
	memset(vis,0,sizeof(vis)); 
	scanf("%d",&n);
	for(int i=0;i<n;i++){
		scanf("%d",&a[i]);
		b[i]=a[i];
	} 
	tot=0;
	sort(a,a+n);
	for(int i=0;i<n;i++){
		if(vis[i]) continue;
		else{
			c[tot]=1;
			vis[i]=1;
			dfs(a[i]);
			tot++;
		} 
	}
	/*
	这部分代码的就是求所有循环节的公倍数
	用l[maxn] 数组来存储 
	*/
	if(tot==1)	ans=c[0];
	else{
		l[0]=c[0];
		for(int i=1;i<tot;i++){
			l[i]=lcm(l[i-1],c[i]);
		}
		ans=l[tot-1];
	}
	printf("%d\n",ans);
	return 0;
} 

AC代码