bzoj4952 [Wf2017]Need for Speed

http://www.elijahqi.win/archives/634
Description
Sheila 是一名学生,她开着一辆经典的学生车:一辆又老,又慢,又锈,还老是崩坏的车。最近,时速表盘的指针还掉
了。她把指针粘了回去,但是她可能没有粘对角度。因此,当表盘读数为s时,她真实的速度可能是s+c,其中c为未知
常数 (可能是负的) 。Sheila 在最近的行程中仔细地做了一些记录,并希望能用这些记录来计算出c的值。行程由n
段组成。在第 i 段中,她匀速行驶了di的距离,表盘对应的读数一直为si。整个行程花费的时间为t。请你帮Sheila
确定c的值。注意即使Sheila的表盘可能有负的读数,她在每段行程的真实速度也是大于零的。
Input
第一行包含两个整数n(1≤n≤1000)和t(1≤t≤10^6),分别表示Sheila的行程段数和总时间。
接下来n行,每行描述了Sheila的一段行程。
第i行包含两个整数di(1≤di≤1000)和si(|si|≤1000),分别表示第 i 段行程的距离和表盘读数。
时间单位是小时,距离单位是英里,速度单位是英里每小时。
Output
输出常数c,其单位是英里每小时。你的答案绝对或相对误差应该小于10^-6。
Sample Input
样例1

3 5

4 -1

4 0

10 3

样例2

4 10

5 3

2 2

3 6

3 1
Sample Output
样例1

3.000000000

样例2

-0.508653377

思路：每次二分答案，验证是否正确

注意每个速度都必须>0 加入速度不>0 那么我们会不单调，不满足二分，可能会出现很多个答案

所以我提前找一下速度最小值，然后取反，作为左端点起始值，再找一下最大值，用可能的最大值去减去仪表盘的最大值

一开始初值给错wa，其实c有可能到1001000这么大

#include<cstdio>
#define N 1100
int d[N],c[N];
inline int read(){
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
inline int min(int x,int y){return x<y?x:y;}
inline int max(int x,int y){return x>y?x:y;}
int n,t;
inline bool check(double x){
    double ans1=0;
    for (int i=1;i<=n;++i){
        if (c[i]+x<0) return true;
        ans1+=(double)d[i]/(c[i]+x)*1.0;
    }
    if(ans1>=t) return true;else return false;
}
int main(){
    freopen("bzoj4952.in","r",stdin);
    n=read();t=read();int min1=0x7fffffff,max1=-0x7f7f7f7f;
    for (int i=1;i<=n;++i) d[i]=read(),c[i]=read(),min1=min(c[i],min1),max1=max(max1,c[i]);
    double l=-min1,r=1e6-max1;
    while (l<=r){
        double mid=(l+r)/2;
        if (check(mid))l=mid+1e-6;else r=mid-1e-6;
    }
    printf("%.9lf",l);
    return 0;
}