【2017 World Final E】Need For Speed（二分）

Sheila is a student and she drives a typical student car: it is old, slow, rusty, and falling apart. Recently, the needle on the speedometer fell off. She glued it back on, but she might have placed it at the wrong angle. Thus, when the speedometer reads v, her true speed is v+c, where c is an unknown constant (possibly negative).

Sheila made a careful record of a recent journey and wants to use this to compute c. The journey consisted of n segments. In the i^th segment she traveled a distance of d_i and the speedometer read v_i for the entire segment. This whole journey took time t. Help Sheila by computing c.

Note that while Sheila’s speedometer might have negative readings, her true speed was greater than zero for each segment of the journey.

Input

The first line of input contains two integers n (1≤n≤1000), the number of sections in Sheila’s journey, and t (1≤t≤10⁶), the total time. This is followed by n lines, each describing one segment of Sheila’s journey. The i^th of these lines contains two integers d_i (1≤d_i≤1000) and v (|v_i|≤1000), the distance and speedometer reading for the i^th segment of the journey. Time is specified in hours, distance in miles, and speed in miles per hour.

Output

Display the constant c in miles per hour. Your answer should have an absolute or relative error of less than 10⁻⁶.

Sample Input 1	Sample Output 1
3 5 4 -1 4 0 10 3	3.000000000

Sample Input 2	Sample Output 2
4 10 5 3 2 2 3 6 3 1	-0.508653377

题意：

该旅程分为n个部分，每个部分的距离为d[i]，速度表所读取的速度为v[i] (v>0)，但真实速度为v[i]+c。已知旅程总时间为t，在总误差不超过1e-6的条件下，求c。

思路：

对于方程SUM(d[i]/s[i]+c)，利用二分求解的方法求c。

#include<bits/stdc++.h>

#define MAX 1050

using namespace std;

int d[MAX],v[MAX],n;

double t;

int ok(double c)

{

    double T=;

    for(int i=;i<n;i++)

    {

        if(c+v[i]<=)return ; //c偏小

        else T+=1.0*d[i]/(v[i]*1.0+c);

    }

    if(T>t)return ;

    else return ;             //c偏大

}

int main()

{

    int i;

    while(scanf("%d%lf",&n,&t)!=EOF)

    {

        for(i=;i<n;i++)

            scanf("%d%d",&d[i],&v[i]);

        double low=-1e9,high=1e9;

        while(high-low>1e-)

        {

            double mid=(high+low)*0.5;

            if(ok(mid))

                high=mid;

            else low=mid;

        }

        printf("%.9f\n",(high+low)*0.5);

    }

    return ;

}

巴特西

【2017 World Final E】Need For Speed（二分）

最新文章

热门文章