題目大意:按順序給出一些點。可以構成一個多邊形,問多邊形是否有核。
解題思路:
半平面交。模版題。
用的是ZZY的O(Nlogn)的算法。
算法大體思路是這樣的:
step1. 將所有半平面按極角排序,對於極角相同的,選擇性的保留一個。 O(nlogn)
step2. 使用一個雙端隊列(deque),加入最開始2個半平面。
step3. 每次考慮一個新的半平面:
a.while deque頂端的兩個半平面的交點在當前半平面外:刪除deque頂端的半平面
b.while deque底部的兩個半平面的交點在當前半平面外:刪除deque底部的半平面
c.將新半平面加入deque頂端
step4.刪除兩端多餘的半平面。
具體方法是:
a.while deque頂端的兩個半平面的交點在底部半平面外:刪除deque頂端的半平面
b.while deque底部的兩個半平面的交點在頂端半平面外:刪除deque底部的半平面
重複a,b直到不能刪除爲止。
step5:計算出deque頂端和底部的交點即可。
下面是代碼:
#include <set>
#include <map>
#include <queue>
#include <math.h>
#include <vector>
#include <string>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
#define eps 1e-8
#define pi acos(-1.0)
#define inf 107374182
#define inf64 1152921504606846976
#define lc l,m,tr<<1
#define rc m + 1,r,tr<<1|1
#define zero(a) fabs(a)<eps
#define iabs(x) ((x) > 0 ? (x) : -(x))
#define clear1(A, X, SIZE) memset(A, X, sizeof(A[0]) * (SIZE))
#define clearall(A, X) memset(A, X, sizeof(A))
#define memcopy1(A , X, SIZE) memcpy(A , X ,sizeof(X[0])*(SIZE))
#define memcopyall(A, X) memcpy(A , X ,sizeof(X))
#define max( x, y ) ( ((x) > (y)) ? (x) : (y) )
#define min( x, y ) ( ((x) < (y)) ? (x) : (y) )
using namespace std;
const int maxn = 55;
int dq[maxn], top, bot, pn, order[maxn], ln;
struct Point
{
double x, y;
} p[maxn];
struct Line
{
Point a, b;
double angle;
} l[maxn];
int dblcmp(double k)
{
if(fabs(k) < eps) return 0;
return k > 0 ? 1 : -1;
}
double multi(Point p0, Point p1, Point p2)
{
return (p1.x-p0.x)*(p2.y-p0.y) - (p1.y-p0.y)*(p2.x-p0.x);
}
bool cmp(int u, int v)
{
int d = dblcmp(l[u].angle-l[v].angle);
if(!d) return dblcmp(multi(l[u].a, l[v].a, l[v].b)) > 0;
//大於0取向量左半部分爲半平面,小於0,取右半部分
return d < 0;
}
void getIntersect(Line l1, Line l2, Point& p)
{
double dot1, dot2;
dot1 = multi(l2.a, l1.b, l1.a);
dot2 = multi(l1.b, l2.b, l1.a);
p.x = (l2.a.x * dot2 + l2.b.x * dot1) / (dot2 + dot1);
p.y = (l2.a.y * dot2 + l2.b.y * dot1) / (dot2 + dot1);
}
bool judge(Line l0, Line l1, Line l2)
{
Point p;
getIntersect(l1, l2, p);
return dblcmp(multi(p, l0.a, l0.b)) < 0;
//大於小於符號與上面cmp()中註釋處相反
}
void addLine(double x1, double y1, double x2, double y2)
{
l[ln].a.x = x1;
l[ln].a.y = y1;
l[ln].b.x = x2;
l[ln].b.y = y2;
l[ln].angle = atan2(y2-y1, x2-x1);
order[ln] = ln;
ln++;
}
void halfPlaneIntersection()
{
int i, j;
sort(order, order+ln, cmp);
for(i = 1, j = 0; i < ln; i++)
if(dblcmp(l[order[i]].angle-l[order[j]].angle) > 0)
order[++j] = order[i];
ln = j + 1;
dq[0] = order[0];
dq[1] = order[1];
bot = 0;
top = 1;
for(i = 2; i < ln; i++)
{
while(bot < top && judge(l[order[i]], l[dq[top-1]], l[dq[top]]))
top--;
while(bot < top && judge(l[order[i]], l[dq[bot+1]], l[dq[bot]]))
bot++;
dq[++top] = order[i];
}
while(bot < top && judge(l[dq[bot]], l[dq[top-1]], l[dq[top]])) top--;
while(bot < top && judge(l[dq[top]], l[dq[bot+1]], l[dq[bot]])) bot++;
}
bool isThereACore()
{
if(top-bot > 1) return true;
return false;
}
int main()
{
int i;
while(scanf ("%d", &pn) && pn)
{
for(i = 0; i < pn; i++)
scanf ("%lf%lf", &p[i].x, &p[i].y);
for(ln = i = 0; i < pn-1; i++)
addLine(p[i].x, p[i].y, p[i+1].x, p[i+1].y);
addLine(p[i].x, p[i].y, p[0].x, p[0].y);
halfPlaneIntersection();
if(isThereACore()) printf ("1\n");
else printf ("0\n");
}
return 0;
}