終於開始接觸一下計算幾何了，剛入門確實很多地方需要數學法去計算，因為畢竟是電腦，如果用方程聯立這種方法 還是比較困難的，多多瞭解 內積和外積性質 十分有利於解決計算幾何這類問題

題目如下 

給出 n 個點座標 每兩個個點qx qy ,  px py 代表一根木棍，現在開始詢問 輸入 x,y 直至輸入 0 0 問 X Y 是否相連 間接也算。

先拋開幾何部分不談，問兩個節點是否相連，還涉及到間接，就很容易讓人想到並查集，雖然感覺書上說的那種 Floyd 連通圖解決方法也很不錯， 所以關鍵點 就在於 判斷兩根木棍是否相連，當然如果不直接相連也沒關係，畢竟我們這裡在維護並查集，最後只需要簡單的判斷是否同根就可以出結果了

以下。。其實並不是AC程式碼。。我自己找了下資料 都過了。。但是 POJ 好像莫名其妙出事了。。莫名502.。。

自信的感覺能AC。。如果不能。。。就DEBUG去了。。

#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 x)
{
    int r=x;
    while(r!=pre[r])
        r=pre[r];
    return r;
}
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);
}
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((p[i]-q[i]).det(p[j]-q[j]) == 0) //平行情況
                {
                     if(on_seg(p[i],q[i],p[j]) || on_seg(p[i],q[i],q[j]) ||
                        on_seg(p[j],q[j],p[i]) || on_seg(p[j],q[j],q[i]))   
                        unite(i, j);
                }
                else           //看看公共點是不是線上段上
                {
                    point temp = inter(p[i],q[i],p[j],q[j]);
                    if(on_seg(p[i], q[i], temp) && on_seg(p[j], q[j], temp))
                        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;
}