hdu/hdoj 1086 You can Solve a Geometry Problem too

最新推荐文章于 2017-09-11 22:25:20 发布

chen_kkkk

最新推荐文章于 2017-09-11 22:25:20 发布

阅读量725

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/chen_kkkk/article/details/10544731

计算几何专栏收录该内容

0 篇文章

订阅专栏

本文提供了一种用于判断两条线段是否相交的有效算法，并附带完整的C++实现代码。该算法通过计算向量乘积来判断线段的位置关系，进而确定是否相交。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

计算两条线段是否相交

#include <iostream>  
using namespace std;  
 
struct line 
{  
    double x1,y1,x2,y2;  
}ln[105];  
 
int n;
 
int mul(double a, double b, double c, double d)  
{  
    return a * d - b * c;  
}
 
bool judge(int i, int j)
{
    double a = mul(ln[i].x1-ln[j].x1, ln[i].y1-ln[j].y1, ln[j].x2-ln[j].x1, ln[j].y2-ln[j].y1);
    double b = mul(ln[i].x2-ln[j].x1, ln[i].y2-ln[j].y1, ln[j].x2-ln[j].x1, ln[j].y2-ln[j].y1);
    double c = a * b;
 
    double d = mul(ln[j].x1-ln[i].x1, ln[j].y1-ln[i].y1, ln[i].x2-ln[i].x1, ln[i].y2-ln[i].y1);
    double e = mul(ln[j].x2-ln[i].x1, ln[j].y2-ln[i].y1, ln[i].x2-ln[i].x1, ln[i].y2-ln[i].y1);
    double f = d * e;
 
    if(c <= 0 && f <= 0)
        return 1;
    else
        return 0;
}
 
 
int main()  
{  
 
    while(scanf("%d", &n) && n)  
    {  
        for(int i=1; i<=n; ++i)  
            scanf("%lf %lf %lf %lf", &ln[i].x1, &ln[i].y1, &ln[i].x2, &ln[i].y2);  
        int count=0;  
        for(int i=1; i<n; ++i)  
        {  
            for(int j=i+1; j<=n; ++j)  
                if(judge(i,j)) 
                    ++count;  
        }  
        cout << count << endl;  
    }  
    return 0;  
}