PKU 1039 Pipe

 

题意理解：

         某公司生产一种弯管（如图），它的材料部不反光也不漫反射更不透明。一束光只能从在弯管中笔直通过，相当于：在一个弯管中从入端用激光照射，当遇到管壁时（不包括转折点），就终止。

         求光线能够到达x最大处的x值；

         其中x1<x2<…<xn;

         每个转折点处上点的纵坐标比下点的纵坐标大1.输入的是上边的点序列。

 

解题思路：

         最有点满足：

                   一定要通过两个转折点，并且1上1下。

         证明（源于黑书）：

                   如果不通过任何转折点，那么可以通过平移使光线使其至少通过一个转折点，并且x增加；

如果只通过一个点，一定可以通过旋转该光线使其通过两个点（一上一下）使x增长。

如果通过的同边（同上边或同下边）的两个，一定可以通过平移时其不通过任何点，使x增加。

        

         有以上性质，那么可以枚举上下各一个点组成一条光线，检验该光线能否通过转折处（从左到右），直到不能通过某个转折点（设为第m个点），说明光线与m-1，到m之间必定与管壁有交点。可以求出交点的x。寻找最大的x即可：#include<iostream> #include<algorithm> #include<stdio.h> #include<cmath> using namespace std; const int N=30; const double esp=1e-6; const int oo=1<<30; typedef struct {double x,y;}point; point up[N],down[N]; double ans; int dblcmp(double x) { if(fabs(x)<esp) return 0; return (x>0)?1:-1; } double crossleft(point a,point b,point c) { double x1,x2,y1,y2; x1=b.x-a.x,y1=b.y-a.y; x2=c.x-a.x,y2=c.y-a.y; return x1*y2-x2*y1; } double intersection(point a,point b,point c,point d) { double k,s1,s2; point ff,tt; ff.x=b.x-(a.x-d.x); ff.y=b.y-(a.y-d.y); s1=crossleft(a,b,c); s2=crossleft(d,ff,c); k=s1/s2; tt.x=c.x+k*(d.x-c.x); return tt.x; } void have(int n) { int i,j,k,s; ans=-oo; for(i=1;i<=n;i++) for(j=1;j<=n;j++) { s=max(i,j); for(k=1;k<=n;k++) if( dblcmp( crossleft(up[i],down[j],up[k]) )* dblcmp( crossleft(up[i],down[j],down[k]) )>0 ) break; if(k>s) { if(k>n) { printf("Through all the pipe./n"); return; } double t; if( dblcmp( crossleft(up[i],down[j],up[k-1]) )* dblcmp( crossleft(up[i],down[j],up[k]) )<0 ) t=intersection(up[i],down[j],up[k-1],up[k]); else t=intersection(up[i],down[j],down[k-1],down[k]); if(ans<t) ans=t; } } printf("%.2lf/n",ans); } int main() { int n; while(scanf("%d",&n)!=EOF&&n) { for(int i=1;i<=n;i++) { scanf("%lf %lf",&up[i].x,&up[i].y); down[i].x=up[i].x; down[i].y=up[i].y-1; } have(n); } return 0; }