最大比例(数论,简化分数)

最大比例

题目链接
X星球的某个大奖赛设了 M 级奖励。
每个级别的奖金是一个正整数。
并且，相邻的两个级别间的比例是个固定值。
也就是说：所有级别的奖金数构成了一个等比数列。
比如：16,24,36,54，其等比值为：3/2。
现在，我们随机调查了一些获奖者的奖金数。
请你据此推算可能的最大的等比值。

输入格式

第一行为数字 N ，表示接下的一行包含 N 个正整数。
第二行 N 个正整数 Xi，用空格分开，每个整数表示调查到的某人的奖金数额。

输出格式

一个形如 A/B 的分数，要求 A、B 互质，表示可能的最大比例系数。

数据范围

0<N<100
0<Xi<10¹²
数据保证一定有解。

输入样例1：

3
1250 200 32

输出样例1：

25/4

输入样例2：

4
3125 32 32 200

输出样例2：

5/2

输入样例3：

3
549755813888 524288 2

输出样例3：

4/1

算法分析

输入数可以转化为a,a*(p/q)¹,a*(p/q)²,…,a*(p/q)^(n-1)
我们一后面每个数都除以第一个数就可以得到
(p/q)¹,(p/q)²,…,(p/q)^(n-1)这样的一些数,我们分别存储分子和分母
存储的时候先进行第一次简化,除以最大公约数
分子就是原数组 a /最大公约数 分母就是原数组a*(p/q)^(n-1)/最大公约数.
剩下的是p/q的次方的一些数组.比如说(p/q)¹,(p/q)³,(p/q)⁷,我们得出来满足最少公比就是gcd((p/q)^3-1,(p/q)^7-3)==(p/q)²,所以我们就是要求指数的最大公约数.这一项求值我们采用更相减损法的变形来计算.指数相减就是互相相除.

代码实现

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=1e3+5;
typedef long long ll;
int n;
ll x[maxn];
ll a[maxn];
ll b[maxn];
ll gcd(ll a,ll b)
{
    return b?gcd(b,a%b):a;
}
ll jgcd(ll a,ll b)
{
    if(a<b)
        swap(a,b);
    if(b==1)
        return a;
    return jgcd(b,a/b);
}
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++)   cin>>x[i];
    sort(x+1,x+n+1);
    int cnt=0;
    for(int i=2;i<=n;i++)
    {
        if(x[i]!=x[i-1])
        {
            ll d=gcd(x[i],x[1]);
            a[cnt]=x[i]/d;
            b[cnt++]=x[1]/d;
        }
    }
    ll up=a[0],down=b[0];
    for(int i=1;i<cnt;i++)
    {
        up=jgcd(up,a[i]);
        down=jgcd(down,b[i]);
    }
    cout<<up<<'/'<<down<<endl;
    return 0;
}