[BZOJ2712][[Violet 2]棒球][类欧几里得算法]

最新推荐文章于 2019-04-19 18:16:25 发布

原创最新推荐文章于 2019-04-19 18:16:25 发布 · 530 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#oi #bzoj #数论 #类欧几里得算法

Bzoj 同时被 3 个专栏收录

54 篇文章

订阅专栏

2017

50 篇文章

订阅专栏

数论

4 篇文章

订阅专栏

本文通过一道具体题目介绍了类欧几里得算法的应用，包括如何将小数转换为分数，并提供了一段完整的C++代码实现。适用于解决涉及分数比较的问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

[BZOJ2712][[Violet 2]棒球][类欧几里得算法]

类似于下面这道题吧，只要把小数转换成分数就好了。

http://blog.youkuaiyun.com/g1n0st/article/details/62044709

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<ll, ll> abcd;
inline ll Min(const ll &a, const ll &b) {
    return a < b ? a : b;
}
inline ll gcd(const ll &a, const ll &b) {
    if (!b) return a;
    return gcd(b, a % b);
}
inline abcd solve(ll p1, ll q1, ll p2, ll q2) {
    ll l = p1 / q1 + 1;
    abcd ret(0,0);
    if (l * q2 < p2) return abcd(l, 1);
    if (p1 == 0) return abcd(1, q2 / p2 + 1);
    if (p1 <= q1 && p2 <= q2) {
        ret = solve(q2, p2, q1, p1);
        return abcd(ret.second, ret.first); 
    }
    ll t = p1 / q1;
    ret = solve(p1 - q1 * t, q1, p2 - q2 * t, q2);
    ret.first += ret.second * t;
    return ret;
}
int main(void) {
    //freopen("in.txt", "r", stdin);
    ll p1, q1, p2, q2, r; int n;
    while (~scanf("%d 0.%lld", &n, &r)) {
        if (r == 0) {
            printf("1\n"); continue;
        }
        ll x = 10; while (n--) x *= 10;
        p1 = r * 10 - 5, q1 = x, p2 = r * 10 + 5, q2 = x;
        ll g = gcd(p1, q1); p1 /= g, q1 /= g;
        g = gcd(p2, q2); p2 /= g, q2 /= g;
        abcd ans = solve(p1, q1, p2, q2);
        printf("%lld\n", Min(ans.second, q1));
    }
    return 0;
}

完。

By g1n0st