POJ 1039 Pipe 简单的叉乘

题意：给你一段弯弯曲曲的管子，问是否存在一条光线可以穿过整条管子，如果存在，则输出“Through all the pipe”,否则输出X轴正方向光线最远能到的X轴坐标。

思路：如果存在这样一束光线，则此束光线必过管子的两个节点。如果不存在，最远的那一丝光线也必过管子的两个的节点，所以剩下的任务就是无厘头的枚举了。

PS:我感觉我已经爱上熬夜了......

#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <cmath>
#include <algorithm>
#include <string>

#define LL long long
#define Max(a,b) (a) > (b) ? (a) : (b)
#define Min(a,b) (a) < (b) ? (a) : (b)
#define EPS (1e-8)
#define Left true
#define Right false

using namespace std;

struct P
{
double x,y;
}p[40][2];

double X_Mul(P a1,P a2,P b1, P b2)
{
P v1 = {a2.x-a1.x,a2.y-a1.y},v2 = {b2.x-b1.x,b2.y-b1.y};
return (v1.x*v2.y - v2.x*v1.y);
}

bool Is_Through(P p1,P p2,int n)
{
int i;
for(i = 0;i < n; ++i)
{
if(X_Mul(p1,p2,p1,p[i][0]) < (-EPS))
return false;
}
for(i = 0;i < n; ++i)
{
if(X_Mul(p1,p2,p1,p[i][1]) > EPS)
return false;
}
return true;
}

double Cal_Position(P a1,P a2,P b1,P b2)
{
double t = fabs(X_Mul(a1,a2,a1,b1))/fabs(X_Mul(a1,a2,b1,b2));

return b1.x + (b2.x - b1.x)*t;
}

double Cal_Len(P p1,P p2,int n,int s,int e)
{
int i;

for(i = 0;i < s; ++i)
{
if(X_Mul(p1,p2,p1,p[i][0]) < (-EPS) || X_Mul(p1,p2,p1,p[i][1]) > EPS )
return p[0][0].x;

//此时表明此光线不存在
}
for(i = s+1;i < e; ++i)
{
if(X_Mul(p1,p2,p1,p[i][0]) < (-EPS) || X_Mul(p1,p2,p1,p[i][1]) > EPS )
return p[0][0].x;

//此时表明此光线不存在
}

for(i = e+1;i < n; ++i)
{
if(X_Mul(p1,p2,p1,p[i][0]) < EPS)
{
return Cal_Position(p1,p2,p[i-1][0],p[i][0]);
}
else if(X_Mul(p1,p2,p1,p[i][1]) > EPS )
{
return Cal_Position(p1,p2,p[i-1][1],p[i][1]);
}
}
return p[0][0].x;
}

void output(int n)
{
int i,j;
for(i = 0;i < n; ++i)
{
for(j = i+1;j < n; ++j)
{
if(Is_Through(p[i][0],p[j][0],n) || Is_Through(p[i][1],p[j][0],n) || Is_Through(p[i][0],p[j][1],n)|| Is_Through(p[i][1],p[j][1],n))
{
cout<<"Through all the pipe."<<endl;
return ;
}
}
}
double Temp,Max = p[0][0].x;
for(i = 0;i < n; ++i)
{
for(j = i+1;j < n; ++j)
{
Temp = Cal_Len(p[i][0],p[j][0],n,i,j);
Max = Max(Temp,Max);
Temp = Cal_Len(p[i][1],p[j][0],n,i,j);
Max = Max(Temp,Max);
Temp = Cal_Len(p[i][0],p[j][1],n,i,j);
Max = Max(Temp,Max);
Temp = Cal_Len(p[i][1],p[j][1],n,i,j);
Max = Max(Temp,Max);
}
printf("%.2f\n",Max);
return ;
}

int main()
{
int n,i;
while(scanf("%d",&n) && n)
{
for(i = 0;i < n; ++i)
{
scanf("%lf %lf",&p[i][0].x,&p[i][0].y);
p[i][1] = p[i][0];
p[i][1].y--;
}

output(n);

}
return 0;
}