hdu 2150- Pipe-解析几何

Pipe

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

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

 

Source

ACM程序设计期末考试——2008-01-02（3 教417）

 

Recommend

lcy   |   We have carefully selected several similar problems for you:   1392  1115  2108  1086  2036 

题意：

题目的意思是给定线段的两个端点的坐标，判断同一组的几条线段是否相交，是就yes，不是就no。

知识点：

我们分两步确定两条线段是否相交：
(1)快速排斥试验
设以线段 P1P2 为对角线的矩形为R， 设以线段 Q1Q2 为对角线的矩形为T，
如果R和T不相交，显然两线段不会相交。
(2)跨立试验
如果两线段相交，则两线段必然相互跨立对方。
若P1P2跨立Q1Q2 ，则矢量 ( P1 - Q1 ) 和( P2 - Q1 )位于矢量( Q2 - Q1 ) 的两侧，
即( P1 - Q1 ) × ( Q2 - Q1 ) * ( P2 - Q1 ) × ( Q2 - Q1 ) < 0。
上式可改写成( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) > 0。
当 ( P1 - Q1 ) × ( Q2 - Q1 ) = 0 时，说明 ( P1 - Q1 ) 和 ( Q2 - Q1 )共线，
但是因为已经通过快速排斥试验，所以 P1 一定在线段 Q1Q2上；
同理，( Q2 - Q1 ) ×(P2 - Q1 ) = 0 说明 P2 一定在线段 Q1Q2上。
所以判断P1P2跨立Q1Q2的依据是：( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) >= 0。
同理判断Q1Q2跨立P1P2的依据是：( Q1 - P1 ) × ( P2 - P1 ) * ( P2 - P1 ) × ( Q2 - P1 ) >= 0。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
#define MAX 100+5
typedef long long LL;
const double pi=3.141592653589793;
const int INF=1e9;
const double inf=1e20;
const double eps=1e-6;
const int mod=1000000007;
struct point
{
    int x,y;
} p[35][105];
double circulation(point a,point b,point c)///计算向量BA、CA的叉积
{
    return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}
bool cross(point a,point b,point c,point d)///判断线段AB和CD是否相交
{
    if(max(a.x,b.x)>=min(c.x,d.x)&&///快速排斥试验
            max(c.x,d.x)>=min(a.x,b.x)&&
            max(a.y,b.y)>=min(c.y,d.y)&&
            max(c.y,d.y)>=min(a.y,b.y)&&
            circulation(a,b,c)*circulation(a,d,b)>=0&&///跨立试验
            circulation(c,a,d)*circulation(c,d,b)>=0)
        return true;
    return false;
}
int main()
{
    int n,k[35];
    while(cin>>n)///n表示管道的数目
    {
        memset(k,0,sizeof(k));
        memset(p,0,sizeof(p));
        for(int i=0; i<n; i++)
        {
            cin>>k[i];///i表示是第i条管道，k[i]表示这条管道是由k[i]个端点组成
            for(int j=0; j<k[i]; j++)
                scanf("%d%d",&p[i][j].x,&p[i][j].y);
        }
        if(n==1)
        {
            printf("No\n");
            continue;
        }
        int flag=0;
        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 h=1; h<k[t]; h++)
                        if(cross(p[i][j-1],p[i][j],p[t][h-1],p[t][h]))
                        {
                            flag=1;
                            break;
                        }
        if(flag) printf("Yes\n");
        else printf("No\n");
    }
    return 0;
}