【POJ 3335】多边形的核

1.题目链接。题目大意：给定一个多边形，是不是在这个多边形的内部存在一个点使得站在这个点可以看到多边形中的任意一个点。

2.一个与多边形核有关的问题。多边形的核是一个点集，在这个点集内部，可以看到边界上的任意一个点。求解方法采用半平面交。

#include<iostream>
#include<cstdio>
#include<stdio.h>
#include<cstring>
using namespace std;
#pragma warning(disable:4996)
struct node
{
	double x;
	double y;
} p[110], temp[110], newp[110];//p是最开始的多边形的每个点，temp是中间过程中临时存的多边形的每个点，newp是切割后的多边形的每个点
int n, newn;//原来的点数，切割后的点数
//直线方程的三个系数
double a, b, c;
void getline(node x, node y)//求x与y两点确定的直线方程ax+by+c=0
{
	a = y.y - x.y;
	b = x.x - y.x;
	c = y.x*x.y - y.y*x.x;
}
node intersect(node x, node y)//求x与y点确定的直线与ax+by+c=0这条直线的交点
{
	double u = a * x.x + b * x.y + c;
	double v = a * y.x + b * y.y + c;
	node t;
	t.x = (x.x*v + y.x*u) / (u + v);//y.y-x.y=u+v;y.y-t.y=v;y.y-x.y=u;
	t.y = (x.y*v + y.y*u) / (u + v);
	return t;
}
void cut()
{
	int cutn = 0;
	for (int i = 1; i <= newn; i++)
	{
		if (a*newp[i].x + b * newp[i].y + c >= 0)//所有的点都大于0，说明所有的点都在这条直线的另一边，所以不用切
			temp[++cutn] = newp[i];
		else
		{
			if (a*newp[i - 1].x + b * newp[i - 1].y + c > 0)
				temp[++cutn] = intersect(newp[i - 1], newp[i]);//把新交点加入
			if (a*newp[i + 1].x + b * newp[i + 1].y + c > 0)
				temp[++cutn] = intersect(newp[i + 1], newp[i]);
		}
	}
	for (int i = 1; i <= cutn; i++)
		newp[i] = temp[i];
	newp[cutn + 1] = temp[1];//能够找出所有点的前驱和后继
	newp[0] = temp[cutn];
	newn = cutn;
}

void solve()
{
	for (int i = 1; i <= n; i++)
	{
		newp[i] = p[i];
	}
	p[n + 1] = p[1];
	newp[n + 1] = newp[1];
	newp[0] = newp[n];
	newn = n;
	for (int i = 1; i <= n; i++)
	{
		getline(p[i], p[i + 1]);//从头开始顺序遍历两个相邻点。
		cut();
	}

}
int main()
{
	int T;
	scanf("%d", &T);
	while (T--)
	{
		scanf("%d", &n);
		for (int i = 1; i <= n; i++)
			scanf("%lf %lf", &p[i].x, &p[i].y);
		solve();
		if (newn == 0) puts("NO");
		else
			puts("YES");
	}
	return 0;
}