POJ 3101 Astronomy 【分数的最小公倍数+物理公式+java大数】

题目链接：http://poj.org/problem?id=3101

思路：

物理公式就不解释了 ：  ;

求 的最小公倍数，因为 可能为分数，所以需要求分数的最小公倍数；

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

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

因为题目的数据范围很大，所以用java处理会方便点；需要注意的是去重，因为速度相同的情况下，周期相同说明在同一轨道上，直接忽略不计就行了；

import java.util.*;
import java.text.*;
import java.math.*;
import java.io.*;

public class Main {
	public static void main(String[] args) {
		Scanner  cin = new Scanner(new BufferedInputStream(System.in));
		int N, a[] = new int [1005], b[] = new int[1005];
		while (cin.hasNext()) {
			N = cin.nextInt();
			for (int i = 0; i < N; ++i) a[i] = cin.nextInt();
			Arrays.sort(a, 0, N); int m = 1; b[0] = a[0]; 
			for (int i = 1; i < N; ++i) { 
				if (a[i] != a[i-1]) b[m++] = a[i]; // 去重
			}
			
			BigInteger lcm = null, g = null, G; // 数据并没有考虑如果N == 2，两者又相同的情况
			for (int i = 1; i < m; ++i) {
				
				BigInteger up = BigInteger.valueOf(b[i]*b[i-1]);
				BigInteger down = BigInteger.valueOf(2*(b[i]-b[i-1]));
				G = up.gcd(down);
				up = up.divide(G); down = down.divide(G);
				if (i == 1) {
					lcm = up; g = down;
				} else {
					lcm = lcm.multiply(up).divide(lcm.gcd(up));
					g = g.gcd(down);
				}
			}
			System.out.println(lcm + " " + g);
		}
	}
}