题目: 传送门

题意:在一个左小角坐标为(0, 0),右上角坐标为(10, 10)的房间里,有 n 堵墙,每堵墙都有两个门。每堵墙的输入方式为 x, y1, y2, y3, y4,x 是墙的横坐标,第一个门的区间为[ (x, y1) ~ (x, y2) ],问你从 (0, 5) 走到 (10, 5) 的最短路径是多少。

0 <= n <= 18

题解:讲每个门的端点存起来,然后,加上(0, 5) 和 (10, 5) 这两个点,暴力判断这些点两两间是否能互相直达,然后再跑一遍最短路就行了。

      你不能直接从 (0, 5) 到 (10, 5) 那你肯定是经过一些门的端点再到终点是最优的。

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <queue>
#include <map>
#include <vector>
#include <set>
#include <string>
#include <math.h>
#define LL long long
#define mem(i, j) memset(i, j, sizeof(i))
#define rep(i, j, k) for(int i = j; i <= k; i++)
#define dep(i, j, k) for(int i = k; i >= j; i--)
#define pb push_back
#define make make_pair
#define INF 1e20
#define inf LLONG_MAX
#define PI acos(-1)
using namespace std; const int N = 1e2 + ;
const double eps = 1e-; struct Point {
double x, y;
Point(double x = , double y = ) : x(x), y(y) { } /// 构造函数
}; /// 向量加减乘除
inline Point operator + (const Point& A, const Point& B) { return Point(A.x + B.x, A.y + B.y); }
inline Point operator - (const Point& A, const Point& B) { return Point(A.x - B.x, A.y - B.y); }
inline Point operator * (const Point& A, const double& p) { return Point(A.x * p, A.y * p); }
inline Point operator / (const Point& A, const double& p) { return Point(A.x / p, A.y / p); } inline int dcmp(const double& x) { if(fabs(x) < eps) return ; else return x < ? - : ; } inline double Cross(const Point& A, const Point& B) { return A.x * B.y - A.y * B.x; } /// 叉积
inline double Dot(const Point& A, const Point& B) { return A.x * B.x + A.y * B.y; } /// 点积
inline double Length(const Point& A) { return sqrt(Dot(A, A)); } /// 向量长度
inline double Angle(const Point& A, const Point& B) { return acos(Dot(A, B) / Length(A) / Length(B)); } /// 向量A,B夹角 inline Point GetLineIntersection(const Point P, const Point v, const Point Q, const Point w) {///求直线p + v*t 和 Q + w*t 的交点,需确保有交点,v和w是方向向量
Point u = P - Q;
double t = Cross(w, u) / Cross(v, w);
return P + v * t;
} inline bool Onsegment(Point p, Point a1, Point a2) { /// 判断点p是否在线段p1p2上
return dcmp(Cross(a1 - p, a2 - p)) == && dcmp(Dot(a1 - p, a2 - p)) <= ;
} inline bool SegmentProperInsection(Point a1, Point a2, Point b1, Point b2) { /// 判断线段是否相交
if(dcmp(Cross(a1 - a2, b1 - b2)) == ) /// 两线段平行
return Onsegment(b1, a1, a2) || Onsegment(b2, a1, a2) || Onsegment(a1, b1, b2) || Onsegment(a2, b1, b2);
Point tmp = GetLineIntersection(a1, a2 - a1, b1, b2 - b1);
return Onsegment(tmp, a1, a2) && Onsegment(tmp, b1, b2);
} Point P[N];
double dis[N][N];
int main() {
int n;
double x, y1, y2, y3, y4;
while(scanf("%d", &n) && n != -) { int cnt = ;
rep(i, , n) {
scanf("%lf %lf %lf %lf %lf", &x, &y1, &y2, &y3, &y4);
P[++cnt] = Point(x, y1);
P[++cnt] = Point(x, y2);
P[++cnt] = Point(x, y3);
P[++cnt] = Point(x, y4);
} rep(i, , cnt + ) rep(j, , cnt + ) if(i == j) dis[i][j] = ; else dis[i][j] = INF; rep(i, , cnt) { /// 判断(0,5)是否能直达P[i],P[i]是否能直达(10,5)
bool flag = ;
int up = ((i + ) / ) * + ;
up -= ;
for(int j = ; j < up; j += ) { /// (0,5)是否能直达P[i]
if(SegmentProperInsection(P[j], P[j + ], Point(, ), P[i]) == false && SegmentProperInsection(P[j + ], P[j + ], Point(, ), P[i]) == false) flag = ;
}
if(!flag) dis[][i] = dis[i][] = Length(Point(, ) - P[i]); flag = ;
for(int j = up + ; j <= cnt; j += ) { /// P[i] 是否能直达(10,5)
if(SegmentProperInsection(P[j], P[j + ], Point(, ), P[i]) == false && SegmentProperInsection(P[j + ], P[j + ], Point(, ), P[i]) == false) flag = ;
}
if(!flag) dis[cnt + ][i] = dis[i][cnt + ] = Length(Point(, ) - P[i]);
} rep(i, , cnt) rep(j, i + , cnt) { /// 枚举两点,判断这两点是否能直达
int st = ((i + ) / ) * + ;
int ed = ((j + ) / ) * ;
bool flag = ;
for(int k = st; k <= ed; k += ) {
if(SegmentProperInsection(P[k], P[k + ], P[i], P[j]) == false && SegmentProperInsection(P[k + ], P[k + ], P[i], P[j]) == false) flag = ;
}
if(!flag) dis[i][j] = dis[j][i] = Length(P[i] - P[j]);
} bool flag = ;
for(int i = ; i <= cnt; i += ) /// 判断(0,5)是否能直达(10,5)
if(SegmentProperInsection(P[i], P[i + ], Point(, ), Point(, )) == false && SegmentProperInsection(P[i + ], P[i + ], Point(, ), Point(, )) == false) flag = ;
if(!flag) dis[][cnt + ] = dis[cnt + ][] = ; rep(k, , cnt + ) rep(i, , cnt + ) rep(j, , cnt + )
if(dis[i][k] + dis[k][j] < dis[i][j]) dis[i][j] = dis[i][k] + dis[k][j];
printf("%.2f\n", dis[][cnt + ]);
}
return ;
}

POJ 2556 (判断线段相交 + 最短路)的更多相关文章

  1. POJ_1556_The Doors_判断线段相交+最短路

    POJ_1556_The Doors_判断线段相交+最短路 Description You are to find the length of the shortest path through a ...

  2. POJ 1556 计算几何 判断线段相交 最短路

    题意: 在一个左下角坐标为(0,0),右上角坐标为(10,10)的矩形内,起点为(0,5),终点为(10,5),中间会有许多扇垂直于x轴的门,求从起点到终点在能走的情况下的最短距离. 分析: 既然是求 ...

  3. C - Segments POJ - 3304 (判断线段相交)

    题目链接:https://vjudge.net/contest/276358#problem/C 题目大意:给你n条线段,问你是否存在一条线段使得所有的线段在这条直线的投影至少具有一个交点? 具体思路 ...

  4. POJ 3449 /// 判断线段相交

    题目大意: 给出多个多边形及其编号 按编号顺序输出每个多边形与其相交的其他多边形编号 注意一个两个多个的不同输出 将每个多边形处理成多条边 然后去判断与其他多边形的边是否相交 计算正方形另外两点的方法 ...

  5. 简单几何(线段相交+最短路) POJ 1556 The Doors

    题目传送门 题意:从(0, 5)走到(10, 5),中间有一些门,走的路是直线,问最短的距离 分析:关键是建图,可以保存所有的点,两点连通的条件是线段和中间的线段都不相交,建立有向图,然后用Dijks ...

  6. POJ 2826 An Easy Problem? 判断线段相交

    POJ 2826 An Easy Problem?! -- 思路来自kuangbin博客 下面三种情况比较特殊,特别是第三种 G++怎么交都是WA,同样的代码C++A了 #include <io ...

  7. 【POJ 2653】Pick-up sticks 判断线段相交

    一定要注意位运算的优先级!!!我被这个卡了好久 判断线段相交模板题. 叉积,点积,规范相交,非规范相交的简单模板 用了“链表”优化之后还是$O(n^2)$的暴力,可是为什么能过$10^5$的数据? # ...

  8. POJ 2653 Pick-up sticks(判断线段相交)

    Pick-up sticks Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 7699   Accepted: 2843 De ...

  9. POJ 1066--Treasure Hunt(判断线段相交)

    Treasure Hunt Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 7857   Accepted: 3247 Des ...

随机推荐

  1. AcWing 787.归并排序

    AcWing 787.归并排序 题目描述 给定你一个长度为n的整数数列. 请你使用归并排序对这个数列按照从小到大进行排序. 并将排好序的数列按顺序输出. 输入格式 输入共两行,第一行包含整数 n. 第 ...

  2. vuex之getter(二)

    说明 使用vue,如果想对data数据派生一些状态,我们就用到计算属性或者侦听器,同样vux想要派生state中的一些状态,可以在store中定义一个getters属性,它相当于state的计算属性. ...

  3. java编写杨辉三角

    import java.util.Scanner; /* *计算杨辉三角: * 规律:两边都是1 * 从第三行开始,上一行的前一个元素+与其并排的元素等于下面的元素 * 例如: * 1 * 11 * ...

  4. 64位Win7下用VS2010编译OpenSSL

    http://blog.csdn.net/henter/article/details/8364532 首先要说明,如果不想自己编译OpenSSL,可以到 http://slproweb.com/pr ...

  5. 二、UDP

    1 端口号 在计算机网络和电磁信号理论中,对共享同一通信信道的多个信号进行区分是个常见的问题.多路复用(multiplexing)就是允许多个会话共享同一介质或机制的一种解决方案. 使用不同的频率来区 ...

  6. HDU_1864_01背包

    http://acm.hdu.edu.cn/showproblem.php?pid=1864 题目好像是输入的数据都是两位小数,先统计能报销的发票,然后把小数*100变成成熟就是01背包问题了. #i ...

  7. [转载]goldendict下优质词典简介及安装

    使用Arch Linux一年以来,如果要问自己最为中意的词典程序是? 当然是Goldendict啦!想详细了解这款瑞士军刀的请猛戳这里. 以前在Win下都是用的lingoes, 感觉还不错,词典库很全 ...

  8. ionic2的返回按钮的编辑问题

    ionic2 返回按钮 首先可以在 app.module.ts 文件中配置. @NgModule 中的 imports 属性的 IonicModule.forRoot 第二个参数,如下: IonicM ...

  9. java设计模式7——桥接模式

    java设计模式7--桥接模式 1.桥接模式介绍 桥接模式是将抽象部分与它的实现部分分离,使他们都可以独立的变化.它是一种对象结构型模式,又称为柄体模式或接口模式. 2.解决问题 2.1.将复杂的组合 ...

  10. layui表格增删改查与上传图片+Api

    API  控制器1 主要用于增删改查已经反填数据查询 using System; using System.Collections.Generic; using System.Data.SqlClie ...