POJ 3304 直线与线段相交的判断

参考:
http://blog.sina.com.cn/s/blog_6635898a0100n2lv.html
的证明过程
可将问题转化为是否存在一条直线经过所有的线段的判断问题

#include<iostream>
#include<cmath>
#include<cstdio> 
using namespace std;

const int MAX = 105;
const double eps = 1e-8;
struct Point{
    double x,y;
}s[MAX],e[MAX];
int n;

double mult(Point sp,Point ep,Point op){
    return (sp.x-op.x)*(ep.y-op.y)-(ep.x-op.x)*(sp.y-op.y);
}

bool solve(Point p1,Point p2){
    if(abs(p1.x-p2.x)<eps && abs(p1.y-p2.y)<eps) return false;
    for(int i=0;i<n;i++)
        if(mult(p1,p2,s[i])*mult(p1,p2,e[i])>eps)   return false;
    return true;
}

int main()
{
    //freopen("input.txt","r",stdin);
    int t;
    cin>>t;
    while(t--){
        cin>>n;
        for(int i=0;i<n;i++)
            cin>>s[i].x>>s[i].y>>e[i].x>>e[i].y;

        bool flag=false;
        if(n<3) flag=true;

        for(int i=0; i<n && !flag ;i++){
            for(int j=i+1; j<n && !flag ;j++){
                if(solve(s[i],s[j])) flag=true;
                else if(solve(s[i],e[j])) flag=true;
                else if(solve(e[i],s[j])) flag=true;
                else if(solve(e[i],e[j])) flag=true;
            }
        }
        if(flag)    cout<<"Yes!"<<endl;
        else        cout<<"No!"<<endl;
    }
}