判断线段与直线交点——noyj 管道问题

管道问题

时间限制： 3000 ms  |  内存限制： 1000 KB

难度： 3

描述

有一宽度为1的折线管道，如图所示，上面的各个定点为 (x0,y0),(x1,y1),(x2,y2) ……(xn,yn),下面各个定点为(x0,y0-1),(x1,y1-1),(x2,y2-1)……(xn,yn-1)，假设管道都是不透明的，不反射的，光线从左边入口的(x0,y0),(x0,y0-1)之间射入，向四面八方直线传播，问光线最远能射到哪里(x坐标)或能穿透整个管道。

输入

第一行有一个整数2<=n<20（n为零结束），表示有n个管道，接下来的n行表示n个管道的上面坐标(x0,y0),(x1,y1),(x2,y2)……(xn,yn)。测试数据组数大于0小于10000，所有坐标的只不大于10000。

输出

如果光线不能穿过管道，输出最远的x坐标，结果保留两位小数。如果有光线可以透过管道，输出“Through all the pipe.”

样例输入

40 12 24 16 460 12 -0.65 -4.457 -5.5712 -10.817 -16.550

样例输出

4.67Through all the pipe.

来源

POJ1039

上传者

iphxer

思路：

这个题N比较小，所以想到暴力枚举所有可能性。方法是，任取两个折点连接成的直线(（a,b）与（b,a）重复)与所有连接管道上顶点和与之对应的下顶点的线段，判断是否都相交，如果都相交，则可以通过管道；反之则与管道壁相交，需要用到算法求出交点X.

下面安利一下求交点的算法3 

http://blog.csdn.net/qq_36738482/article/details/73252913

#include<bits/stdc++.h>
#define E 1e-10
using namespace std;
struct zz
{
	double x;
	double y;
}p[25][2];
double judge(zz a,zz b,zz c)
{
	return ((b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y));
}
double check(zz a1,zz a2,zz b1,zz b2)//判断是否相交
{
	return judge(a1,a2,b1)*(judge(a1,a2,b2));// A(x1, y1), B(x2, y2)的直线方程为：
                                            // f(x, y) =  (y - y1) * (x1 - x2) - (x - x1) * (y1 - y2) = 0
}
double zhexian(zz a,zz b,zz c)//计算三角形面积公式已知三角形各顶点A（Xi，Yi）（ i=1，2，3)
                              //S=((x2-x1)*(y3-y1)-(x3-x1)*(y2-y1))或 海伦公式
{
    /*double ab,ac,bc,q;
    ab=sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
    ac=sqrt((a.x-c.x)*(a.x-c.x)+(a.y-c.y)*(a.y-c.y));
    bc=sqrt((b.x-c.x)*(b.x-c.x)+(b.y-c.y)*(b.y-c.y));
    q=(ab+bc+ac)/2;
    return sqrt(q*(q-ab)*(q-ac)*(q-bc));*/
	return fabs((b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y))/2;
}
double point(zz a1,zz a2,zz b1,zz b2)//计算交点坐标
{
	double s1,s2;
	s1=zhexian(a1,a2,b1);
	s2=zhexian(a1,a2,b2);
	return (s1*b2.x+s2*b1.x)/(s1+s2);
}
int main()
{
	int n;
	while(scanf("%d",&n)!=EOF)
	{
		int i,j;
		if(n==0) break;
		for(i=0;i<n;i++)
		{
			scanf("%lf%lf",&p[i][0].x,&p[i][0].y);
			p[i][1].x=p[i][0].x;
			p[i][1].y=p[i][0].y-1;
		}
		double mx=-0x3f3f3f3f;//0x7fffffff
		bool flag=false;
		zz a,b;
		for(i=0;i<n&&flag==false;i++)//n组折点
		{
			int x1,x2;
			int tmp=0;
			for(tmp=0;tmp<2;tmp++)//上下
			{
				a=p[i][tmp];
				for(x1=i+1;x1<n&&flag==false;x1++)//该点与任意的折点组合成直线
				{
					for(x2=0;x2<2&&flag==false;x2++)
					{
						b=p[x1][x2];
						if(check(a,b,p[0][0],p[0][1])<E)
						{
							for(j=1;j<n;j++)
							{
								if(check(a,b,p[j][0],p[j][1])>E)//如果a,b两点所连直线在p[j][1]、p[j][0],这个折点处没有交点，那么光线在这附近会有交点。
								{
									double x;
									if(judge(a,b,p[j][0])>0)
										x=point(a,b,p[j-1][1],p[j][1]);
									else
										x=point(a,b,p[j-1][0],p[j][0]);
									if(x>mx)
										mx=x;
									break;
								}
							}
							if(j==n)
								flag=true;
						}
					}
				}
			}
		}
		if(flag==true)
			printf("Through all the pipe.\n");
		else
			printf("%.2lf\n",mx);
	}
	return 0;
}