zoj 2551 判线段相交

最新推荐文章于 2021-03-16 19:46:49 发布

leolin_

最新推荐文章于 2021-03-16 19:46:49 发布

阅读量450

点赞数

分类专栏：计算几何文章标签： struct

本文链接：https://blog.csdn.net/leolin_/article/details/6624619

版权

计算几何专栏收录该内容

19 篇文章 0 订阅

订阅专栏

本文详细介绍了如何使用C++语言编写一个算法来判断二维平面上的线段是否相交。通过四次叉乘判断两线段是否相交，确保了算法的准确性和效率。实例代码包含输入读取、线段定义、相交判断函数及输出结果的部分，适合初学者理解和实践。

摘要由CSDN通过智能技术生成

#include<iostream>
#include<vector>
#include<map>
#include<stack>
#include<algorithm>
#include<queue>
#include<list>
#include<set>
#include<string.h>
#include<stdlib.h>
#include<math.h>
#include<stdio.h>
#include<ctype.h>
#include<iomanip>

using namespace std;

#define LL long long
#define pi acos(-1)
#define N 50010
#define INF 9999999999
#define eps 1e-8

struct line
{
    double x1,y1,x2,y2;
    bool ok;
}p[100010];
int ans[1100];

bool fuck(line a,line b)///判两线段是否相交，4次叉乘///
{
    double s1=(a.x2-a.x1)*(b.y1-a.y1)-(a.y2-a.y1)*(b.x1-a.x1);double s2=(a.x2-a.x1)*(b.y2-a.y1)-(a.y2-a.y1)*(b.x2-a.x1);
    double s3=(b.x2-b.x1)*(a.y1-b.y1)-(b.y2-b.y1)*(a.x1-b.x1);double s4=(b.x2-b.x1)*(a.y2-b.y1)-(b.y2-b.y1)*(a.x2-b.x1);
    if(s1*s2<=0 && s3*s4<=0)
    return false;
    else
    return true;
}
int main()
{
//  freopen("a.txt","r",stdin);
    int n;
    while(scanf("%d",&n) && n)
    {
        int i,j,k;
        for(i=1;i<=n;i++)
            scanf("%lf%lf%lf%lf",&p[i].x1,&p[i].y1,&p[i].x2,&p[i].y2);
        p[n].ok=true;
        int cnt=1;
        for(i=n-1;i>=1;i--)
        {
                p[i].ok=true;
                for(j=i+1;j<=n;j++)
                if(!fuck(p[i],p[j]))
                {

                    p[i].ok=false;
                    break;
                }
                if(p[i].ok)
                ans[cnt++]=i;
        }
        ans[0]=n;
        printf("Top sticks:");
        for(i=cnt-1;i>=0;i--)
        {
                printf(" %d",ans[i]);
                if(i!=0)
                printf(",");
                if(i==0)
                printf(".\n");
        }
    }
    return 0;
}