题目链接:

id=3101">点击打开链接

题目大意:有n个行星,给出每个行星的旋转的周期。问最少多少时间后n个行星会在一条直线上,初始点在一起,不存在全部的行星都有同一个周期

如果A行星的周期是t1。B行星的周期是t2(t2>t1),要在一条直线上,一定会执行的相差半个周期的倍数,时间(t/t2 - t/t1) % (1/2) = 0。也就是t*(t1-t2)/(t1*t2)%(1/2) = 0,要是时间最小。所以也就是差出一个半周期。也就是t = (t2-t1)/(t2*t1*2)这个t也就是A。B执行到一条直线上的最小时间,我们能够求出其它全部行星和A行星的在一条直线的最小时间,然后求出这个时间的最小公倍数。也就是整体的时间。

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

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

import java.util.*;
import java.math.* ;
public class Main {
public static void main(String[] args) {
Scanner cin = new Scanner(System.in) ;
int n , i , cnt = 0 ;
int a[] = new int[1010];
int b[] = new int[1010] ;
BigInteger x , y , temp , u , v ;
n = cin.nextInt() ;
for(i = 0 ; i < n ; i++)
a[i] = cin.nextInt() ;
Arrays.sort(a,0,n) ;
b[cnt++] = a[0] ;
for(i = 1 ; i < n ; i++)
if( a[i] != b[cnt-1] )
b[cnt++] = a[i] ;
x = BigInteger.valueOf(b[1]*b[0]) ;
y = BigInteger.valueOf((b[1]-b[0])*2) ;
temp = x.gcd(y) ;
x = x.divide(temp) ;
y = y.divide(temp) ;
for(i = 2 ; i < cnt ; i++) {
u = BigInteger.valueOf(b[i]*b[0]) ;
v = BigInteger.valueOf((b[i]-b[0])*2) ;
temp = u.gcd(v) ;
u = u.divide(temp) ; v = v.divide(temp) ;
temp = x.gcd(u) ;
x = x.multiply(u).divide(temp) ;
y = y.gcd(v) ;
}
System.out.println(x + " " + y) ;
}
}

poj3101--Astronomy(分数的最小公倍数)的更多相关文章

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

    http://poj.org/problem? id=3101 大致题意:求n个运动周期不全然同样的天体在一条直线上的周期. 这题我是看解题报告写的,没想到选用參照物,用到了物理中的角速度什么的. 由 ...

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

    题目链接 要用大数,看了别人的博客,用java写的. 题意:求n个运动周期不完全相同的天体在一条直线上的周期. 分析:两个星球周期为a,b.则相差半周的长度为a*b/(2*abs(a-b)),对于n个 ...

  3. hdu 1713求分数的最小公倍数

    题意中的圈数和天数说反了 #include<stdio.h> __int64 gcd(__int64 a,__int64 b) {/* 比如4/3 3/5 通分20/15 9/15 所以这 ...

  4. poj 3101Astronomy(圆周追击+分数最小公倍数)

    /* 本题属于圆周追击问题: 假设已知两个圆周运动的物体的周期分别是a ,b, 设每隔时间t就会在同一条直线上 在同一条直线上的条件是 角度之差为 PI ! 那么就有方程 (2PI/a - 2PI/b ...

  5. poj 3101 Astronomy

    2个星球周期为a,b.则相差半周的长度为a*b/(2*abs(a-b)),对于n个只需求这n个 分数的最小公倍数即可! 公式: 分数的最小公倍数 = 分子的最小公倍数/分母的最大公约数 由于涉及到大数 ...

  6. HDU 1713 相遇周期 (最小公倍数)

    题意:... 析:求周期就是这两个分数的最小公倍数,可以先通分,再计算分子的最小倍数. 代码如下: #pragma comment(linker, "/STACK:1024000000,10 ...

  7. acm数学(转)

    这个东西先放在这吧.做过的以后会用#号标示出来 1.burnside定理,polya计数法    这个大家可以看brudildi的<组合数学>,那本书的这一章写的很详细也很容易理解.最好能 ...

  8. [转] POJ数学问题

    转自:http://blog.sina.com.cn/s/blog_6635898a0100magq.html 1.burnside定理,polya计数法 这个大家可以看brudildi的<组合 ...

  9. ACM数学

     1.burnside定理,polya计数法 这个专题我单独写了个小结,大家可以简单参考一下:polya 计数法,burnside定理小结 2.置换,置换的运算 置换的概念还是比较好理解的,< ...

随机推荐

  1. 串口通讯 ADC0804 数码管

    #include<reg52.h> #include<stdio.h> #include<intrins.h> #define uchar unsigned cha ...

  2. Arduino Yun高速新手教程(大学霸内部资料)

    Arduino Yun高速新手教程(大学霸内部资料) 本资料为国内第一本Arduino Yun教程.具体解说Arduino Yun的基本结构.开发环境.系统配置.并着力解说关键功能--Bridge.最 ...

  3. Nginx 做系统的前端反向proxy

    Nginx是一款很优秀的基于event的webserver.吞吐量大.占用资源少,只是文档就很让人郁闷了,免费的Nginx和收费的Nginx+的文档共用一份,配置完之后才发现免费的Nginx启动某些命 ...

  4. windows、linux劫持技术

    windows系统以下能够利用detours劫持 realse  模式劫持,调试的程序不能够 函数劫持能够实现的效果. 函数的劫持原理 我们怎样实现-detours detours是微软亚洲研究院出品 ...

  5. Android - 使用JD-GUI反编译Android代码

    使用JD-GUI反编译Android代码 本文地址: http://blog.csdn.net/caroline_wendy Android程序出现Bug时,须要依据Bug寻找问题出错的地方; 须要使 ...

  6. 在eclipse中关联android源代码

    1打包源代码成jar: 1 新建一个java项目 2  import  想打包的源代码文件 3 export 这个文件 : 选择java->jar file .  这里会让你选择输出路径 2 加 ...

  7. MFC exe使用C++ dll中的std::string 崩溃

    VC6中 MFC exe中 new  纯C++ dll dll 崩溃 我把纯C++的 dll,用/MTd 换成/MDd.就能够了

  8. 51NOD 1686 第K大区间 二分

    第k大区间   定义一个区间的值为其众数出现的次数.现给出n个数,求将所有区间的值排序后,第K大的值为多少. Input   第一行两个数n和k(1<=n<=100000,k<=n* ...

  9. hdoj--3440--House Man(差分约束)

    House Man Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total ...

  10. Oracle 查看锁定表用户的SQL

    #描述:当执行如下SQL select * From XXX for Updata 发现表“XXX”被锁定 #只需要执行如下SQL 查看锁定该表的用户 SELECT object_name, mach ...