POJ 3101 Astronomy

本文解析了POJ3101题目中的卫星同步算法,通过数学推导得出解决该问题的关键步骤,包括大整数乘法及素因子分解等,并提供了完整的C++实现代码。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

//poj 3101
//公式推导
//0.5*L/abs(L/ti-L/tj)--->相差半周的时间

//(ti*tj) / (2*abs(ti-tj))
//令 (ti*tj)=bi; (2*abs(ti-tj))=ai;

//相邻两个卫星平行的时间间隔di=bi/ai,问题转化为求这n-1个分数的最小公倍数


//分母p=gcd(a1,a2,...,an-1)
//分子q=lcm(b1,b2,...,bn-1)

//q/p 约分即得最终答案

//知识点:大整数乘法+素因子分解

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;

const int N=1000;
const int M=10000;

int n;
int t[N],r[N];

int c[M];

int gcd(int a,int b)
{
	if(!b) return a;
	return gcd(b,a % b);
}

void run()
{
	int i,j;

    for(i=0;i<n;i++) scanf("%d",&t[i]);

    r[0]=1;
    for(i=1;i<N;i++) r[i]=0;

    int a,b,g,d,k;

    d=0;
    for(i=1;i<n;i++) if(t[i]!=t[0])
    {
        b = t[i] * t[0];
        a = abs(t[i] - t[0])<<1;
        g = gcd(a,b);

        a/=g;
        b/=g;
        d=gcd(a,d);

		//分解素因子
        for(j = 2; b > 1; j++ ) if(b % j == 0)
        {
			k=0;
			while(b % j == 0)
			{
				b /= j;
				k++;
			}
			if(k > c[j]) c[j] = k;
        }
    }

	//lrj结论:  lcm(a,b)=p1^max(a1,b1)*p2^max(a2,b2)……pn^max(an,bn)
	//here:		预先计算出 所有bi的素因子个数 max 存放在 c[i] 数组中
	//大整数乘法 r=∏(i^c[i]) {c[i]!=0}
	int tmp;
	for(i = 0; i < M; i++)
	{
		for(j=0; j<c[i]; j++)
		{
			tmp = 0;
			for(k = 0; k < N; k++)
			{
				r[k] = r[k]*i+tmp;
				tmp = r[k] / 10000;
				r[k] %= 10000;
			}
		}
	}

	//ans : (r[i] / d)
	i=999;
	while(i && r[i] == 0) i--;

	printf("%d",r[i]);
	for(--i;i>=0;--i) printf("%04d",r[i]);

	printf(" %d\n",d);	
}

int main()
{
    while(scanf("%d",&n)!=EOF) run();
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值