POJ1039 Pipe

最新推荐文章于 2013-09-15 19:14:30 发布

二十一画生

最新推荐文章于 2013-09-15 19:14:30 发布

阅读量863

点赞数 1

分类专栏：计算几何文章标签： POJ

本文链接：https://blog.csdn.net/qhc_20130731/article/details/8941421

版权

计算几何专栏收录该内容

37 篇文章 0 订阅

订阅专栏

//Time 47ms,Memory 248k

#include<iostream>
#include<cmath>
using namespace std;
double d[22];
typedef struct point 
{
	double x,y;
	point(double x=0,double y=0):x(x),y(y){}
}vector;
point p[22];
double cross(vector a,vector b)
{
	return a.x*b.y-a.y*b.x;
}
int f(int a,double s,point i,point j,double &right)
{
	vector v(j.x-i.x,j.y-i.y),w(1,s);
	vector u(i.x-p[a].x,i.y-p[a].y-d[a]);
	double t=cross(v,w),l;
	if(!t) 
		return 0;
	t=cross(w,u)/t;
	if(t>=0 && t<=1) 
	{
		l=i.x+v.x*t;
		if(l>right) right=l;
		return 1;
	}
	return 0;
}
int f1(point a,point b,point c)
{
	vector v1(b.x-a.x,b.y-a.y),v2(c.x-b.x,c.y-b.y);
	double t=cross(v1,v2);
	if(t>0) 
		return 1;
	return 0;
}
void f2(int i,int j,double m,double &right)
{
	int ang=0;
	if(j!=i+1 || d[i])
		ang=f(i,m,p[j-1],p[j],right);
	if(!ang && (j!=i+1 || !d[i]))
		ang=f(i,m,point(p[j-1].x,p[j-1].y-1),point(p[j].x,p[j].y-1),right);
}
int main()
{
	int i,j,n;
	double ang,ang1,ang2,max,min,right;
	while(cin>>n && n)
	{
		for(i=0;i<n;i++)
			cin>>p[i].x>>p[i].y;
		for(i=1;i<n-1;i++)
		{
			if(f1(p[i-1],p[i],p[i+1])) d[i]=0;
			else d[i]=-1;
		}
		if(p[n-1].y<p[n-2].y) d[n-1]=0;
		else d[n-1]=-1;
		right=p[0].x;
		for(i=1;i<n;i++)
		{
			max=(p[i].y+d[i]-p[0].y)/(p[i].x-p[0].x);
			min=(p[i].y+d[i]-p[0].y+1)/(p[i].x-p[0].x);
			for(j=1;j<i;j++)
			{
				ang=(p[i].y+d[i]-p[j].y-d[j])/(p[i].x-p[j].x);
				if(d[j])
				{
					if(min>ang) min=ang;
				}
				else if(max<ang) max=ang;
				if(max>min) break;
			}
			if(max<=min)
			{
				if(i==n-1) 
				{
					right=p[i].x;break;
				}
				ang1=ang2=1;
				for(j=i+1;j<n;j++)
				{
					if(ang1)
					{
						ang1=f(i,max,p[j],point(p[j].x,p[j].y-1),right);
						if(!ang1)
							f2(i,j,max,right);
					}
					if(ang2)
					{
						ang2=f(i,min,p[j],point(p[j].x,p[j].y-1),right);
						if(!ang2)
							f2(i,j,min,right);
					}
					if(right==p[n-1].x) break;
					if(!ang1 && !ang2) break;
				}
				if(j==n && (ang1 || ang2))
				{
					right=p[j-1].x;break;
				}
			}
			if(right==p[n-1].x) break;
		}
		if(right==p[n-1].x) printf("Through all the pipe.\n");
		else printf("%.2lf\n",right);
	}
	return 0;
}