【计算几何初步-代码好看了点线段相交】【HDU2150】Pipe

题目没什么 

只是线段相交稍微写的好看了点

Pipe

Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 1252    Accepted Submission(s): 481

Problem Description

经过激烈的争夺，Lele终于把那块地从Yueyue的手里抢了回来。接下来，Lele要开始建造他的灌溉系统。

通过咨询Lele的好友——化学系的TT，Lele决定在田里挖出N条沟渠，每条沟渠输送一种肥料。

每条沟渠可以看作是一条折线，也就是一系列线段首尾连接而成(除了第一条线段开头和最后一条线段的结尾)。由于沟渠很细，你可以忽略掉它的宽度。

由于不同的肥料之间混合会发生化学反应，所以修建的沟渠与沟渠之间不能相交。

现在TT给Lele画了一些设计图，Lele请你判断一下设计图中的沟渠与沟渠之间是否有相交。

 

Input

本题目包含多组测试，请处理到文件结束(EOF)。
每组测试的第一行有一个正整数N(0<N<30)，表示管道的数目。接下来给出这N条管道的信息。
对于每条管道，第一行是一个正整数K(0<K<100)，表示这条管道是由K个端点组成。
接下来的K行给出这K个端点信息。每个端点占一行，用两个整数X和Y(0<X,Y<1000)分别表示这个端点的横坐标和纵坐标的值。

 

Output

对于每组测试，如果该测试管道与管道之间有相交的话，输出"Yes"，否则输出"No"。

 

Sample Input

2
2
0 0
1 1
2
0 1
1 0
2
2
0 0
1 1
2
1 0
2 1
2
3
0 0
1 1
2 1
2
2 0
3 0

 

Sample Output

Yes
No
No

 

Author

Linle

 

#include <cstdio>  
#include <cstdlib>  
#include <cmath>  
#include <cstring>  
#include <ctime>  
#include <algorithm>  
#include <iostream>
#include <sstream>
#include <string>
#define oo 0x13131313   
using namespace std;
struct point{
	long long x,y;
};
int OK=1;
long long crossdet(long long x1,long long y1,long long x2,long long y2)
{
	return x1*y2-x2*y1;
}
long long cross(point a,point b,point c)
{
	return crossdet(b.x-a.x,b.y-a.y,c.x-a.x,c.y-a.y);
}
long long dotdet(long long x1,long long y1,long long x2,long long y2)
{
	return x1*x2+y1*y2;
}
long long dot(point a,point b,point c)
{
	return dotdet(b.x-a.x,b.y-a.y,c.x-a.x,c.y-a.y);
}
point A[40][110];
int N;
int K[40];
void input()
{
	OK=1;
	for(int i=1;i<=N;i++)
	{
		cin>>K[i];
		for(int j=1;j<=K[i];j++)
		{
			scanf("%I64d%I64d",&A[i][j].x,&A[i][j].y);
		}
	}
}
int pan(point &s1,point &e1,point &s2,point &e2)
{
	//跨立实验
	long long a1=cross(s1,e1,s2);                //s1-e1 s2 
	long long a2=cross(s1,e1,e2);                //s1-e1 e2
	long long a3=cross(s2,e2,s1); 				//s2-e2 s1
	long long a4=cross(s2,e2,e1);               //s2-e2 s2
	if(a1*a2<0&&a3*a4<0)   return 1;    //均在两端
	if((a1==0&&dot(s2,s1,e1)<0)||
	   (a2==0&&dot(e2,s1,e1)<0)||
	   (a3==0&&dot(s1,s2,e2)<0)||
	   (a4==0&&dot(e1,s2,e2)<0))
	   return 1;
	return 0; 
}
int solve()
{
	for(int i=1;i<=N;i++)
	 for(int j=i+1;j<=N;j++)
	   for(int k=1;k<=K[i]-1;k++)
	    for(int t=1;t<=K[j]-1;t++)
	    {
	    	if(pan(A[i][k],A[i][k+1],A[j][t],A[j][t+1])==1) {OK=0;return 0;}
		}
}
void init()
{
	freopen("a.in","r",stdin);
	freopen("a.out","w",stdout); 
}
int main()
{
//	init();
	while(cin>>N)
	{
		input();
		solve();
		if(OK) cout<<"No"<<endl;
		else cout<<"Yes"<<endl;
	}
	return 0;
}