bzoj 2187: fraction 类欧几里德算法

最新推荐文章于 2020-11-06 15:44:17 发布

_beginend

最新推荐文章于 2020-11-06 15:44:17 发布

阅读量573

点赞数

分类专栏：类欧几里德算法

本文链接：https://blog.csdn.net/qq_33229466/article/details/76975258

版权

类欧几里德算法专栏收录该内容

2 篇文章 0 订阅

订阅专栏

题意

分数给你4个正整数a，b，c，d，求一个最简分数 p / q满足 a / b < p / q < c / d，若有多组解，输出q最小的一组，若仍有多组解，输出p最小的一组。
1 < = a，b，c，d < = 10^9

分析

又是一道经典的类欧题目。推一波做法：

若 a b 与 c d 之 间 存 在 一 个 整 数 ， 则 返 回 这 个 整 数 即 可 。

$若\frac{a}{b}与\frac{c}{d}之间存在一个整数，则返回这个整数即可。$

若 a = 0 ， 那 么 p q < c d = > q > p d c

$若a=0，那么\frac{p}{q}<\frac{c}{d}=>q>\frac{pd}{c}$

由 于 q 要 最 小 ， 所 以 p = 1, q = ⌊ d c ⌋ + 1

$由于q要最小，所以p=1,q=\lfloor\frac{d}{c}\rfloor+1$

若 a < = b 且 c < = d 则 d c < q p < b a

$若a<=b且c<=d则\frac{d}{c}<\frac{q}{p}<\frac{b}{a}$

否 则 的 话 ， a % b b < p q - ⌊ a b ⌋ < c d - ⌊ a b ⌋

$否则的话，\frac{a\%b}{b}<\frac{p}{q}-\lfloor\frac{a}{b}\rfloor<\frac{c}{d}-\lfloor\frac{a}{b}\rfloor$

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long LL;

LL gcd(LL a,LL b)
{
    if (!b) return a;
    else return gcd(b,a%b);
}

void sim(LL &a,LL &b)
{
    LL d=gcd(a,b);a/=d;b/=d;
}

void solve(LL a,LL b,LL c,LL d,LL &p,LL &q)
{
    sim(a,b);sim(c,d);
    LL x=a/b+1,y=c/d+(c%d>0)-1;
    if (x<=y) p=x,q=1;
    else if (!a) p=1,q=d/c+1;
    else if (a<=b&&c<=d) solve(d,c,b,a,p,q),swap(p,q);
    else solve(a%b,b,c-d*(a/b),d,p,q),p+=q*(a/b);
}

int main()
{
    LL a,b,c,d;
    while (scanf("%lld%lld%lld%lld",&a,&b,&c,&d)!=EOF)
    {
        LL p,q;
        solve(a,b,c,d,p,q);
        printf("%lld/%lld\n",p,q);
    }
}