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

  1 /*
  2    本题属于圆周追击问题:
  3      假设已知两个圆周运动的物体的周期分别是a ,b, 设每隔时间t就会在同一条直线上 
  4      在同一条直线上的条件是 角度之差为 PI !
  5      那么就有方程 (2PI/a - 2PI/b)* t=PI  所以就有 t=ab/(2|a-b|);
  6      如果有多个物体, 就会有多个t值,所以每隔 所有 t值的最小公倍数的时间所有的物体就会在同一直线上!
  7      
  8      另外:如果分数的分子分别是 a1, a2, ...., 和 b1, b2, ....
  9      那么所有分数的最小公倍数就是lcm(a1, a2, ...)/gcd(b1, b2,....);
 10      
 11      再有:如何求多个数的最小公倍数呢?
 12      根据数论,每一个数都可以表示成素数的乘积的形式! 
 13      令p[i]存储素数,将a1,a2,...分别整除以p[i],直到除尽!并记录除以每个p[i]时的个数temp;
 14      并更新该个数的最大值cnt[i]=max(temp, cnt[i]); 
 15      
 16      最后cnt[i]个p[i]分别相乘得到最终的结果就是所有数的最小公倍数! 
 17 */
 18 #include<iostream>
 19 #include<cstring>
 20 #include<cmath>
 21 #include<cstdio>
 22 #define M 10005
 23 #define N 1005
 24 using namespace std;
 25 typedef long long LL;
 26 LL p[M];
 27 bool isP[M];
 28 LL cnt[M];
 29 LL q[N];
 30 LL ans[N], endx;
 31 LL top;
 32 
 33 void bigN(){//大数据的处理
 34    LL c=0;
 35    endx=0;
 36    ans[0]=1;
 37    for(LL i=0; i<top; ++i)
 38       for(LL j=0; j<cnt[i]; ++j){
 39           for(LL k=0; k<=endx; ++k){
 40              ans[k]=ans[k]*p[i] + c;
 41              c=ans[k]/10000;
 42              ans[k]%=10000;
 43           }
 44           if(c>0){
 45              ans[++endx]=c;
 46              c=0;
 47           }
 48       }
 49 }
 50 
 51 void isPrime(){
 52    LL i, j;
 53    isP[1]=1;
 54    for(i=2; i<M; ++i){
 55        if(!isP[i])  p[top++]=i;
 56        for(j=0; j<top && i*p[j]<M; ++j){
 57            isP[i*p[j]]=1;
 58            if(i%p[j]==0) break;
 59        }
 60    }
 61 }
 62 
 63 void solve(LL k){
 64    for(LL i=0; i<top && p[i]<=k; ++i){
 65        LL tmp=0;
 66        while(k%p[i]==0){
 67           ++tmp;
 68           k/=p[i];
 69        }
 70        
 71        if(tmp>cnt[i])
 72           cnt[i]=tmp;
 73    }
 74 }
 75 
 76 LL gcd(LL a, LL b){
 77    while(b){
 78       LL r=a%b;
 79       a=b;
 80       b=r;
 81    }
 82    return a;
 83 }
 84 
 85 int main(){
 86    LL n;
 87    isPrime();
 88    while(scanf("%lld", &n)!=EOF){
 89           memset(cnt, 0, sizeof(cnt));
 90           scanf("%lld", &q[0]); 
 91        for(LL i=1; i<n; ++i){
 92            scanf("%lld", &q[i]);
 93            LL tmp=q[0]-q[i]>0 ? q[0]-q[i] : q[i]-q[0];
 94            if(tmp!=0){
 95               LL GCD=gcd(tmp, q[0]*q[i]);
 96               solve(q[0]*q[i]/GCD);
 97               q[i]=tmp/GCD;
 98            }
 99            else q[i]=0;
100        } 
101        
102        LL ans2=0;
103        for(LL i=1; i<n; ++i)
104           ans2=gcd(ans2, q[i]);
105        if(cnt[0]>0)//除以2 
106            --cnt[0];
107        else ans2*=2; 
108      
109        bigN();
110        if(ans2==0){
111            endx=0;
112            ans[endx]=0;
113        }
114        printf("%lld", ans[endx]);
115        for(int i=endx-1; i>=0; --i)
116           printf("%04lld", ans[i]);
117        printf(" %lld\n", ans2);
118    }
119    return 0;
120 } 

 1 //用java爽一下,处理大数
 2 import java.util.Scanner;
 3 import java.util.Arrays;
 4 import java.math.*;
 5 import java.io.BufferedInputStream;
 6 class Main{
 7    static int[] tt = new int[1005];
 8    static int n;
 9    static int top=0;
10    static boolean[] flag = new boolean[10005];
11    static int[] p = new int[10005];
12    static int[] q = new int[10005];
13    static int[] aa = new int[10005];
14    public static void isprime(){
15        int i, j;
16        Arrays.fill(flag, false);
17        for(i=2; i<=10000; ++i){
18            if(!flag[i]) p[top++]=i;
19            for(j=0; j<top && i*p[j]<=10000; ++j){
20               flag[i*p[j]]=true;
21               if(i%p[j]==0) break;
22            }     
23        }
24        --top;
25        flag[1]=true;
26    }
27    public static void solve(int k){
28         int i, cnt;
29         for(i=0; i<=top && p[i]<=k; ++i){
30            cnt=0;
31            while(k%p[i]==0){
32                 ++cnt;
33                 k=k/p[i];
34              }
35            if(cnt>aa[i])
36               aa[i]=cnt;
37         }
38    }
39 
40    public static int gcd(int a, int b){
41        while(b!=0){
42           int r=a%b;
43           a=b;
44           b=r;
45        }
46        return a;
47    }
48    
49    
50    public static  void main(String[] args){
51        isprime();
52        Scanner input = new Scanner(new BufferedInputStream(System.in));
53        n=input.nextInt();
54        q[0]=input.nextInt();
55        for(int i=1; i<n; ++i){
56            q[i]=input.nextInt();
57            int temp=Math.abs(q[0]-q[i]);
58            if(temp!=0){
59               int GCD=gcd(temp, q[0]*q[i]);
60               solve(q[0]*q[i]/GCD);
61               q[i]=temp/GCD;
62            }
63            else q[i]=0;
64        }
65        
66        BigInteger bigN = BigInteger.ONE;
67        for(int i=0; i<=top; ++i){
68            for(int j=0; j<aa[i]; ++j)
69               bigN=bigN.multiply(BigInteger.valueOf(p[i]));
70        }
71        for(int i=0; i<=top; ++i)
72           if(aa[i]!=0)
73             System.out.println(p[i]+" "+aa[i]);
74        int ans=0;
75        for(int i=1; i<n; ++i){
76            ans=gcd(ans, q[i]);
77        }
78        if(aa[0]>0)
79             bigN=bigN.divide(BigInteger.valueOf(2));
80        else ans*=2;
81        if(ans==0)
82             bigN=BigInteger.ZERO;
83        System.out.println(bigN+" "+ans);
84    }
85 }

 



 

转载于:https://www.cnblogs.com/hujunzheng/p/3891180.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/531807.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

广东省计算机应用考试试题,2015广东省计算机等级考试试题 二级C试题最新考试试题库...

1、在计算机的应用中&#xff0c;“DSS”表示( B )A、管理信息系统 B、决策支持系统C、办公自动化 D、人工智能2、计算机病毒是可以造成机器故障的( D )A、一种计算机设备 B、一种计算机芯片C、一种计算机部件 D、一种计算机程序3、防范病毒的有效手段&#xff0c;不正确的是( …

mysql查找大小写_mysql查询不区分大小写

摘自&#xff1a;http://www.jb51.net/article/70884.htm当我们输入不管大小写都能查询到数据&#xff0c;例如&#xff1a;输入 aaa 或者aaA ,AAA都能查询同样的结果&#xff0c;说明查询条件对大小写不敏感。解决方案一&#xff1a;于是怀疑Mysql的问题。做个实验&#xff1a…

poj 2492A Bug's Life(并查集)

1 /*2 目大意&#xff1a;输入一个数t&#xff0c;表示测试组数。然后每组第一行两个数字n,m&#xff0c;n表示有n只昆虫&#xff0c;编号从1—n,m表示下面要输入m行交配情况&#xff0c;每行两个整数&#xff0c;表示这两个编号的昆虫为异性&#xff0c;要交配。3 要求统计交配…

mysql 有外键 怎么插入数据_外键约束的表怎么插入数据

有外键的情况应该先添加主表数据&#xff0c;再添加副表数据。如&#xff1a;有以下两张表班级表&#xff1a;CLASSID NAME1 一班2 二班学生表&#xff1a;SID NAME CLASSID1 张三 12 李四 13 王五 2其中学生表中的CLASSID是班级表CLASSID的外键。现在要求在学生表中添加一条SI…

计算机在岗位上的应用,计算机岗位应用论文.doc

计算机岗位应用论文计算机岗位应用论文计算机岗位应用论文计算机岗位应用论文计算机岗位应用论文进入到新世纪以来,随着我国国民经济水平的提升,我国的计算机也得到了迅速的发展,计算机应用技术也已经广泛的应用到了各个行业中,并且计算机应用技术对于促进这些行业的快速发展也…

poj1703Find them, Catch them(并查集以及路径压缩)

1 /*2 题目大意&#xff1a;有两个不同的黑帮&#xff0c;开始的时候不清楚每个人是属于哪个的&#xff01;3 执行两个操作4 A a, b回答a, b两个人是否在同一帮派&#xff0c;或者不确定5 D a, b表示a, b两个人不在同一个帮派6 7 思路&#xff1a;利用并查集将相同帮…

计算机课程word教学,Word教学方法及使用技巧

Word教学方法及使用技巧来源:用户上传作者: 李志爱【摘 要】Word是一款实用的、简单的具有编辑、排版、处理各种表格与图片等功能的文字处理软件。本文就教学中会采用的一些Word的教学方法进行了简单的介绍与分析&#xff0c;并以泰山版初中计算机的教学计划为例&#xff0c;重…

poj2186Popular Cows(Kosaraju算法--有向图的强连通分量的分解)

1 /*2 题目大意&#xff1a;有N个cows, M个关系3 a->b 表示 a认为b popular&#xff1b;如果还有b->c, 那么就会有a->c 4 问最终有多少个cows被其他所有cows认为是popular&#xff01;5 6 思路&#xff1a;强连通分量中每两个节点都是可达的&#xff…

yum删除mysql数据库_MySQL数据库之Centos中彻底删除Mysql(rpm、yum安装的情况)

本文主要向大家介绍了MySQL数据库之Centos中彻底删除Mysql(rpm、yum安装的情况) &#xff0c;通过具体的内容向大家展现&#xff0c;希望对大家学习MySQL数据库有所帮助。我用的centos6&#xff0c;mysql让我整出了各种问题&#xff0c;我想重装一个全新的mysql&#xff0c;yum…

计算机视觉领域常见期刊和会议,计算机视觉领域常见期刊和会议

会议&#xff1a;ICCV&#xff1a; IEEE International Conference on Computer Vision(每两年举办一次&#xff0c;由IEEE主办&#xff0c;百度百科&#xff1a;https://baike.baidu.com/item/iccv/7054436?fraladdin&#xff0c;ICCV 2017&#xff1a;http://iccv2017.thecv…

mysql怎么新增_mysql怎么新增用户

匿名用户1级2018-01-27 回答展开全部首先以root身份登录到MySQL服务器中。$ mysql -u root -p当验证提示出现的时候&#xff0c;输入MySQL的root帐号的密码。创建一个MySQL用户使用如下命令创建一个用户名和密码分别为"myuser"和"mypassword"的用户。mysql…

江西财经大学计算机排名2019,2019年全国商科院校评价报告出炉 江西财经大学排名第七...

中国江西网/江西头条新闻客户端讯 记者郑周贇报道&#xff1a;6月5日&#xff0c;江西财经大学应用统计研究中心发布了2019年全国商科院校评价报告&#xff0c;江西财经大学在46所指标数据完整的商科院校中综合得分排名第七。该报告权衡了指标的完备性和可获得性&#xff0c;确…

Uvaoj10054 - The Necklace

1 /*2 题意&#xff1a;打印欧拉回路&#xff01;3 思路&#xff1a;开始时不明白&#xff0c;dfs为什么是后序遍历&#xff1f; 4 因为欧拉回路本身是一条回路&#xff0c;那么我们在dfs时&#xff0c;可能存在提前找到回路&#xff0c;这条回路可能不是欧拉回路&…

w7计算机应用放大按键,Win7窗口最大化和最小化快捷键是什么

Win7窗口最大化和最小化快捷键是什么Win7有很多快捷键&#xff0c;你都知道多少呢&#xff1f;今天小编给大家科普的是win7窗口最大化和最小化快捷键&#xff0c;下面就一起来了解看看吧&#xff01;Windows 键 方向键“↑”使当前使用的窗口最大化。Windows 键 方向键“↓”…

s7-1200跟mysql_让西门子S7-1200直接连接MySQL数据库!!!

最近项目上有个需求&#xff0c;要把采集的数据存储到数据库中&#xff0c;当前西门子有很多方法&#xff0c;必读IDB&#xff0c;还有通过WINCC的脚本&#xff0c;第三方的软件等等&#xff0c;但是随着发展&#xff0c;有些需求希望设备直接到数据库&#xff0c;比如云端的RD…

uva oj 567 - Risk(Floyd算法)

1 /*2 一张有20个顶点的图上。3 依次输入每个点与哪些点直接相连。4 并且多次询问两点间&#xff0c;最短需要经过几条路才能从一点到达另一点。5 6 bfs 水过7 */8 #include<iostream>9 #include<cstring> 10 #include<vector> 11 #include<cstdio> 12…

10034 - Freckles 克鲁斯克尔最小生成树!~

1 /*2 10034 - Freckles3 克鲁斯克尔最小生成树&#xff01;&#xff5e; 4 */5 #include<iostream>6 #include<cstdio>7 #include<cmath>8 #include<algorithm>9 using namespace std; 10 11 struct node{ 12 double x, y; 13 }; 14 15 struct t…

win7个人计算机的ip地址,win7计算机ip地址查询_win7本机ip地址查询

2016-12-09 11:40:21查找计算机的ip地址的方法&#xff1a;点击你的电脑桌面左下角的“开始”找到“运行”点击运行, 在出现的对话框里面输入“cmd” 点击确定然后就会出现一个黑色的命令行窗口,你会看到“&...2016-12-19 15:59:44手机设置静态IP 1、点设置-线网络 2、WLAN…

最全的mysql 5.7.13_最全的mysql 5.7.13 安装配置方法图文教程(linux) 强烈推荐!

linux环境Mysql 5.7.13安装教程分享给大家&#xff0c;供大家参考&#xff0c;具体内容如下1系统约定安装文件下载目录&#xff1a;/data/softwareMysql目录安装位置&#xff1a;/usr/local/mysql数据库保存位置&#xff1a;/data/mysql日志保存位置&#xff1a;/data/log/mysq…

iis 日志 post数据_云原生日志的趋势(1):logscape和logiq

作为日志产品的PM&#xff0c;跟进国内外日志产品动向是个长期工作。这几天翻新一些历史记录&#xff0c;发现logscape自2017年开源以来&#xff0c;突然2019年10月又更新了一会。于是顺着翻翻logscape的github账号&#xff0c;起了兴致来写点文字。https://github.com/logscap…