POJ - 1066 Treasure Hunt

传送门：Treasure Hunt

题意

有一个左下角为（0,0），右上角为（100,100）的正方形，给出n条线段，线段的起点和终点都在正方形上，给定一个宝藏的坐标，问从正方形外走到宝藏至少要经过多少条线段。每次只能走线段的中点，也就是每两个交点的中点。 

题解

因为肯定是从正方形的边开始走，所以枚举每条边上的相邻两个点的中点，连接宝藏所在的点，看这条线段和几条线段有交点，并更新最小值。（为什么这么做是对的呢？因为我没找到反例（逃

存每条边上的点时记得排序 : ) ，忘记排序花式wa。

代码

#include<stdio.h>
#include<algorithm>
#include<iostream>
#include<string.h>
#include<math.h>
#define eps 1e-6
using namespace std;

struct point
{
    double x,y;
};

struct line
{
    point s,t;
}l[1000];

double a[10010];

double cross(point a,point b)
{
    return a.x*b.y-a.y*b.x;
}

point operator - (point a,point b)
{
    return {a.x-b.x,a.y-b.y};
}

bool intersection(point a,point b,point c,point d)
{
    //快速排斥实验
    if(max(c.x,d.x)<min(a.x,b.x)||max(a.x,b.x)<min(c.x,d.x)||max(c.y,d.y)<min(a.y,b.y)||max(a.y,b.y)<min(c.y,d.y)){
        return false;
    }
    //跨立实验
    if(cross(a-d,c-d)*cross(b-d,c-d)>0||cross(d-b,a-b)*cross(c-b,a-b)>0){
        return  false;
    }
    return true;
}

int dcmp(double x)
{
    return fabs(x)<eps?0:x<0?-1:1;
}

int main()
{
    int n;
    scanf("%d",&n);
    for(int i=0;i<n;i++){
        scanf("%lf%lf%lf%lf",&l[i].s.x,&l[i].s.y,&l[i].t.x,&l[i].t.y);
    }
    l[n++]={{0,0},{0,100}};
    l[n++]={{0,0},{100,0}};
    l[n++]={{100,100},{0,100}};
    l[n++]={{100,100},{100,0}};
    point q;
    scanf("%lf%lf",&q.x,&q.y);
    line p;
    int ans=0x3f3f3f3f;
    int k=0;
    for(int i=0;i<n;i++){
        if(dcmp(l[i].s.x)==0){
            a[k++]=l[i].s.y;
        }
        if(dcmp(l[i].t.x)==0){
            a[k++]=l[i].t.y;
        }
    }
    sort(a,a+k);
    for(int i=1;i<k;i++){
        p.s={0,(a[i]+a[i-1])/2.0};
        p.t=q;
        int cnt=0;
        for(int j=0;j<n;j++){
            if(intersection(p.s,p.t,l[j].s,l[j].t)) cnt++;
        }
        ans=min(ans,cnt);
    }
    k=0;
    for(int i=0;i<n;i++){
        if(dcmp(l[i].s.y)==0){
            a[k++]=l[i].s.x;
        }
        if(dcmp(l[i].t.y)==0){
            a[k++]=l[i].t.x;
        }
    }
    sort(a,a+k);
    for(int i=1;i<k;i++){
        p.s={(a[i]+a[i-1])/2.0,0.0};
        p.t=q;
        int cnt=0;
        for(int j=0;j<n;j++){
            if(intersection(p.s,p.t,l[j].s,l[j].t)) cnt++;
        }
        ans=min(ans,cnt);
    }
    k=0;
    for(int i=0;i<n;i++){
        if(dcmp(l[i].s.x-100)==0){
            a[k++]=l[i].s.y;
        }
        if(dcmp(l[i].t.x-100)==0){
            a[k++]=l[i].t.y;
        }
    }
    sort(a,a+k);
    for(int i=1;i<k;i++){
        p.s={100,(a[i]+a[i-1])/2.0};
        p.t=q;
        int cnt=0;
        for(int j=0;j<n;j++){
            if(intersection(p.s,p.t,l[j].s,l[j].t)) cnt++;
        }
        ans=min(ans,cnt);
    }

    k=0;
    for(int i=0;i<n;i++){
        if(dcmp(l[i].s.y-100)==0){
            a[k++]=l[i].s.x;
        }
        if(dcmp(l[i].t.y-100)==0){
            a[k++]=l[i].t.x;
        }
    }
    sort(a,a+k);
    for(int i=1;i<k;i++){
        p.s={(a[i]+a[i-1])/2.0,100};
        p.t=q;
        int cnt=0;
        for(int j=0;j<n;j++){
            if(intersection(p.s,p.t,l[j].s,l[j].t)) cnt++;
        }
        ans=min(ans,cnt);
    }
    if(q.x<0||q.x>100||q.y<0||q.y>100) ans=0;
    printf("Number of doors = %d\n",ans);
    return 0;
}