HDU-1086 You can Solve a Geometry Problem too

最新推荐文章于 2019-11-04 21:42:16 发布

Mosu_

最新推荐文章于 2019-11-04 21:42:16 发布

阅读量370

点赞数

分类专栏：基础题目

本文链接：https://blog.csdn.net/Mosu_/article/details/45484779

版权

基础题目专栏收录该内容

22 篇文章 0 订阅

订阅专栏

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <math.h>
#include <algorithm>
using namespace std;

struct jd
{
    double x1,y1,x2,y2;
};
int main()
{
    int n,i,cnt,k,q;
    double a,b,mi1,ma1,mi2,ma2;
    struct jd p[105];
    while (scanf("%d",&n),n!=0)
    {
        k=0;cnt=0;
        scanf("%lf%lf%lf%lf",&p[k].x1,&p[k].y1,&p[k].x2,&p[k].y2);
        k++;
        while (k<n)
        {
            q=0;
            scanf("%lf%lf%lf%lf",&p[k].x1,&p[k].y1,&p[k].x2,&p[k].y2);
            for (i=0;i<k;i++)
            {
                a=((p[k].x1-p[i].x1)*(p[i].y2-p[i].y1)-(p[k].y1-p[i].y1)*(p[i].x2-p[i].x1))*
                    ((p[i].x2-p[i].x1)*(p[k].y2-p[i].y1)-(p[i].y2-p[i].y1)*(p[k].x2-p[i].x1));
                b=((p[i].x1-p[k].x1)*(p[k].y2-p[k].y1)-(p[i].y1-p[k].y1)*(p[k].x2-p[k].x1))*
                    ((p[k].x2-p[k].x1)*(p[i].y2-p[k].y1)-(p[k].y2-p[k].y1)*(p[i].x2-p[k].x1));
                if (a>=0&&b>=0)
                    cnt++;
            }
            k++;
        }
        printf("%d\n",cnt);
    }
    return 0;
}

判断两线段是否相交，可以把两条线段换成向量，如果线段P1P2和直线Q1Q2相交，则P1P2跨立Q1Q2，即：( P1 - Q1 ) × ( Q2 - Q1 ) * ( Q2 - Q1 ) × ( P2 - Q1 ) >= 0，
如果线段Q1Q2和直线P1P2相交，则Q1Q2跨立P1P2，即：( Q1 - P1 ) × ( P2 - P1 ) * ( P2 - P1 ) × ( Q2 - P1 ) >= 0。