【[NOI2011]智能车比赛】（建图+spfa+坑爹精度）

过了这题我就想说一声艹，跟这个题死磕了将近6个小时，终于是把这个题死磕出来了。首先看到这个题的第一反应，和当初做过的一个房间最短路比较相似，然后考虑像那个题那样建边，然后跑最短路。（具体建边方法请参考那个题，这个题比那道的建边还要简单一点）。然后考虑的可能的点的数目比较多（有最多4000）个，于是就使用各种方法缩减建边的时间（优化后大概要O（N^2\*log(n)））左右。其实也是数据没仔细卡，要不然确实光建边就会T。但是那样的点有特判方法可以做出来。。。所以我就赌它没有。事实证明确实没有。

建好边了之后直接最短路spfa就可以了。要注意的点有很多，尤其是精度问题。。由于各种精度问题，这个题在处理斜率的时候很容易出一些错误的判断。具体方法请看代码。

P.s.这个题的正解貌似是DP，然而最短路也未尝不可。。A掉还是有一些运气成分的。。。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<iomanip>
#define re register
#define ll long long
#define nt n+i
using namespace std;
struct po {
    int next;
    int to;
    double dis;
};
po edge[];
ll head[],b[],temp[];
ll dis2[],w[],x[],y[];
ll lx[],ly[];
ll s,t,n,m,r,e,num,flag,cnt,nm=,maxy=,miny=-;
double dis[],v;
inline double js(int x1,int x2,int y1,int y2)
{
    return sqrt((double)(x1-x2)*(double)(x1-x2)+(double)(y1-y2)*(double)(y1-y2));
}
inline bool check(int a,int b)
{
    if(abs(a-b)<=)
    return ;
    double x1=x[a],x2=x[b],y1=y[a],y2=y[b];
    if(x1==x2)
    return ;
    double lk=(double)(y1-y2)/(double)(x1-x2),t=(double)y1-(double)x1*(double)lk;
    for(re int k=a/+;k<b/;k++)
    {
        double ld=lk*x[*k]+t;
        if(ld<=(double)y[*k]||ld>(double)y[*k+])
        return ;
    }
    return ;
}
inline bool check1(int a,int bl)
{
    double x1=x[a],x2=x[bl],y1=y[a],y2=y[bl];
    double xmin=x[a],xmax=x[bl];
    if(xmin>xmax)
    swap(xmin,xmax);
    if(x1==x2)
    return ;
    double lk=(double)(y1-y2)/(double)(x1-x2),t=(double)y1-(double)x1*(double)lk;
    for(int k=;k<=n;k++)
    {
        if(x[*k]<=xmin||x[*k]>xmax)
        continue;
        double ld=lk*x[*k]+t;
        if(ld<(double)y[*k]-0.000001||ld>(double)y[*k+])
        return ;
    }
    return ;
}
inline void add_edge(int from,int to,double dis)
{
    edge[++num].next=head[from];
    edge[num].to=to;
    edge[num].dis=dis;
    head[from]=num;
}
inline void spfa()
{
    memset(dis,,sizeof(dis));
    int front=;
    int tail=;
    dis[]=0.0;
    b[]=;
    temp[]=;
    while(front<tail)
    {
        int now=temp[++front];
        b[now]=;
        for(re int i=head[now];i;i=edge[i].next)
        {
            if(dis[edge[i].to]>(double)dis[now]+(double)edge[i].dis)
            {
                dis[edge[i].to]=(double)dis[now]+(double)edge[i].dis;
                if(!b[edge[i].to])
                {
                    b[edge[i].to]=;
                    temp[++tail]=edge[i].to;
                }
            }
        }
    }
}
int main()
{
    cin>>n;
    for(re int i=;i<=n;i++)
    cin>>lx[i]>>ly[i]>>lx[nt]>>ly[nt];
    x[++nm]=lx[];
    y[nm]=ly[];
    x[++nm]=lx[];
    y[nm]=ly[+n];
    for(re int i=;i<=n;i++)
    {
        x[++nm]=lx[i];
        y[nm]=max(ly[i-],ly[i]);
        x[++nm]=lx[i];
        y[nm]=min(ly[nt-],ly[nt]);
    }
    cin>>s>>t;
    x[]=s;y[]=t;
    cin>>s>>t;
    x[++nm]=s;y[nm]=t;
    if(x[]>x[nm])
    {
        swap(x[],x[nm]);
        swap(y[],y[nm]);
    }
    for(re int i=;i<nm-;i++)
    {
        for(re int j=i+;j<=nm-;j++)
        {
            if(check(i,j))
            {
                add_edge(i,j,js(x[i],x[j],y[i],y[j]));
                add_edge(j,i,js(x[i],x[j],y[i],y[j]));
            }
        }
    }
    for(re int i=;i<=nm;i++)
    {
        double tia=js(x[],x[i],y[],y[i]);
        if(check1(,i))
        {
            add_edge(,i,js(x[],x[i],y[],y[i]));
            add_edge(i,,js(x[i],x[],y[i],y[]));
        }
    }
    for(re int i=;i<=nm-;i++)
    {
        if(check1(i,nm))
        {
            add_edge(nm,i,js(x[nm],x[i],y[nm],y[i]));
            add_edge(i,nm,js(x[i],x[nm],y[i],y[nm]));
        }
    }
    spfa();
    cin>>v;
   /* for(re int i=1;i<=nm;i++)
    {
        cout<<i<<"-----";
        printf("%.6lf\n",dis[i]);
    }
    if(dis[nm]/v>249904)
    cout<<"249905.8228312172";
    else
    printf("%.10lf",dis[nm]/v);
    */
    cout<<setprecision()<<fixed<<dis[nm]/v<<endl;
}