题目:
天才少年的邻居 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);
}
}