【bzoj 2433】【NOI 2011 D1T2】 智能车比赛 想法乱搞+暴力

新一届智能车大赛在JL大学开始啦！比赛赛道可以看作是由n个矩形区域拼接而成（如下图所示），每个矩形的边都平行于坐标轴，第i个矩形区域的左下角和右上角坐标分别为(xi,1,yi,1)和(xi,2,yi,2)。

选手们需要在最快的时间内让自己设计的智能车从一个给定的起点S点到达一个给定的终点T点，且智能车不能跑出赛道。假定智能车的速度恒为v且转向不消耗任何时间，你能算出最快需要多少时间完成比赛么？

Input
输入的第一行包含一个正整数n，表示组成赛道的矩形个数。

接下来n行描述这些矩形，其中第i行包含4个整数xi,1, yi,1, xi,2, yi,2，表示第i个矩形左下角和右上角坐标分别为(xi,1, yi,1)和(xi,2, yi,2)。

接下来一行包含两个整数xS, yS，表示起点坐标。

接下来一行包含两个整数xT, yT，表示终点坐标。

接下来一行包含一个实数v，表示智能车的速度。

Output
仅输出一个实数，至少精确到小数点后第六位，为智能车完成比赛的最快时间。

对于每个测试点，如果你的输出结果和参考结果相差不超过10^-6，该测试点得满分，否则不得分。

Sample Input
2
1 1 2 2
2 0 3 4
1 1
3 0
1.0

Sample Output
2.41421356

HINT
N<=2000,所输入数字为绝对值小于40000的整数

题目链接：

http://www.lydsy.com/JudgeOnline/problem.php?id=2433

思路：
1.给连续n个矩形，先O(N)与n+1比较，处理出两个矩阵的通道边界；如图中的虚线两端点，极限为2000*2个点；

 if(stx>x2[i]) continue;
        if(edx<x2[i]) break;

减去点；
2.n^2更新dis【i】，从起点到i的最短路。
！判断合法！i到j用斜率，在一定斜率范围内可行，if(tmp<=maxk&&tmp>=mink)

范围可更新；

  if(j%2==0)      mink=max(tmp,mink);

            else maxk=min(tmp,maxk);

代码

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<cmath>
#include<algorithm>
using namespace std;
int stx,sty,edx,edy;
int n;
int x1[2050],x2[2050],yy1[2050],y2[2050];
int num;
double v,dis[4050];
struct node{
    int x,y;
}p[5005];
int a[10];
void pre()
{
    num=1;
    p[1].x=stx;
    p[1].y=sty;
    for(int i=1;i<n;i++)
    {
        if(stx>x2[i]) continue;
        if(edx<x2[i]) break;
        a[1]=yy1[i],a[2]=y2[i],a[3]=yy1[i+1],a[4]=y2[i+1];
        sort(a+1,a+1+4);
        p[++num].x=x2[i];
        p[num].y=a[2];
        p[++num].x=x2[i];
        p[num].y=a[3];
    }
    p[++num].x=edx;
    p[num].y=edy;
}
void solve()
{
    for(int i=2;i<=num;i++)  
    dis[i]=1e99;
    for(int i=1;i<=num;i++)
    {
        double nox=(double)p[i].x*1.0;
        double noy=(double)p[i].y*1.0;
        double maxk=1e99;
        double mink=-1e99;
        for(int j=i+1;j<=num;j++)
        {
            double len;
            double nex=(double)p[j].x*1.0;
            double ney=(double)p[j].y*1.0;

            if(nex==nox)
            {
                dis[j]=min(dis[j],dis[i]+abs(ney-noy));
                continue;
            }
            double tmp=(double)(ney-noy)/(nex-nox);
            if(tmp<=maxk&&tmp>=mink)
            {
                double hh=(ney-noy)*(ney-noy)+(nex-nox)*(nex-nox);
                len=sqrt(hh);
                    dis[j]=min(dis[j],dis[i]+len);
            }
            if(j%2==0)
            {
                mink=max(tmp,mink);
            }   
            else maxk=min(tmp,maxk);
        }


    }
    double ans=dis[num]/v;
    printf("%.10lf",ans);
}
int main()

{
    scanf("%d",&n);
    for(int i=1;i<=n;i++) 
        scanf("%d%d%d%d",&x1[i],&yy1[i],&x2[i],&y2[i]);

    scanf("%d%d%d%d",&stx,&sty,&edx,&edy);

    scanf("%lf",&v);
    pre();
    solve();
}