POJ 2653 计算几何 判断线段相交

题意：

按顺序给出一些木棍，输出在最上面的木棍标号

题解：

一开始看到n=100000还在想怎么优化，没想出来，最后看讨论，原来暴力可以水过~

嘿嘿~用了下并查集压缩路径，在刷数组的时候会快一些~

#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <algorithm>

#define N 110000
#define EPS 1e-8

using namespace std;

struct P
{
    double x,y;
}li[N][2];

bool vis[N];
int n,fa[N];

inline void read()
{
    for(int i=1;i<=n;i++)
    {
        scanf("%lf%lf%lf%lf",&li[i][0].x,&li[i][0].y,&li[i][1].x,&li[i][1].y);
        fa[i]=i;
    }
}

inline double cross(const P &o,const P &p,const P &q)
{
    double fx=p.x-o.x,fy=p.y-o.y,sx=q.x-o.x,sy=q.y-o.y;
    return fx*sy-fy*sx;
}
inline bool check(int a,int b)
{
    double fg1=cross(li[a][0],li[b][0],li[a][1])*cross(li[a][0],li[b][1],li[a][1]);
    double fg2=cross(li[b][0],li[a][0],li[b][1])*cross(li[b][0],li[a][1],li[b][1]);
    if(fg1<EPS&&fg2<EPS) return true;
    return false;
}

inline int findfa(int u)
{
    if(fa[u]!=u) fa[u]=findfa(fa[u]);
    return fa[u];
}

inline void go()
{
    memset(vis,0,sizeof vis);
    for(int i=1;i<=n;i++)
        for(int j=findfa(1);j<i;j=findfa(j+1))
            if(check(i,j)) fa[j]=j+1,vis[j]=true;
    printf("Top sticks:");
    for(int i=1,fg=0;i<=n;i++)
        if(!vis[i])
        {
            if(fg) printf(", %d",i);
            else printf(" %d",i),fg=true;
        }
    puts(".");
}

int main()
{
    while(scanf("%d",&n),n) read(),go();
    return 0;
}