描述

在一个长宽均为10,入口出口分别为(0,5)、(10,5)的房间里,有几堵墙,每堵墙上有两个缺口,求入口到出口的最短路经。

格式

输入格式

第一排为n(n<=20),墙的数目。

接下来n排,每排5个实数x,a1,b1,a2,b2。

x表示墙的横坐标(所有墙都是竖直的),a1-b1和a2-b2之间为空缺。

a1、b1、a2、b2保持递增,x1-xn也是递增的。

输出格式

输出最短距离,保留2位小数。

样例1

样例输入1[复制]

2
4 2 7 8 9
7 3 4.5 6 7

样例输出1[复制]

10.06


【题解】

这是一道最短路问题。

这样做。

给每个墙上的两个“裂缝”的端点标号。则每个墙有4个端点(起点和终点除外)。

则第i个墙的第j个端点它的编号为(i - 1) + j;

然后起点的编号为0,终点的编号为4*n + 1;

最后输出0到4*n+1的最短距离;

然后建图的过程如下

枚举第ii个墙和第jj个墙

然后第ii个墙枚举的是第i个端点,第jj个墙枚举的是第j个端点。

设两个端点的坐标为(x1,y1),(x2,y2);

然后我们就要看看ii+1..jj-1这些墙里面是否有会挡住这两点所构成的直线的墙。

如果有这样的墙,则这两个点不能连线,否则在它们两个点之间建立一条边。

这条边的权值就是两点之间的距离。

可以用两点之间的距离公式求得。

sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));

然后要判断会不会有墙挡住这条直线则需要用到直线方程的两点式;

已知两个点(x1,y1),(x2,y2);

则它们所构成的直线方程为

(x-x1)/(x2-x1)=(y-y1)/(y2-y1)

然后两边消一下,换一下。

得到y=(y2-y1)*(x-x1)/(x2-x1) + y1

然后假设我们枚举到ii+1..jj-1中的墙k。

则x=墙k的横坐标

代入之后得到y0;

看看这个y0是不是在k墙的两个窟窿中。如果不是则会撞墙返回false跳过这条边。

如果ii+1..jj-1这些墙都不会阻碍这条直线。则可以在这两个点之间建一条边。

然后还要单独处理起点到其他墙的、其他墙到终点的以及起点到终点的这3种情况。

最后选择用floyd来求最短路即可(简单!)

【代码】

#include <cstdio>
#include <cmath>
#include <cstring> struct point
{
double x;
double height[5];
}; int n, t[200][200] = { 0 };
double w[200][200];
point a[200]; int main()
{
for (int i = 0; i <= 199; i++)
for (int j = 0; j <= 199; j++) //这是做floyd算法用的数组。一开始赋值为无穷大。
w[i][j] = 2100000000 / 3;
scanf("%d", &n);
for (int i = 1; i <= n; i++) //输入n个墙的信息。
{
double x, a1, b1, a2, b2;
scanf("%lf%lf%lf%lf%lf", &x, &a1, &b1, &a2, &b2);
a[i].x = x, a[i].height[1] = a1, a[i].height[2] = b1, a[i].height[3] = a2, a[i].height[4] = b2;
} for (int ii = 2; ii <= n; ii++) //第ii列
{
for (int i = 1;i <= 4;i++)
for (int jj = 1; jj <= ii - 1; jj++) //第jj列
{
for (int j = 1; j <= 4; j++) //这层要先求出i,j两个点构成的直线方程。
{
double x1, y1, x2, y2;
x1 = a[ii].x;y1 = a[ii].height[i]; //这是所枚举的两个点的坐标。
x2 = a[jj].x;y2 = a[jj].height[j];
bool judge = true;
for (int kk = jj + 1; kk <= ii - 1; kk++)//它们中间的墙。
{
double x = a[kk].x;
double y0 = ((y2 - y1)*(x - x1) / (x2 - x1)) + y1; //这是这两列的相应的点的直线方程,把第kk列的点的信息代入。
if ((y0 >= a[kk].height[1] && y0 <= a[kk].height[2]) || (y0 >= a[kk].height[3] && y0 <= a[kk].height[4]))
judge = true; //如果在两个窟窿里面则是可以添加的边
else//否则只要有一个不行就判错
{
judge = false;
break;//结束循环了。
}
}
if (judge) //如果可以添加这条边
{
int fr = (jj - 1) * 4 + j, to = (ii - 1) * 4 + i;//这是两个点的编号。
t[fr][0]++; //起点的出度增加。
t[fr][t[fr][0]] = to;
w[fr][to] = sqrt((x1 - x2)*(x1 - x2) + (y1 - y2)*(y1 - y2));//权值为两点之间的距离。
}
}
}
} //其实下面这两段直接复制上面的同时枚举两个点的片段即可。 for (int jj = 1; jj <= n; jj++) //这是起点到其他墙的点。(不包括终点)
{
for (int j = 1; j <= 4; j++)
{
double x1, y1, x2, y2;
x1 = 0; y1 = 5; //这下其中的一个点变成是固定的了。
x2 = a[jj].x; y2 = a[jj].height[j];//现在只需要枚举一个点了。
bool judge = true;
for (int kk = 1; kk <= jj- 1; kk++)//然后枚举它们中间的墙,如果它们中间没有墙,则说明一定可以添加,所以初值是true;
{
double x = a[kk].x;
double y0 = ((y2 - y1)*(x - x1) / (x2 - x1)) + y1;//用同样的办法判断是否可以添加这条边。
if ((y0 >= a[kk].height[1] && y0 <= a[kk].height[2]) || (y0 >= a[kk].height[3] && y0 <= a[kk].height[4]))
judge = true;
else
{
judge = false;
break;
}
}
if (judge)
{
int fr = 0, to = (jj - 1) * 4 + j;//如果可以就添加这条边。
t[fr][0]++;
t[fr][t[fr][0]] = to;
w[fr][to] = sqrt((x1 - x2)*(x1 - x2) + (y1 - y2)*(y1 - y2));
}
}
} for (int jj = 1; jj <= n; jj++) //其他墙上的点(不包括起点)到终点(4*n+1)
{
for (int j = 1; j <= 4; j++)
{
double x1, y1, x2, y2;
x1 = a[jj].x; y1 = a[jj].height[j];
x2 = 10; y2 = 5;//这下其中一个固定的点变成了终点
bool judge = true;
for (int kk = jj+1; kk <= n; kk++)
{
double x = a[kk].x;
double y0 = ((y2 - y1)*(x - x1) / (x2 - x1)) + y1;
if ((y0 >= a[kk].height[1] && y0 <= a[kk].height[2]) || (y0 >= a[kk].height[3] && y0 <= a[kk].height[4]))
judge = true;
else
{
judge = false;
break;
}
}
if (judge) //如果可以则添加这条边。
{
int fr = (jj - 1) * 4 + j, to = 4*n+1;
t[fr][0]++;
t[fr][t[fr][0]] = to;
w[fr][to] = sqrt((x1 - x2)*(x1 - x2) + (y1 - y2)*(y1 - y2));
}
}
} //起点到终点,因为可能可以直接从起点画一条直线就到终点了。所以一定要判断这个。不然会WA!
double x1, y1, x2, y2;
x1 = 0; y1 = 5;
x2 = 10; y2 = 5; //两个点都是固定的了。
bool judge = true;
for (int kk = 1; kk <= n; kk++)//看看起点和终点直接的墙会不会阻碍它们俩画一条直线。
{
double x = a[kk].x;
double y0 = ((y2 - y1)*(x - x1) / (x2 - x1)) + y1;
if ((y0 >= a[kk].height[1] && y0 <= a[kk].height[2]) || (y0 >= a[kk].height[3] && y0 <= a[kk].height[4]))
judge = true;
else
{
judge = false;
break;
}
}
if (judge)//如果不会被阻碍,就画一条边。
{
int fr = 0, to = 4 * n + 1;
t[fr][0]++;
t[fr][t[fr][0]] = to;
w[fr][to] = sqrt((x1 - x2)*(x1 - x2) + (y1 - y2)*(y1 - y2));
} n = 4 * n + 1;//让n变成节点的编号。
for (int k = 0; k <= n; k++)//从0-n做floyd;
for (int i = 0; i <= n; i++)
for (int j = 0; j <= n; j++)
if (w[i][j] > w[i][k] + w[k][j])
w[i][j] = w[i][k] + w[k][j];
printf("%.2lf", w[0][n]); //保留两位小数输出即可。
return 0;
}

最新文章

  1. Scalaz(3)- 基础篇:函数概括化-Generalizing Functions
  2. yum服务器设置
  3. ORACLE CHECK CONSTRAINT使用示例(转载) .
  4. ChartType属性
  5. hdu1005 Number Sequence
  6. python字符串中的中文处理
  7. codevs 3342 绿色通道
  8. AngularJS2+调用原有的js脚本(AngularJS脚本跟本地原有脚本之间的关系)
  9. 51行代码实现简单的PHP区块链
  10. 数据结构与算法(C#)入门 --- 序
  11. 阻止事件冒泡传播stopPropagation() 阻止自身默认行为preventdefault()
  12. css实现文本溢出显示省略号
  13. emWin 文字图形同时刷新导致图形显示异常
  14. Mysql数据备份与还原
  15. 研究jenkins集成unittest成图
  16. python之numpy的基本使用
  17. DHTMLX学习总结
  18. Mina的服务器
  19. unity 在Game视图中显示Gizmos
  20. jQueryValidate实现重复性验证

热门文章

  1. [置顶] MVC三层架构在各框架中的特征
  2. 《WPF》Expander控件简单美化
  3. SQLite-SQLiteDatabase 数据库实例练习
  4. java list 容器的ConcurrentModificationException
  5. Spark 概念学习系列之Spark存储管理机制
  6. Hbase技术详细学习笔记
  7. 接口如何使用(以笑话大全api为例)
  8. EXT gridPanel 添加图片
  9. jquery设置attr属性值
  10. 2013腾讯编程马拉松初赛第〇场(HDU 4503) 湫湫系列故事——植树节