[HDU 6209] The Intersection

最新推荐文章于 2019-07-28 16:23:32 发布

zhysora

最新推荐文章于 2019-07-28 16:23:32 发布

阅读量662

点赞数

分类专栏：数学-数学原理文章标签：数学分数逼近

本文链接：https://blog.csdn.net/u013578420/article/details/78035107

版权

数学-数学原理专栏收录该内容

12 篇文章 0 订阅

订阅专栏

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6209

题目大意：T组数据，输入一个整数k，求最接近 $k^{\frac{2}{3}}$ 的分母不超过 $10^5$ 的有理数。( $T, k \leq 10^5$ )

思路：设 $x = k^{\frac{2}{3}}$ ，若x落在 $[\frac{a}{b}, \frac{c}{d}]$ 内，考虑新的分界点 $\frac{a + c}{b + d}$ ，每次都用新的边界更新答案，直到 $b + d > 10^5$ 。初始区间可以认为是 $[\lfloor x \rfloor, \lfloor x \rfloor + 1 ]$ 。

PS: 这个题要到处用long double

#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>

#define ld long double

using namespace std;

const int lim = (int)1e5;

int T, k;

int gcd(int a, int b){
    if (b == 0) return a;
    return gcd(b, a % b);
}
void simplify(int &a, int &b){
    int g = gcd(a, b);
    a /= g, b /= g;
}

int main(){
    scanf("%d", &T);
    while (T --){
        scanf("%d", &k);
        ld x = pow((ld)k, (ld)2.0 / 3);
        int a = (int)floor(x);
        int b = 1, c = a + 1, d = 1;

        if (1LL * a * a * a == 1LL * k * k){
            printf("%d/1\n", a); continue;
        }

        int u, v;
        if (fabs((ld)a - x) < fabs((ld)c - x)) u = a, v = 1;
        else u = c, v = 1;

        while(1){
            int i = a + c, j = b + d;
            simplify(i, j);

            if (j > lim) break;

            if (fabs((ld)1.0 * i / j - x) < fabs((ld)1.0 * u / v - x))
                u = i, v = j;

                if ((ld)1.0 * i / j < x) a = i, b = j;
            else c = i, d = j;
        }

        printf("%d/%d\n", u, v);
    }   

    return 0;
}