2016年蓝桥杯省赛C/C++ A组 最大比例题解(100分)

10. 最大比例

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

输入格式:
第一行为数字N(n<100),表示接下的一行包含N个正整数
第二行N个正整数Xi(Xi<1 000 000 000 000),用空格分开。每个整数表示调查到的某人的奖金数额
要求输出:
一个形如A/B的分数,要求A、B互质。表示可能的最大比例系数
测试数据保证了输入格式正确,并且最大比例是存在的。

例如,输入:
3
1250 200 32
程序应该输出:
25/4
再例如,输入:
4
3125 32 32 200
程序应该输出:
5/2
再例如,输入:
3
549755813888 524288 2
程序应该输出:
4/1

资源约定:
峰值内存消耗 < 256M
CPU消耗 < 3000ms

请严格按要求输出,不要画蛇添足地打印类似:“请您输入…” 的多余内容。 所有代码放在同一个源文件中,调试通过后,拷贝提交该源码。
注意: main函数需要返回0 注意: 只使用ANSI C/ANSI C++ 标准,不要调用依赖于编译环境或操作系统的特殊函数。 注意:
所有依赖的函数必须明确地在源文件中 #include , 不能通过工程设置而省略常用头文件。
提交时,注意选择所期望的编译器类型。

这题比较复杂,我们需要分析一下。
首先比值都是相邻级别的,所以我们可以将这n个数排序再去重,
(如果只剩一个数的话,那答案应该是无穷大,这里就不考虑了)
然后把相邻的比值都算出来,为了方便讨论,再进行排序去重。
得到一个序列, f 0 = q k 0 , f 1 = q k 1 , … , f n − 1 = q k n − 1 f_0=q^{k_0}, f_1=q^{k_1}, \dots, f_{n-1}=q^{k_{n-1}} f0=qk0,f1=qk1,,fn1=qkn1 k i k_i ki 严格递增。

如果我们真的知道 q q q和这些整数 k k k,那么答案就是 q g q^g qg,其中 g g g是所有 k k k的最大公约数。
为什么呢?
因为我们要找的是一个 x x x,使得 log ⁡ q x f i \log_{q^x}f_i logqxfi 均为整数,即 k i x log ⁡ q q = k i x \frac{k_i}{x}\log_qq=\frac{k_i}{x} xkilogqq=xki均为整数,
显然 x x x是所有 k k k的约数,且 x x x最大,所以所求的 x x x就是 g g g

当然问题在于我们根本不知道 q q q和这些 k k k,但是我们很惊喜的看到了最大公约数。
这里我们可以借鉴辗转相减法的操作 辗转相减法资料
不妨取两个数,定义 Q ( a , b ) Q(a,b) Q(a,b)是由 a , b a,b a,b唯一确定的最大比例,
Q Q Q定义在无序对上,即 Q ( b , a ) = Q ( a , b ) Q(b,a)=Q(a,b) Q(b,a)=Q(a,b)
那么 Q ( a , b ) = Q ( q x , q y ) = q g c d ( x , y ) = Q ( q x , q y − x ) = Q ( a , b a ) a < b Q(a,b)=Q(q^x,q^y)=q^{gcd(x,y)}=Q(q^x,q^{y-x})=Q(a,\frac{b}{a})\quad a<b Q(a,b)=Q(qx,qy)=qgcd(x,y)=Q(qx,qyx)=Q(a,ab)a<b
所以我们只需要重复这个“辗转相减法”就可以得到最后的答案。

具体实现的时候可以每次只操作相邻的两个比值。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

const int N=110;
const ll MAX=1000000000000;
ll a[N];

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

struct Frac
{
    ll up,dw;
    Frac() {}
    Frac(ll a,ll b) {
        if (a==0) {
            up=0;dw=1;
            return;
        }
        ll g=gcd(a,b);
        up=a/g;dw=b/g;
    }
    bool operator <(const Frac &b) const {
        return up*b.dw<dw*b.up;
    }
    bool operator ==(const Frac &b) const {
        return !((*this)<b)&&!(b<(*this));
    }
    Frac operator /(const Frac &b) {
        return Frac(up*b.dw,dw*b.up);
    }
    void print() {
        printf("%I64d/%I64d\n",up,dw);
    }
} f[N];

int main()
{
    int n;
    scanf("%d",&n);
    for (int i=0;i<n;i++) {
        scanf("%I64d",&a[i]);
    }
    sort(a,a+n);
    n=unique(a,a+n)-a-1;
    for (int i=0;i<n;i++) {
        f[i]=Frac(a[i+1],a[i]);
    }
    Frac ans(MAX,1);
    while (n>1) {
        sort(f,f+n);
        if (f[0]<ans) ans=f[0];
        n=unique(f,f+n)-f-1;
        for (int i=0;i<n;i++) {
            f[i]=f[i+1]/f[i];
        }
    }
    if (f[0]<ans) ans=f[0];
    ans.print();
    return 0;
}

评测记录截图

  • 10
    点赞
  • 14
    收藏
    觉得还不错? 一键收藏
  • 5
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值