tzoj 3348: 线段相交Ⅲ

描述

线段相交有两种情形:一种是“规范相交”,另一种是“非规范相交”。规范相交是指两条线段恰有唯一一个不是端点的公共点。即如果一条线段的端点在另一条线段上则不视为相交。如果两条线段有部分重合,也不视为相交。而非规范相交则把以上两种情况都视为相交。如下图所示:

规范相交认为a,b两种情况都是不相交的,而非规范相交认为a,b两种情况都是相交的。

本题要求判断两条线段是否相交。如果是规范相交则输出YES,并输出交点坐标,如果是非规范相交则只需输出YES,如果不相交则输出NO。

输入

输入有多组数据,T表示输入数据的组数。每组测试数据有两行第一行输入一条线段的两个端点的坐标,第二行输入另一个线段的两个端点的坐标。

输出

对于每组测试数据,输出一行。如果是规范相交则输出YES,并输出交点坐标(小数点后面保留3位),如果是非规范相交则只需输出YES,如果不相交则输出NO。

样例输入

4 0 0 1 1 0 1 1 0 0 0 2 2 2 2 3 3 0 0 2 2 1.5 1.5 3 3 0 0 1 1 2 2 3 3 

样例输出

YES (0.500,0.500)
YES
YES
NO

题解

先进行快速排斥判断两个线段形成的最小矩形是否相交

 再进行跨立实验,利用叉乘判断c和d是否在ab的两侧即(\overrightarrow {ba}\times\overrightarrow {bc})*(\overrightarrow {ba}\times\overrightarrow {bd})<0

同理另一条线段也要进行这样的判断,当结果等于零时两条线段共线,其中一个等于零时,有一个端点在另一条线段上

可以利用面积求交点,Sabc/Sabd就是交点在线段cd上的比例

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef pair<int,int> PII;
typedef unsigned long long ULL;
    
mt19937 rnd(random_device{}());
uniform_int_distribution<int> dist(0,INT_MAX);
//#define int long long
const int N = 400010,M = 1000010,mod = 998244353,INF = 1e8;

struct Node 
{
    double x,y;
};

double cross(Node a,Node b,Node c)
{
    double x1 = a.x-b.x;
    double y1 = a.y-b.y;
    double x2 = c.x-b.x;
    double y2 = c.y-b.y;
    return x1*y2-x2*y1;
}
int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        Node a,b,c,d;
        cin>>a.x>>a.y>>b.x>>b.y;
        cin>>c.x>>c.y>>d.x>>d.y;
        if(min(a.x,b.x)>max(c.x,d.x)||max(a.x,b.x)<min(c.x,d.x)||min(a.y,b.y)>max(c.y,d.y)||max(a.y,b.y)<min(c.y,d.y))
        {
            puts("NO");
            continue;
        }
        double t1 = cross(a,b,c)*cross(a,b,d);
        double t2 = cross(c,d,a)*cross(c,d,b);
        if(t1<=0&&t2<=0)
        {
            printf("YES");
            if(t1!=0&&t2!=0) 
            {
                double s1 = cross(a,b,c),s2 = cross(a,b,d);
                s1 = fabs(s1),s2 = fabs(s2);
                double temp = s1/(s1+s2);
                double x = c.x+temp*(d.x-c.x);
                double y = c.y+temp*(d.y-c.y);
                printf(" (%.3lf,%.3lf)",x,y);
            }
            cout<<endl;
        }
        else 
            puts("NO");
    }
}

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值