问题描述
在日常生活中,我们在收集数据之后,习惯用比例描述两数之间的比例关系,但是如果数据量过大,如566487:3842145,我们很难看出两者之间的实际关系,但是如果把数据范围压缩,变成3:20,,虽然数值上造成了一定损失,但是依旧能看出两者之间的关系,并更容易理解。
所以有以下的算法问题:对于任意的输入n:m,在最大范围不超过l下,我们要算出在l的范围内,最接近原始比值的数n':m'。其中n'和m'为最简比例,且互质。
输入格式
n m l
输入描述
输入有若干行,1<=n,m<=2147483647;1<=l<=n,m;1<=l<=50000。
输出格式
n' m'
输出描述
对于任意的输入n,m,均要算出对应的n',m',按行显示,其中n':m'>=n:m;
样例输入
566487 3842145 20
1415926 8584074 100
样例输出
3 20
16 97
实现代码
我的代码如下,利用了区域优化的思想,即先把输入的比例压缩为在l下的最大比例。因为任意一个比这个设置的数更大的,会让结果离真实结果更远
另外,利用状态的判断,算法可以更清晰地知道,下一步应该给左边的数自减,还是给右边的数自减,可以把时间复杂度从O(n^2)压缩到O(l+min(n,m)/max(n,m)*l)
#include<iostream>
#include<cmath>
using namespace std;
int main()
{
double n,m,l,min,max,result,temp;
int left,right;
while(cin>>n>>m>>l)
{
min=left=right=l;
result=n/m;
max=n>m?n:m;
n=ceil(n/max*l);
m=ceil(m/max*l);
//以上三行代码,就是区域优化的实现
while(n&&m)
{
temp=n/m-result;
if(temp>=0&&temp<=min)
{
min=temp;
left=n;
right=m;
}
if(n/m>=result)n--;else m--;
}
cout<<left<<" "<<right<<endl;
}
return 0;
}
总结
该算法的代码精简,当l设置为10000000时,依旧可以秒出结果