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

时间:2022-08-16 00:56:19

题目没什么

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

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