最大比例
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 , 不能通过工程设置而省略常用头文件。
思路:我们假设这个最大的公比是q=(q1/q2) //q1,q2均为正整数且互质
我们将这个序列按照升序排列之后,a1,a2,a3…an,可以得出a2/a1=q1x1/q2x1, a3/a2=q1x2/q2x2…由于这个q一定是存在的,所以现在的问题就是分别找到分子和分母,
既然前提假设了q是确定的,那么q1,q2也一定是确定的,要想从 q1x1,q1x2,…q1xn中找到q1,我们可以利用类似于求gcd的辗转相除法,我们将所有的ai/ai-1进行约分确立分
子和分母,然后分别对所有的分子和分母进行辗转相除,最后约分输出就好了。还要注意的一点就是一定要对数组去重= =,不然会认为存在一个qx==1
注意点:
1.先对数据排序并去重
2.两两保存约分后的比值(把a1约掉)
3.求n个数的qgcd时,ans用第一个值赋初值,用后面的数更新就好,ans=qgcd(ans,a[i]);
4.qgcd模仿gcd算法
gcd算法:
a=k*b+r;
a和b的公因子c = k和r的公因子c
但是r和b是不断变小的,可以求出c
最后一定会出现整除的情况,所以r会等于0,此时gcd就是除数b。比如12,8。
综上
ll gcd(ll a,ll b)
{
if(b==0)
return a;
return gcd(b,a%b);
}
qgcd算法:
要想从 q1x1,q1x2,…q1xn中找到q1,就是q的1次方,q的2次方,和1q,2q的思路不一样。
我们发现,qn/qm=q(n-m),qn和qm的公因子q也就是qm和q(n-m)的公因子q
所以,
qn和qm的q = qm和q(n-m)的q
最后,一定会出现除数qn和商q(n-m)相等的情况,此时就是q。注意被除数要大,除数要小,b/a,所以每次a要选qm和q(n-m)中小的,b要选qm和q(n-m)中大的。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll a[110],q1[110],q2[110];
ll t1,t2;
ll gcd(ll a,ll b)
{
if(b==0)
return a;
return gcd(b,a%b);
}
ll qgcd(ll a,ll b)
{
if(a==b)
return a;
return qgcd(min(b/a,a),max(b/a,a));
}
int main()
{
int n;
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%lld",&a[i]);
sort(a+1,a+1+n);
int N = unique(a+1,a+1+n) - (a+1);
for(int i=1;i<=N-1;i++)
{
ll j = i+1;
int x=gcd(a[i],a[j]);
q1[++t1]=a[j]/x;
q2[++t2]=a[i]/x;
}
ll ans1=q1[1],ans2=q2[1];
for(int i=2;i<=t1;i++){
ans1=qgcd(ans1,q1[i]);
ans2=qgcd(ans2,q2[i]);
}
printf("%lld/%lld\n",ans1,ans2);
return 0;
}
几乎完全参考这个博客的思路,感谢大佬!
https://www.cnblogs.com/zzqc/p/8550086.html