蓝桥真题-幂一矩阵

题目：

天才少年的邻居 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

得知题目给的起始矩阵即幂一矩阵除了1其余都是0，那么一个矩阵要通过自乘变成单位矩阵(第i行第i行的值为1)，通过手动运算几次会发现，每进行一次自乘，原矩阵有1的位置为2 3和3 1，那么自乘后有1的位置就是2 1，也就是说我们需要遍历每一行，查找哪些行进行运算后有1的位置会变成i i，并统计次数，则说明这些参与运算的行经过统计次数后是会变成i i位置为1的，因为要使所有行同时i i位置变为1，因此要算所有统计次数的最小公倍数，这就是答案了。

注意：光开long是不够的，算最小公倍数时两数相乘是会溢出long的，因此对于涉及最小公倍数运算的变量应该都开BigInteger才够。
 

如何用辗转相除法求最小公倍数？？？
被除数 / 除数 = 商 ...... 余数
6497 / 3869 = 1 ...... 2628
3869 / 2628 = 1 ...... 1241
2628 / 1241 = 2 ...... 146
1241 / 146 = 8 ...... 73
146 / 73 = 2 ...... 0
因此最大公约数为:73
最小公倍数=两数之积/最大公约数=6497*3869/73=25136893

import java.math.BigInteger;
import java.util.Scanner;

public class 幂一矩阵 {

    public static void main(String[] args) {
        Scanner in = new Scanner(System.in);
        int n = in.nextInt();
        //创建索引数组num[n]
        int[] num = new int[n];
        
        for(int i = 0;i<n;i++){
            //输入数组，注意要-1，因为题目说输入的行和列都是以1开始
            num[in.nextInt()-1] = in.nextInt()-1;
        }
        //这里vis数组为了加快循环
        boolean[] vis = new boolean[n];
        //max存储最终结果
        BigInteger max = BigInteger.ONE;
        //从原数组的第0行开始，寻找搜索
        for(int i = 0;i<n;i++){
            //如果当前行已经被搜索过，直接continue
            if(vis[i]){
                continue;
            }
            //count记录搜索次数
            BigInteger count = BigInteger.ONE;
            //s为每次搜索的数组下标索引
            int s = i;
            //当前第s行被搜索了，记录vis数组
            vis[s] = true;
            //当索引与值不相等时就继续循环
            while(num[s] != i){
                //每次循环一次count+1
                count = count.add(BigInteger.ONE);
                //当前索引对应的值是下一次的索引
                s = num[s];
                //当前索引对应的行被搜索过，记录vis数组
                vis[s] = true;
            }
            //搜索结束，判断搜索结果是否能被max整除
           
            max = (max.multiply(count)).divide(max.gcd(count));
        }
        //输出最终结果
        System.out.println(max);
        in.close();
    }
    private static long gcd(long max, long count) {
        //辗转相除法求最大公约数
        if(max%count==0){
            return count;
        }
        return gcd(count,max%count);
    }
    
   


}