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

http://poj.org/problem?id=3101

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

这题我是看解题报告写的,没想到选用参照物，用到了物理中的角速度什么的。

因为n个天体的周期已知，那么它们的角速度为vi = 2*pi/Ti，若统一选第0个天体为参照物，那么其余天体的相对速度vi' = 

2*pi*(T0-Ti)/(T0*Ti)(把周期T相同的天体合为一个天体)。则与第0个天体角度相差180度的时间为ti = (T0*Ti)/((T0-Ti)*2)。那么求得所有ti的最小公倍数就是答案。

终于到重点了。ti作为分数，它们的最小公倍数定义为 : 所有分子的最小公倍数/所有分母的最大公约数数。

由于N太大，需要用大数处理。又各种百度java，终于A啦。 两点了，洗洗睡吧。

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

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){
		if(b == 0)
			return a;
		return 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); //java中对数组的排序方法
		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] * tt[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 = 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);
	}
	
}