[FT][2]codevs 2190 有理逼近
来源:互联网 发布:java写界面 编辑:程序博客网 时间:2024/05/15 02:28
2190 有理逼近
时间限制: 1 s
空间限制: 32000 KB
题目等级 : 黄金 Gold
题解
查看运行结果
题目描述 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。
样例输入 Sample Input
样例1:
2 5
样例2:
5 100
样例输出 Sample Output
样例1:
4/3 3/2
样例2:
38/17 85/38
数据范围及提示 Data Size & Hint
P、N<30000
【分析】
枚举分母,二分分子…
【代码】
//有理逼近 #include<iostream>#include<cmath>#include<cstdio>#include<cstring>#include<algorithm>#define M(a) memset(a,0,sizeof a)#define fo(i,j,k) for(i=j;i<=k;i++)using namespace std;const int mxn=1e6;int n,p;double d,x;int a[mxn];inline int read(){ int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f;}inline int gcd(int x,int y) {if(x%y==0) return y;return gcd(y,x%y);}inline void yue(int &zi,int &mu){ int tmp=gcd(zi,mu); zi/=tmp,mu/=tmp;}int main(){ int i,j; p=read();n=read(); double d=sqrt(p); int x1=1,x2=1,y1=30000,y2=1,mu; fo(mu,1,n) //枚举分母 { int l=1,r=n; while(l<r) { int mid=(l+r)/2; double tmp=(double)mid/(double)mu; if(tmp>d) r=mid; else l=mid+1; } double tmp=(double)l/(double)mu; if(tmp<(double)y1/(double)y2 && tmp>d) y1=l,y2=mu; if(tmp>(double)x1/(double)x2 && tmp<d) x1=l,x2=mu; tmp=(double)(l-1)/(double)mu; if(tmp>(double)x1/(double)x2 && tmp<d) x1=l-1,x2=mu; if(tmp<(double)y1/(double)y2 && tmp>d) y1=l,y2=mu; yue(x1,x2);yue(y1,y2); } printf("%d/%d %d/%d\n",x1,x2,y1,y2); return 0;}
1 0
- [FT][2]codevs 2190 有理逼近
- codevs 2190 有理逼近 题解
- Tyvj P1077 CODE[VS] 2190 有理逼近
- TYVJ 1077 有理逼近
- codevs2190 有理逼近
- 【数学{枚举策略}】有理逼近
- 有理逼近的误差函数
- ft...
- 数值逼近学习笔记2
- [FT][2] 洛谷 P2365 任务安排
- [FT][2]洛谷 P1072 Hankson 的趣味题
- FT Introduction
- Codevs P1653 种树 2
- 线段覆盖2 CODEVS
- Codevs
- 【牛顿迭代逼近】求根号2的快速方法
- 逼近阶
- 连分数逼近。
- 在Qsys上运行程序时出现make elf文件失败的原因
- Android面试分享
- android与html5交互
- 马云:“电商”一词将被淘汰 新型娱乐购物崛起
- leetcode Longest Increasing Subsequence 最长递增子序列
- [FT][2]codevs 2190 有理逼近
- HTML 上机4
- ART运行时为新创建对象分配内存的过程分析
- bootstrap3-dialog 使用说明
- XCode 7.2 swift2.x Editor placeholder in source code
- android 状态栏颜色的两种设置模式
- java Excel 数据到导入在tomcat下正常,weblogic下异常原因
- HTML 第4章上机5
- openfalcon 咋就不发送邮件了来