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 s s, her true speed is s+c s+c, where c c is an unknown constant (possibly negative).
Sheila made a careful record of a recent journey and wants to use this to compute c c. The journey consisted of n n segments. In the ith ith segment she traveled a distance of di di and the speedometer read si si for the entire segment. This whole journey took time t t. Help Sheila by computing c 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 n ( 1≤n≤1000 1≤n≤1000), the number of sections in Sheila’s journey, and t t ( 1≤t≤106 1≤t≤106), the total time. This is followed by n n lines, each describing one segment of Sheila’s journey. The ith ith of these lines contains two integers di di ( 1≤di≤1000 1≤di≤1000) and si si ( |si|≤1000 |si|≤1000), the distance and speedometer reading for the ith ith segment of the journey. Time is specified in hours, distance in miles, and speed in miles per hour.
Output
Display the constant c c in miles per hour. Your answer should have an absolute or relative error of less than 10−6 10−6.
| 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
``` | |
\
\
题意:
二分查找结果。
```cpp
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
int v[1100], d[1100], n;
double t;
int slove(double c){
double te = 0;
for(int i =0; i<n; i++){
if ((c + v[i]) <= 0)
return -1;
else
te += 1.0*d[i] /(c + v[i]);
}
if( te < t) // 返回1表示c取大了
return 1;
else return -1;
}
int main(){
while(scanf("%d%lf", &n,&t) != EOF){
for(int i =0; i < n; i++){
scanf("%d%d", &d[i], &v[i]);
}
double le = -1e9, ri = 1e9, mid = (le + ri) / 2;
int flag, cnt = 62;
while((ri - le) > 1e-6){
flag = slove(mid);
if(flag > 0){
ri = mid;
mid = (le + ri) /2;
}
else {
le = mid;
mid = (le + ri) /2;
}
}
printf("%0.6lf\n", mid);
}
return 0;
}
\
本文已参与「新人创作礼」活动,一起开启掘金创作之路