【POJ3335】博物馆守卫Rotating Scoreboard 半平面交

题目描述

  现在有一个博物馆,俯瞰图是一个多边形。现在想要在博物馆内部设立一个守卫,但要求他可以看见博物馆的每一个角落。现在想知道是否存在这样一个点,使得安置其上的守卫可以完成任务。

数据范围

点数小于100

样例输入

2 4 0 0 0 1 1 1 1 0 8 0 0 0 2 1 2 1 1 2 1 2 2 3 2 3 0

样例输出

YES NO

解题思路

和POJ1474的监控摄像头几乎完全没有区别

代码

#include <algorithm> #include <iostream> #include <cstring> #include <cstdlib> #include <cstdio> #include <cmath> #include <ctime> #define Maxn 233333 using namespace std; int n,nn,m; struct Point{ double x,y; }a[Maxn],b[Maxn],c[Maxn]; void Init(){ scanf("%d",&m); for(int i=1;i<=m;i++) scanf("%lf%lf",&b[i].x,&b[i].y); n=m; memcpy(a,b,sizeof(b)); } double Cross(Point a,Point b,Point c){ return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y); } void AddCross(Point a,Point b,Point C,Point d){ double c1=Cross(a,d,b),c2=Cross(a,C,b); c[++nn]=(Point){(c1*C.x-c2*d.x)/(c1-c2),(c1*C.y-c2*d.y)/(c1-c2)}; } void Cut(Point A,Point B){ nn=0; a[n+1]=a[1]; for(int i=1;i<=n;i++){ if(Cross(a[i],B,A)>=0){ c[++nn]=a[i]; if(Cross(a[i+1],B,A)<0)AddCross(A,B,a[i+1],a[i]); }else if(Cross(a[i+1],B,A)>0)AddCross(A,B,a[i+1],a[i]); } for(int i=1;i<=nn;i++)a[i]=c[i]; n=nn; } int main(){ int Case=0; scanf("%d",&Case); while(Case--){ Init(); b[m+1]=b[1]; for(int i=1;i<=m;i++) Cut(b[i],b[i+1]); if(n)cout<<"YES\n"; else cout<<"NO\n"; } }