poj 3304 Segments（计算几何基础）

 

Segments

Time Limit: 1000MS		Memory Limit: 65536K
Total Submissions: 11593		Accepted: 3657

Description

Given n segments in the two dimensional space, write a program, which determines if there exists a line such that after projecting these segments on it, all projected segments have at least one point in common.

Input

Input begins with a number T showing the number of test cases and then, T test cases follow. Each test case begins with a line containing a positive integer n ≤ 100 showing the number of segments. After that, n lines containing four real numbers x₁ y₁ x₂ y₂ follow, in which (x₁, y₁) and (x₂, y₂) are the coordinates of the two endpoints for one of the segments.

Output

For each test case, your program must output "Yes!", if a line with desired property exists and must output "No!" otherwise. You must assume that two floating point numbers a and b are equal if |a - b| < 10^-8.

Sample Input

3
2
1.0 2.0 3.0 4.0
4.0 5.0 6.0 7.0
3
0.0 0.0 0.0 1.0
0.0 1.0 0.0 2.0
1.0 1.0 2.0 1.0
3
0.0 0.0 0.0 1.0
0.0 2.0 0.0 3.0
1.0 1.0 2.0 1.0

Sample Output

Yes!
Yes!
No!

Source

Amirkabir University of Technology Local Contest 2006

【题意】

　　给出n条线段两个端点的坐标，问所有线段投影到一条直线上，如果这些所有投影至少相交于一点就输出Yes！，否则输出No！。
【思路】

　　如果有存在这样的直线，过投影相交区域作直线的垂线，该垂线必定与每条线段相交，问题转化为问是否存在一条线和所有线段相交
　　若存在一条直线与所有线段相机相交，此时该直线必定经过这些线段的某两个端点，所以枚举任意两个端点即可

【代码】

#include<cmath>
#include<cstdio>
#include<cstring>
#define FOR(a,b,c) for(int a=(b);a<(c);a++)
using namespace std;

const int N = 100+10;
const double eps = 1e-8;

int dcmp(double x) {
    if(x<eps) return 0; else return x<0? -1:1;
}

struct Line {
    double x1,y1,x2,y2;
}L[N];

int n;

double cross(double x1,double y1,double x2,double y2,double x,double y) {
    return (x2-x1)*(y-y1)-(x-x1)*(y2-y1);
}
bool judge(double x1,double y1,double x2,double y2) {
    if(!dcmp(x2-x1) && !dcmp(y2-y1)) return 0;
    FOR(i,0,n) {
        if(cross(x1,y1,x2,y2,L[i].x1,L[i].y1)*
            cross(x1,y1,x2,y2,L[i].x2,L[i].y2)>eps) return 0;
    }
    return 1;
}

int main() {
    int T;
    scanf("%d",&T);
    while(T--) {
        scanf("%d",&n);
        FOR(i,0,n) 
            scanf("%lf%lf%lf%lf",&L[i].x1,&L[i].y1,&L[i].x2,&L[i].y2);
        if(n==1) { puts("Yes!"); continue; }
        bool ans=0;
        FOR(i,0,n) {
            FOR(j,i+1,n) {
                if(judge(L[i].x1,L[i].y1,L[j].x1,L[j].y1)) ans=1;
                if(judge(L[i].x1,L[i].y1,L[j].x2,L[j].y2)) ans=1;
                if(judge(L[i].x2,L[i].y2,L[j].x1,L[j].y1)) ans=1;
                if(judge(L[i].x2,L[i].y2,L[j].x2,L[j].y2)) ans=1;
                if(ans) break;
            }
            if(ans) break;
        }
        if(ans) puts("Yes!");
        else puts("No!");
    }
    return 0;
}