poj 3101 Astronomy (java 分数的最小公倍数 gcd)

题目链接

要用大数，看了别人的博客，用java写的。

题意：求n个运动周期不完全相同的天体在一条直线上的周期。

分析：两个星球周期为a,b。则相差半周的长度为a*b/(2*abs(a-b)),对于n个只需求这n个

分数的最小公倍数即可。

分数的最小公倍数 = 分子的最小公倍数/分母的最大公约数

 

import java.util.*;
import java.math.*;
public class Main {
    public static int [] t = new int [1200];
    public static int [] tt = new int [1200];
    public static BigInteger [] fz = new BigInteger [1200];
    public static BigInteger [] fm = new BigInteger [1200];
    
    public static int gcd(int a, int b){
        return b==0?a:gcd(b, a%b);
    }
    public static void main(String args[]){
    int n, m;
     Scanner cin=new Scanner(System.in);
     n = cin.nextInt();
     for(int i = 0; i < n; i++)
         t[i] = cin.nextInt();
     Arrays.sort(t, 0, n);   //对数组排序
     m = 0;
     tt[m++] = t[0];
     for(int i = 1; i < n; i++)
         if(t[i]!=t[i-1])
             tt[m++] = t[i];
     for(int i = 1; i < m; i++){
         int a = tt[i]*t[0];
         int b = 2*(tt[i]-tt[0]);
         int g = gcd(a, b);
         fz[i] = BigInteger.valueOf(a/g);
         fm[i] = BigInteger.valueOf(b/g);
     }
     BigInteger t1, t2;
     t1 = fz[1]; t2 = fm[1];
     for(int i = 2; i < m; i++)
     {
         BigInteger aa = t1.multiply(fz[i]);
         BigInteger gg = t1.gcd(fz[i]);
         t1 = aa.divide(gg);
         t2 = t2.gcd(fm[i]);
     }
     System.out.println(t1+" "+t2);
    }
}