本文主要是介绍codevs2190 有理逼近,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
题目描述 Description
对于一个素数P,我们可以用一系列有理分数(分子、分母都是不大于N的自然数)来逼近sqrt(p),例如P=2,N=5的时候:1/1<5/4<4/3< sqrt(2)<3/2<5/3<2/1。
任 务 :
给定P、N(N>sqrt(p)),求X、Y、U、V,使x/y< sqrt(p)< u/v且x/y与sqrt(p)之间、sqrt(p)与u/v之间都不能再插入满足题意的有理分数。
输入描述 Input Description输入文件的第一行为P、N 输出描述 Output Description
输出文件只有一行,格式为“X/Y U/V”。注意,答案必须是既约的,也就是说分子、分母的最大公约数必须等于1。
枚举分母,对于这个分母找到最接近的分子,二分查找。
#include<cstdio>
#include<cstring>
#define LL long long
int gcd(int x,int y)
{return y?gcd(y,x%y):x;
}
int main()
{int i,j,k,m,n,p,q,x,y,z,l,r,mid,la=-1,lb,ua=-1,ub;scanf("%d%d",&p,&n);for (i=1;i<=n;i++){l=1;r=n;while (l<=r){mid=(l+r)/2;if (mid*mid>(LL)p*i*i){if (gcd(i,mid)==1&&(ua==-1||mid*ua<i*ub)){ua=i;ub=mid;}r=mid-1;}else{if (gcd(i,mid)==1&&(la==-1||mid*la>i*lb)){la=i;lb=mid;}l=mid+1;}}}printf("%d/%d %d/%d\n",lb,la,ub,ua);
}
这篇关于codevs2190 有理逼近的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!