问题:黑书P359页。
有一个宽度为1的折线管道,上面顶点为(x0,y0).......(xn,yn),下面各顶点为(x0,y0-1)......(xn,yn-1),假设管壁不透明,不反射,光线从最左边的(x0,y0) (x0,y0-1)之间射入,向四面八方传播,问光线最远能射到哪里或着能穿透整个管道。
1.最优光线必须是擦到一个上顶点一个下顶点。
2.枚举算法,任取一个上顶点和一个下顶点,形成直线L,若L能射穿入口,则是一条可行的光线,在从左向右判断是否相交的管壁交点,并不断更新保存,如果和所有管壁都不与L相交,说明L射穿了整个管道。
#include <iostream>
#include <cstring>
#include <string>
#include <cmath>
#include <stdio.h>
#include <vector>
using namespace std;
#define max(a,b) ((a)>(b)?(a):(b))
const double INF = 999999999.0;
int dblcmp(double d)//控制精度
{
if(fabs(d) < 1e-11)
return 0;
return (d>0)?1:-1;
}
struct point
{
double x,y;
};
double det(double x1,double y1,double x2,double y2)
{
return x1*y2-x2*y1;
}
double cross(point a,point b,point c)
{
return det(b.x-a.x,b.y-a.y,c.x-a.x,c.y-a.y);
}
int segcross(point a,point b,point c,point d,point &p)//判断AB这条直线是否和CD线段相交
{
double s1,s2;
int d1,d2;
d1 = dblcmp(s1=cross(a,b,c));
d2 = dblcmp(s2=cross(a,b,d));
if(d1*d2<0)
{
p.x = (c.x*s2-d.x*s1)/(s2-s1);
return 1;
}
if(d1==0)//经过拐点
{
p.x = c.x;
return 1;
}
if(d2==0)//经过拐点
{
p.x = d.x;
return 1;
}
return 0;
}
int main()
{
point record;
int n,i,j,k;
while(scanf("%d",&n) == 1 && n)
{
vector<point> up(n);
vector<point> down(n);
up.clear();
down.clear();
for(i = 0; 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;
}
double res = -INF;
for(i = 0; i < n; i++)
{
for(j = 0; j < n; j++)
{
for(k = 0; k < n; k++)
{
if( !segcross(up[i],down[j],up[k],down[k],record) )//未穿过( up[k], down[k] )直线
{
break;
}
}
if(k == n)//穿过了所有( up[k], down[k] ),则穿过了整个管道
res = -1;
else if(k > max(i,j))//未穿过,但是是合法的光线,寻找最远交点
{
record.x = -INF;
segcross(up[i],down[j],up[k-1],up[k],record);//判断是否与上壁相交。
if(record.x > res) res = record.x;
segcross(up[i],down[j],down[k-1],down[k],record);//判断是否与下壁相交。
if(record.x > res) res = record.x;
}
}
}
if(res == -1)
printf("Through all the pipe.\n");
else
printf("%.2lf\n",res);
}
return 0;
}