POJ 1127 跨立实验和快速排斥实验

还是那道题 但是这次有一个全新的解法了，涉及到的跨立实验和快速排斥实验，在学习部分我有提及过。

不多说直接上一下核心代码

 

直接进行双实验，如果快速排斥不通过 直接PASS 之后再进行跨立实验，如果成了 那么两直线相交

bool judge(point p1,point p2,point q1,point q2)
{
    if( max(p1.x,p2.x)<min(q1.x,q2.x)|| //快速排斥实验
        max(q1.x,q2.x)<min(p1.x,p2.x)||
        max(p1.y,p2.y)<min(q1.y,q2.y)||
        max(q1.y,q2.y)<min(p1.y,p2.y))
        return 0;
    if((q1-p1).det(q2-p1) * (q1-p2).det(q2-p2) <=0 && (p1-q1).det(p2-q1) * (p1-q2).det(p2-q2)<=0) //跨立实验
        return 1;
    else
        return 0;
}

之后。。虽然和本题无太大关系。。但是我又学会一种新的并查集路径压缩。。网上的大牛真实厉害

int findd(int v)
{
   if(pre[v]==v)
        return v;
    else
    {
        pre[v]=findd(pre[v]);//路径压缩
        return pre[v];
    }
}

 

 

以下是完整代码。。不得不说，用了这两个实验进行相交判定，整串代码都显得清爽了不少

（原来的函数并没有删除）

#include<iostream>
#include<cmath>
#include<cstdio>
using namespace std;
const double eps = 1e-10;
int pre[30];
double add(double a,double b)
{
    if(abs(a+b)<eps*(abs(a)+abs(b)))
        return 0;
    else
        return a+b;
}
struct point
{
    double x,y;
    point(){}
    point(double x,double y): x(x),y(y){}
    point operator- (point p)
    {
        return point(add(x, -p.x),add(y, -p.y));
    }
    point operator+ (point p)
    {
        return point(add(x, p.x),add(y, p.y));
    }
    point operator* (double d)
    {
        return point(x*d, y*d);
    }
    double dot(point p)
    {
        return add(x*p.x,y*p.y);
    }
    double det(point p)
    {
        return add(x*p.y,-y*p.x);
    }
};
point p[30],q[30];
bool on_seg(point p1, point p2, point q)
{
    return (p1 - q).det(p2 - q) == 0 && (p1 - q).dot(p2 - q) <= 0;
}
point inter(point p1, point p2, point q1, point q2)
{
    return p1 + (p2 - p1) * ((q2 - q1).det(q1 - p1) / (q2 - q1).det(p2 - p1));
}
int findd(int v)
{
   if(pre[v]==v)
        return v;
    else
    {
        pre[v]=findd(pre[v]);//路径压缩
        return pre[v];
    }

}
void unite(int x,int y)
{
    int fx=findd(x);
    int fy=findd(y);
    if(fx!=fy)
        pre[fy]=fx;
}
bool same(int x,int y)
{
    return findd(x)==findd(y);
}
bool judge(point p1,point p2,point q1,point q2)
{
    if( max(p1.x,p2.x)<min(q1.x,q2.x)|| //快速排斥实验
        max(q1.x,q2.x)<min(p1.x,p2.x)||
        max(p1.y,p2.y)<min(q1.y,q2.y)||
        max(q1.y,q2.y)<min(p1.y,p2.y))
        return 0;
    if((q1-p1).det(q2-p1) * (q1-p2).det(q2-p2) <=0 && (p1-q1).det(p2-q1) * (p1-q2).det(p2-q2)<=0) //跨立实验
        return 1;
    else
        return 0;
}
int main()
{
    int n;
    while(~scanf("%d",&n)&&n)
    {
        for(int i=0;i<=n;i++)
            pre[i]=i;
        for(int i=1;i<=n;i++)
        {
            cin>>p[i].x>>p[i].y;
            cin>>q[i].x>>q[i].y;
        }
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<i;j++)
            {
                if(judge(p[i],q[i],p[j],q[j]))
                    unite(i,j);
            }
        }
        int x,y;
        while(~scanf("%d%d",&x,&y)&&x&&y)
        {
            if(same(x, y))
                printf("CONNECTED\n");
            else
                printf("NOT CONNECTED\n");
        }
    }
    return 0;
}