幂一矩阵（大数）

幂一矩阵

天才少年的邻居 atm 最近学习了线性代数相关的理论，他对“矩阵”这个概念特别感兴趣。矩阵中有个概念叫做幂零矩阵。对于一个方阵 M ，如果存在一个正整数 k 满足 M^k = 0 ，那么 M 就是一个幂零矩阵。（^ 表示乘方）

atm 不满足幂零矩阵，他自己设想了一个幂一矩阵：对于一个方阵 M ，如果存在一个正整数 k 满足 M^k = I ，其中 I 是单位矩阵，那么 M 就是一个幂一矩阵。

atm 特别钟情于这样一种方阵：每行每列有且仅有一个 1 。经过 atm 不断实验，他发现这种矩阵都是幂一矩阵。

现在，他的问题是，给定一个满足以上条件的方阵，他想求最小的 k 是多少。

【输入格式】
第一行一个正整数 n ，表示矩阵大小是 n * n 。
接下来 n 行，每行两个正整数 i j 表示方阵的第 i 行第 j 列为 1。
1 <= i, j <= n 。
行号，列号都从1开始。

【输出格式】
一行。一个正整数，即题目中所说最小的 k 。

【样例输入】
5
3 1
1 2
4 4
2 3
5 5

【样例输出】
3

【数据范围】
对于 30% 的数据满足 n <= 10
对于 60% 的数据答案不超过 10^18
对于 100% 的数据满足 n <= 10000

资源约定：
峰值内存消耗（含虚拟机） < 256M
CPU消耗 < 1000ms

大数做法暂且不会，日后再补，这里贴一个能过50%测试点的方法。

思路：假设题目给的的矩阵为A，A*A得到的矩阵为B，设A第一行的1在第y1列（y1!=1），而y1行1的位置为y2，那么得到的B的第一行的1的位置为y2，此时B的第y2行1的位置为y3，那么A和B相乘得到的矩阵第一行1的位置为y3，依次类推，y1,y2…yn中一定有一个为1，此时进行的阶乘数即为第一行复原所需最小的k。
其他行也以此类推，有一个ki，那么最后结果就是这些k的最小公倍数。

#include<iostream>
using namespace std;
typedef long long ll;
#define N 10005
int a[N];    //a[i]代表，第i行的1的位置（在哪一列）
int n;
int gcd(int x,int y)
{
	int r=x%y;
	while(r)
	{
		x=y;y=r;r=x%y;
	}
	return y;
}
int main()
{
	int x,y;
	ll k=0,ans=1;
	cin>>n;
	for(int i=0;i<n;i++)
	{
		scanf("%d%d",&x,&y);
		a[x]=y;
	}
	for(int i=1;i<=n;i++)  //找出每一个ki
	{
		k=1,x=a[i];
		while(x!=i)
		{
			x=a[x];
			k++;
		}
		ans=ans*k/gcd(ans,k);
	}
	cout<<ans<<endl;
	return 0;
}