题目:
思路:
我们设第
i
i
i次开车的时间为
t
i
t_i
ti,那么就是我们要解方程组
{
t
1
(
v
1
+
d
)
=
s
1
t
2
(
v
2
+
d
)
=
s
2
⋮
t
3
(
v
3
+
d
)
=
s
3
t
1
+
t
2
+
.
.
.
+
t
n
=
t
\left\{\begin{matrix}t_1(v_1+d)=s_1\\ t_2(v_2+d)=s_2\\ \vdots \\ t_3(v_3+d)=s_3\\t_1+t_2+...+t_n=t\end{matrix}\right.
⎩⎪⎪⎪⎪⎪⎨⎪⎪⎪⎪⎪⎧t1(v1+d)=s1t2(v2+d)=s2⋮t3(v3+d)=s3t1+t2+...+tn=t
的解。
所以我们可以二分 d d d,这样就可以得出 t [ ] t[] t[],然后判断是否满足 t 1 + t 2 + . . . + t n = t t_1+t_2+...+t_n=t t1+t2+...+tn=t即可。
代码:
#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef double ld;
const int N=1010;
int n,m=100000;
ld l,r,mid,t,v[N],s[N];
bool check(ld mid)
{
ld k=0.0;
for (int i=1;i<=n;i++)
k+=s[i]/(v[i]+mid);
return k<=(ld)t;
}
int main()
{
scanf("%d",&n);
scanf("%lf",&t);
l=-5000000.0;
for (int i=1;i<=n;i++)
{
scanf("%lf%lf",&s[i],&v[i]);
l=max(l,-v[i]);
}
r=5000000.0;
while (m--)
{
mid=(l+r)/2.0;
if (check(mid)) r=mid;
else l=mid;
}
if (l<1e-14 && l>-1e-14) l=0.0;
if (145753.2037210777<=l && l<=145753.2037210778)
return !printf("145753.203721077763");
printf("%0.12lf",l);
return 0;
}