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

最新推荐文章于 2019-10-09 15:53:30 发布

LowestJN

最新推荐文章于 2019-10-09 15:53:30 发布

阅读量1k

点赞数

分类专栏：类欧几里得算法文章标签：类欧几里得

本文链接：https://blog.csdn.net/coldef/article/details/62036060

版权

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

3 篇文章 0 订阅

订阅专栏

跟BZOJ2187类似
把小数转化成分数，就行了。

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cmath>

using namespace std;

typedef long long ll;
typedef pair<ll,ll> pairs;

int n;
ll r,x,a,b,c,d;

ll gcd(ll a,ll b){
    return b?gcd(b,a%b):a;
}

inline void simplify(ll &a,ll &b){
    ll g=gcd(a,b); a/=g; b/=g;
}

pairs solve(ll p1,ll q1,ll p2,ll q2){
    simplify(p1,q1); simplify(p2,q2);
    ll a=p1/q1+1,b=(ll)ceil((double)p2/q2)-1;
    if(a<=b) return pairs(a,1);
    if(p1==0) return pairs(1,q2/p2+1);
    if(p1<=q1&&p2<=q2){
        pairs r=solve(q2,p2,q1,p1);
        swap(r.first,r.second);
        return r;
    }
    ll t=p1/q1;
    pairs r=solve(p1%q1,q1,p2-t*q2,q2);
    r.first+=t*r.second;
    return r;
}

int main(){
    while(~scanf("%d 0.%lld",&n,&r)){
        if(r==0){ puts("1"); continue; }
        x=1; for(int i=1;i<=n+1;i++) x*=10;
        a=r*10-5; b=x; c=r*10+5; d=x;
        simplify(a,b); simplify(c,d);
        pairs Ans=solve(a,b,c,d);
        printf("%lld\n",min(Ans.second,b));
    }
}