Codeforces Round #365 (Div. 2) C. Chris and Road 【数学】

如果按照题意思考，很难搞清楚情况。所以把问题转换一下，假设只有人在走，水平方向速度为v，竖直方向速度为u。根据物理知识，人在问题中的停止可以当场先在x轴上走一段距离，然后再同时朝上和朝右运动。又由于u和v已知，所以同时运动时的运动轨迹是确定的，而且是一条直线，假设这个斜率是a，又假设现在x轴上走了b米，那么运动方程是y=ax-ab。最后要求的就是一个最小的b满足y[i]>=a*x[i]-a*b，即所有点都在直线上方。另外特殊考虑b=0时是否所有点都在直线下方，如果是那么这个是最优解。

1

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
//#include <bits/stdc++.h>
#define rep(i,a,n) for(int i=a;i<n;i++)
#define per(i,a,n) for(int i=n-1;i>=a;i--)
#define pb push_back
using namespace std;
typedef vector<int> VI;
typedef long long ll;
const ll mod=1000000007;
const int N=1e5+10;
const double eps=1e-8;
double x[N],y[N],w,v,u,t,a,b,l,r,m,ans;
int n;

int get(double q)
{
    int f1,f2;
    f1=f2=0;
    for(int i=0;i<n;i++)
    {
        if(fabs(x[i]*a-q*a-y[i])<eps) continue;
        if(x[i]*a-q*a>y[i]) f1++;
        else f2++;
        if(f1&&f2) return 0;
    }
    return 1;
}

int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    //scanf("%d%lflf%lf",&n,&w,&v,&u);
    cin>>n>>w>>v>>u;
    r=0;
    for(int i=0;i<n;i++)
    {
        scanf("%lf%lf",&x[i],&y[i]);
        r=max(x[i],r);
    }
    t=w/u;
    a=w/(v*t);
    int f1,f2;
    f1=f2=0;
    for(int i=0;i<n;i++)
    {
        if(fabs(x[i]*a-y[i])<eps) continue;
        if(x[i]*a>y[i]) f1++;
        else f2++;
    }
    //printf("%d %d\n",f1,f2);
    if(f1!=0&&f2!=0)
    {
        m=0;
        for(int i=0;i<n;i++)
        {
            m=max(m,x[i]-y[i]/a);
        }
        ans=t+m/v;
        printf("%lf\n",ans);
    }
    else
    {
        printf("%lf\n",t);
    }
    return 0;
}

2

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
//#include <bits/stdc++.h>
#define rep(i,a,n) for(int i=a;i<n;i++)
#define per(i,a,n) for(int i=n-1;i>=a;i--)
#define pb push_back
using namespace std;
typedef vector<int> VI;
typedef long long ll;
const ll mod=1000000007;
const int N=1e5+10;
const double eps=1e-8;
double x[N],y[N],w,v,u,t,a,l,r,m,ans;
int n;

int get(double q)
{
    int f1,f2;
    f1=f2=0;
    for(int i=0;i<n;i++)
    {
        if(fabs(x[i]*a-q*a-y[i])<eps) continue;
        if(x[i]*a-q*a>y[i]) f1++;
        else f2++;
        if(f1&&f2) return 0;
    }
    return 1;
}

int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    //scanf("%d%lflf%lf",&n,&w,&v,&u);
    cin>>n>>w>>v>>u;
    r=0;
    for(int i=0;i<n;i++)
    {
        scanf("%lf%lf",&x[i],&y[i]);
        r=max(x[i],r);
    }
    t=w/u;
    a=w/(v*t);
    int f1,f2;
    f1=f2=0;
    for(int i=0;i<n;i++)
    {
        if(fabs(x[i]*a-y[i])<eps) continue;
        if(x[i]*a>y[i]) f1++;
        else f2++;
    }
    //printf("%d %d\n",f1,f2);
    if(f1!=0&&f2!=0)
    {
        l=0;
        while(fabs(r-l)>eps)
        {
            m=(l+r)/2;
            if(get(m)==1)
            {
                ans=m/v+t;
                r=m;
            }
            else
            {
                l=m;
            }
        }
        printf("%lf\n",ans);
    }
    else
    {
        printf("%lf\n",t);
    }
    return 0;
}