TYVJ P1077 有理逼近 Label：坑，tle的好帮手 不懂

描述

对于一个素数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之间都不能再插入满足题意的有理分数。

输入格式

输入文件的第一行为P、N，其中 P、N<30000。

输出格式

输出文件只有一行，格式为“X/Y U/V”。注意，答案必须是既约的，也就是说分子、分母的最大公约数必须等于1。

测试样例1

输入

样例1： 
2 5 
样例2： 
5 100

输出

样例1： 
4/3 3/2

样例2： 
38/17 85/38

代码

 #include<iostream>
 #include<cstdio>
 #include<vector>
 #include<cmath>
 #include<queue>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 int des_int,P,N,a,b,ans[]={,,,,,,,,,};
 double des;

 int check(int i,int j,double x){
      return (fabs( des-double(ans[i])/double(ans[j]) ) > fabs( des-x ));
 }

 void xiao(){//ans1~2
     for(int i=des_int-;i<=N;i++){
         int flag=;
         for(int j=;j<=i;j++){
             double x=double(i)/double(j);
             if(x>des) continue;
             if(check(,,x)){
                 double k1=(double(ans[])/double(ans[])),k2=x;
                 bool fl=(k1==k2);
                 if(fl) continue;
                 ans[]=i;ans[]=j;
             }
             else{
                 flag=;break;
             }
         }
         if(flag==) continue;
     }
 }

 void da(){//ans3~4
     for(int i=des_int+;i<=N;i++){
         int flag=;
         for(int j=i-;j>=;j--){
             double x=double(i)/double(j);
             if(x<des) continue;
             if(check(,,x)){
                 if(fabs( des-double(ans[])/double(ans[]) )== fabs( des-x )) continue;
                 ans[]=i;ans[]=j;
             }
             else{
                 flag=;break;
             }
         }
         if(flag==) continue;
     }
 }
 int main(){
 //    freopen("01.txt","r",stdin);
     scanf("%d %d",&P,&N);
     des=sqrt(double(P));
     des_int=des;
     xiao();
     da();

     printf("%d/%d %d/%d\n",ans[],ans[],ans[],ans[]);
     return ;
 }

穷尽今生所学，T掉2个点

上看不懂的别人代码：

 #include <math.h>
 #include <stdio.h>
 #include <stdlib.h>

 int main(int argc, char **argv)
 {
     int i, j;
     int p, n;
     int a, b, c, d;
     int beg, end;
     double s, max = , min = , t;
     scanf("%d%d", &p, &n);
     s = sqrt(p);
     for(i = ; i <= n; i++){
         beg = i * s - ;
         end = beg + ;
         if(beg <= ){
             beg = ;
         }
         if(end > n){
             end = n;
         }
         for(j = beg; j <= end; j++){
             t = (double)(j) / i;
             if(t < s && t > max){
                 a = j, b = i;
                 max = t;
             }
         }
     }
     for(i = ; i <= n; i++){
         beg = i * s;
         end = beg + ;
         if(beg <= ){
             beg = ;
         }
         if(end > n){
             end = n;
         }
         for(j = beg; j <= end; j++){
             t = (double)(j) / i;
             if(t > s && min > t){
                 c = j, d = i;
                 min = t;
             }
         }
     }
     printf("%d/%d %d/%d\n", a, b, c, d);
     return ;
 }

---------------我是猥琐的分割线---------------

 #include<iostream>
 #include<cstdio>
 #include<vector>
 #include<cmath>
 #include<queue>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 int des_int,P,N,a,b,ans[]={,,,,,,,,,};
 double des;

 int check(int i,int j,double x){
      return (fabs( des-double(ans[i])/double(ans[j]) ) > fabs( des-x ));
 }

 void xiao(){//ans1~2
     for(int i=des_int-;i<=N;i++){
         int flag=;
         for(int j=;j<=i;j++){
             double x=double(i)/double(j);
             if(x>des) continue;
             if(check(,,x)){
                 double k1=(double(ans[])/double(ans[])),k2=x;
                 bool fl=(k1==k2);
                 if(fl) continue;
                 ans[]=i;ans[]=j;
             }
             else{
                 flag=;break;
             }
         }
         if(flag==) continue;
     }
 }

 void da(){//ans3~4
     for(int i=des_int+;i<=N;i++){
         int flag=;
         for(int j=des_int/i+;j>=;j--){
             double x=double(i)/double(j);
             if(x<des) continue;
             if(check(,,x)){
                 if(fabs( des-double(ans[])/double(ans[]) )== fabs( des-x )) continue;
                 ans[]=i;ans[]=j;
             }
             else{
                 flag=;break;
             }
         }
         if(flag==) continue;
     }
 }
 int main(){
 //    freopen("01.txt","r",stdin);
     scanf("%d %d",&P,&N);
     des=sqrt(double(P));
     des_int=des;
     xiao();
     da();

     printf("%d/%d %d/%d\n",ans[],ans[],ans[],ans[]);
     return ;
 }

39行调了个阈值成功骗分，主要都是da（）函数耗的时间，也就没多想，+500有点危险，+100绝对过不了

想法来源于（能跑这么久的数肯定很大，分母也就没必要那么大了！！！）