poj 3304 线段与直线相交判断

别人的题解很详细，我就不献丑了，下面题解的连接： http://hi.baidu.com/cloudayc/item/243f3c4c2a687aaadf2a9fb7

证明：坐标系内若干线段，是否存在一条直线与每条线段都有交点

转成叉积计算线段和直线是否相交,注意重点

转自discuss上的证明：from hanjialong

首先题中的要求等价于：存在一条直线l和所有的线段都相交。
证明：若存在一条直线l和所有线段相交，作一条直线m和l垂直，则m就是题中要求的直线，所有线段投影的一个公共点即为垂足。（l和每条线段的交点沿l投影到m上的垂足处）
反过来，若存在m，所有线段在m上的投影有公共点，则过这点垂直于m作直线l，l一定和所有线段相交。

然后证存在l和所有线段相交等价于存在l过某两条线段的各一个端点且和所有线段相交。
充分性显然。必要性：若有l和所有线段相交，则可保持l和所有线段相交，左右平移l到和某一线段交于端点停止（“移不动了”）。然后绕这个交点旋转。也是转到“转不动了”（和另一线段交于其一个端点）为止。这样就找到了一个新的l。

于是本题可归结为枚举两两线段的各一个端点，连一条直线，再判断剩下的线段是否都和这条直线有交点。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const int mx = 200;
const double eps = 1e-8;
struct po
{
    double x,y;
} point1[mx],point2[mx];
int lineintersect(po p1, po p2, po p3, po p4)
{
    po tp1, tp2, tp3;
    tp1.x = p1.x - p3.x;
    tp1.y = p1.y - p3.y;
    tp2.x = p4.x - p3.x;
    tp2.y = p4.y - p3.y;
    tp3.x = p2.x - p3.x;
    tp3.y = p2.y - p3.y;
//    printf("(tp1.x * tp2.y - tp1.y * tp2.x) * (tp2.x * tp3.y - tp2.y * tp3.x):%lf\n",(tp1.x * tp2.y - tp1.y * tp2.x) * (tp2.x * tp3.y - tp2.y * tp3.x));
    if((tp1.x * tp2.y - tp1.y * tp2.x) * (tp2.x * tp3.y - tp2.y * tp3.x) >= 0)
        return 1;//intersect
    else
        return 0;//not intersect
}

bool check(po p, po q, int n, int a, int b)
{
    if(fabs(p.x - q.x) < eps && fabs(p.y - q.y) < eps) return false;
    for(int i=1; i<=n; i++)
    {
        if(i == a || i == b) continue;
        if(!lineintersect(point1[i],point2[i],p,q))
            return false;
    }
    return true;
}

int main()
{
//    freopen("3304.txt","r",stdin);
    int t,cases = 0;
    scanf("%d",&t);
    while(t--)
    {
        cases++;
        int n;
        scanf("%d",&n);
        for(int i=1; i<=n; i++)
            scanf("%lf %lf %lf %lf",&point1[i].x,&point1[i].y,&point2[i].x,&point2[i].y);
//        for(int i=1; i<=n; i++)
//        {
//            printf("x1 %lf y1 %lf x2 %lf y2 %lf\n",point1[i].x,point1[i].y,point2[i].x,point2[i].y);
//        }
        if(n==1 || n==2)
        {
            printf("Yes!\n");
            continue;
        }

        int vis = 0;
        po p1,p2,q1,q2;
        for(int i=1; i<=n; i++)
        {
            for(int j=i+1; j<=n; j++)
            {
                p1 = point1[i];
                p2 = point2[i];
                q1 = point1[j];
                q2 = point2[j];
                if(check(p1,q1,n,i,j) || check(p1,q2,n,i,j) || check(p2,q1,n,i,j) || check(p2,q2,n,i,j))
                {
                    vis = 1;
                    printf("Yes!\n");
                    break;
                }
            }
            if(vis == 1)    break;
        }
        if(!vis)
            printf("No!\n");
    }
    return 0;
}