hdu2150 Pipe 计算几何 判断相交

链接：http://acm.hdu.edu.cn/showproblem.php?pid=2150

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
  
  

思路：挨个判断两条线段是否相交，有个很坑的是int一直过不了，换成double就行了，而且叉积是<0即可，测试例子里应该没有<=0的情况....

代码：

 

#include<iostream>
#include<stdlib.h>
#include<stdio.h>
#include<cmath>
#include<algorithm>
#include<string>
#include<string.h>
#include<set>
#include<queue>
#include<stack>
#include<vector>
#include<functional> 
#include<map>
using namespace std;
const int maxn = 100 + 10;
const int INF = (int)1e9;
const double eps = 1e-9;

struct point {
	double x, y;
};

point P[35][105];//P[i][j]表示第i根管的第j个点
int N, K;
int k[35];

double cross(point a, point b, point c) {
	return (a.x - c.x)*(b.y - c.y) - (b.x - c.x)*(a.y - c.y);
}

bool intersect(point aa, point bb, point cc, point dd) {
	if (max(aa.x, bb.x) >= min(cc.x, dd.x)
		&& max(cc.x, dd.x) >= min(aa.x, bb.x)
		&& max(aa.y, bb.y) >= min(cc.y, dd.y)
		&& max(cc.y, dd.y) >= min(aa.y, bb.y)
		&& cross(cc, bb, aa)*cross(dd, bb, aa) < 0 
		&& cross(aa, dd, cc)*cross(bb, dd, cc) < 0)
		return true;//有交点
	return false;
}

bool fun() {//很好很暴力2333
	for (int i = 0; i < N - 1; i++) {//确定第一个向量
		for (int j = 1; j < k[i]; j++) {
			for (int t = i + 1; t < N; t++) {//确定第二个向量
				for (int m = 1; m < k[t]; m++) {
					if (intersect(P[i][j - 1], P[i][j], P[t][m - 1], P[t][m])) {
						return false;//有交点
					}
				}
			}
		}
	}
	return true;
}

int main()
{
	while (scanf("%d", &N) != EOF) {			
		memset(P, 0, sizeof(P));
		memset(k, 0, sizeof(k));
		for (int i = 0; i < N;i++) {
			scanf("%d", &K);
			k[i] = K;
			for (int j = 0; j < k[i]; j++)
				scanf("%lf %lf", &P[i][j].x, &P[i][j].y);
		}
		//全部录入完毕
		if (N == 1) {
			printf("No\n");
		}
		else {
			bool flag = fun();
			if (flag)printf("No\n");
			else printf("Yes\n");
		}
	}
	//system("pause");
	return 0;
}