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

题目链接

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

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

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

分数的最小公倍数即可。

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

 1 import java.util.*;
 2 import java.math.*;
 3 public class Main {
 4     public static int [] t = new int [1200];
 5     public static int [] tt = new int [1200];
 6     public static BigInteger [] fz = new BigInteger [1200];
 7     public static BigInteger [] fm = new BigInteger [1200];
 8
 9     public static int gcd(int a, int b){
10         return b==0?a:gcd(b, a%b);
11     }
12     public static void main(String args[]){
13     int n, m;
14      Scanner cin=new Scanner(System.in);
15      n = cin.nextInt();
16      for(int i = 0; i < n; i++)
17          t[i] = cin.nextInt();
18      Arrays.sort(t, 0, n);   //对数组排序
19      m = 0;
20      tt[m++] = t[0];
21      for(int i = 1; i < n; i++)
22          if(t[i]!=t[i-1])
23              tt[m++] = t[i];
24      for(int i = 1; i < m; i++){
25          int a = tt[i]*t[0];
26          int b = 2*(tt[i]-tt[0]);
27          int g = gcd(a, b);
28          fz[i] = BigInteger.valueOf(a/g);
29          fm[i] = BigInteger.valueOf(b/g);
30      }
31      BigInteger t1, t2;
32      t1 = fz[1]; t2 = fm[1];
33      for(int i = 2; i < m; i++)
34      {
35          BigInteger aa = t1.multiply(fz[i]);
36          BigInteger gg = t1.gcd(fz[i]);
37          t1 = aa.divide(gg);
38          t2 = t2.gcd(fm[i]);
39      }
40      System.out.println(t1+" "+t2);
41     }
42 }

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