POJ 3304 Segments

//按题目意思，所有线段按某方向进行投影，有一个公共点，即为找一条直线过所有线段//我们可以先找到一条过所有线段的直线，再将该直线进行偏移直到该直线经过任两个端点，则问题转化为对任两个端点形成的直线是否过所有线段 //像这种先找到我们需要的答案，再将答案进行偏移的思想是计算几何中一个常见的思想。 #include<stdio.h>#include<math.h>const double eps=1e-8;struct Point {double x,y;Point(){}Point(double xx,double yy){x=xx;y=yy;}}p[205];int T,n;int cross(Point a,Point b,Point c,Point d){double tmp1=(c.x-a.x)*(b.y-a.y)-(c.y-a.y)*(b.x-a.x);double tmp2=(d.x-a.x)*(b.y-a.y)-(d.y-a.y)*(b.x-a.x);if(fabs(tmp1*tmp2)<eps) return 0;else if(tmp1*tmp2>eps) return -1;else return 1;}int judge(Point a,Point b){if(fabs(a.x-b.x)<eps&& fabs(a.y-b.y)<eps) return 0;//这个一定要加 for(int i=0;i<n;i++){if(cross(a,b,p[i*2],p[i*2+1])==-1) return 0;}return 1;}int main() {#ifndef ONLINE_JUDGEfreopen("in.txt","r",stdin);#endifscanf("%d",&T);while(T--){scanf("%d",&n);for(int i=0;i<n;i++){scanf("%lf%lf%lf%lf",&p[i*2].x,&p[i*2].y,&p[i*2+1].x,&p[i*2+1].y);}int flag=0;for(int i=0;i<2*n;i++){for(int j=i+1;j<2*n;j++){if(judge(p[i],p[j])) flag=1;}}if(flag) printf("Yes!\n");else printf("No!\n");}}