ACM刷题之HDU————相遇周期(真·打表(素数筛选))

相遇周期

Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 3136 Accepted Submission(s): 1465
 
Problem Description
2007年3月26日,在中俄两国元首的见证下,中国国家航天局局长孙来燕与俄罗斯联邦航天局局长别尔米诺夫共同签署了《中国国家航天局和俄罗斯联邦航天局关于联合探测火星-火卫一合作的协议》,确定中俄双方将于2009年联合对火星及其卫星“火卫一”进行探测。

而卫星是进行这些探测的重要工具,我们的问题是已知两颗卫星的运行周期,求它们的相遇周期。
 
Input
输入数据的第一行为一个正整数T, 表示测试数据的组数. 然后是T组测试数据. 每组测试数据包含两组正整数,用空格隔开。每组包含两个正整数,表示转n圈需要的天数(26501/6335,表示转26501圈要6335天),用\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\'/\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\'隔开。
 
Output

            对于每组测试数据, 输出它们的相遇周期,如果相遇周期是整数则用整数表示,否则用最简分数表示。
 
Sample Input
2
26501/6335 18468/42
29359/11479 15725/19170
 
Sample Output
81570078/7
5431415
 


关于分数的最小公倍数算法是:先通分,然后求出分子的最小公倍数,然后最小公倍数除以通分后的分母就是答案了(最后别忘化简)

没错这题的思路就是这么简单,但是分数的化简有点麻烦。 
下面是ac代码
#include<stdio.h>
__int64 p[20000]={0};
__int64 a[20000]={0};
int main()
{
	__int64 zu,zi1,zi2,mu1,mu2,zi,mu,m,i,c,num;
	 __int64 h,k=0;
	 //素数表 
	 for(i=1;i<20000;i++)
     a[i]=i+1;
   for(i=0;i<20000;i++)
     if(a[i]!=0)
       {
         p[k]=a[i];
         for(h=i;h<20000;h++)
           if(a[h]%p[k]==0)
            a[h]=0;
         k++;
       }
        
	scanf("%d",&zu);
	while(zu--)
	{
		scanf("%I64d/%I64d %I64d/%I64d",&zi1,&mu1,&zi2,&mu2);
		zi1=zi1*mu2;
		zi2=zi2*mu1;
		mu=mu1*mu2;		
		//求 公倍数
		i=zi1*zi2;
		if(zi1<zi2){
			c=zi1;
			zi1=zi2;
			zi2=c;
		}
		while(zi1%zi2!=0){
			c=(zi1%zi2);
			zi1=zi2;
			zi2=c;
		}
		zi=i/zi2;	
				
		if(zi%mu==0)
		{
			printf("%I64d\n",zi/mu);
			continue;
		}
		else
		{
			for(i=0;i<2000;)
			{
				
				if(zi%p[i]==0&&mu%p[i]==0)
				{
					zi/=p[i];
					mu/=p[i];
					continue;
				}
				i++;
			}
			printf("%I64d/%I64d\n",zi,mu);
		}
	}
}

化简分子分母,其实就是除以一个同时能除的数。这个数就是素数(其中的原理百度下吧。。我也是猜的)
所以要很多的素数一个个试过来才行。
这里就要用到打表了。
打表的
#include<stdio.h>
__int64 p[20000]={0};
__int64 a[20000]={0};
int main()
{
	 __int64 i;
	 __int64 h,k=0;
	 //素数表 
	 for(i=1;i<20000;i++)
     a[i]=i+1;
   for(i=0;i<20000;i++)
     if(a[i]!=0)   
         p[k]=a[i];
         for(h=i;h<20000;h++)
           if(a[h]%p[k]==0)
            a[h]=0;
         k++;
       }
}
一层一层去掉无用数据,就是方法。

  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值