hdu2150 Pipe 计算几何 判断相交

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

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

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

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

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

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

 

[align=left]Input[/align]
本题目包含多组测试，请处理到文件结束(EOF)。

每组测试的第一行有一个正整数N(0<N<30)，表示管道的数目。接下来给出这N条管道的信息。

对于每条管道，第一行是一个正整数K(0<K<100)，表示这条管道是由K个端点组成。

接下来的K行给出这K个端点信息。每个端点占一行，用两个整数X和Y(0<X,Y<1000)分别表示这个端点的横坐标和纵坐标的值。

 

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

[align=left]Sample Input[/align]

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

 

[align=left]Sample Output[/align]

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;
}