luogu P3743 kotori的设备（二分小数，小数高精度）

题目链接：https://www.luogu.com.cn/problem/P3743

题意：输入n,p（n为设备个数，p为每分钟可以填补的能量）

n组a,b（a为每分钟消耗的能量，b为最初始的能量）

注意，可以再任意时刻给任意一个设备补能。问在一个设备能量为0之前最多可以使用多少分钟（保留9位小数），如果能无限使用，输出-1。

题解：洛谷二分答案专题做的，其实知道是二分答案的话会好做很多。就二分时间，然后看是否所有的设备都还没到0.注意精度问题（要是再codeforces做这题就凉了，eps最初0.0000000001,r=1e9,超时。发现样例测试数据t有超过r的，那就r=1e10，超时再eps=0.000001才全过）。

总结呢：这种能量填补题，或者设计到需要用到eps的，可能会是要用到二分答案。

代码：

#include <bits/stdc++.h>

#define ll long long
#define pi acos(-1)
#define pb push_back
#define mst(a, i) memset(a, i, sizeof(a))
#define pll pair<ll, ll>
#define fi first
#define se second
#define dbg(x) cout << #x << "===" << x << endl
#define dbgg(i, x) cout << #x << "[" << i << "]===" << x[i] << endl
using namespace std;
template<class T>void read(T &x){T res=0,f=1;char c=getchar();while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}while(isdigit(c)){res=(res<<3)+(res<<1)+c-'0';c=getchar();}x=res*f;}
void print(ll x){if(x<0){putchar('-');x=-x;}if(x>9)print(x/10);putchar(x%10+'0');}
const ll maxn = 1e5+ 10;
const ll mod = 1e9 + 7;
const double eps = 0.000001;

ll n;
double p;
struct node{
    double a,b;
}s[maxn];
bool judge(double x){
    double t=0;
    for(ll i=1;i<=n;i++){
        if(s[i].a*x-s[i].b<=0) continue;
        t+=(s[i].a*x-s[i].b)/p;
        //dbg(t);
    }
    if(t<=x) return true;
    else return  false;
}
ll gcd(ll a,ll b){return (b==0)?a:gcd(b,a%b);}
ll qpow(ll a,ll p,ll mod){ll ans=1;a=a%mod;while(p){if(p&1)ans=(ans*a)%mod;p>>=1;a=(a*a)%mod;}return ans;}
int main() {
    ll _s = 1;
    //read(_s);
    for (ll _=1;_<=_s;_++) {
        cin>>n>>p;
        ll cnt=0;
        double sum=0;
        for(ll i=1;i<=n;i++){
            cin>>s[i].a>>s[i].b;
            sum+=s[i].a;
        } 
        if(p>=sum){
            cout<<-1<<endl;
            break;
        }

        double l=0,r=1e10;
        //printf("%.9lf\n",r);
        while((r-l)>eps){
            double mid=(l+r)/2;
            if(judge(mid)) l=mid;
            else r=mid-eps;
        }
       printf("%.9lf",l);
    }
    return 0;
}