poj 3130

题目概述

按逆时针序给定一个简单N边形所有顶点的坐标x,y，问其是否存在核

时限

5000ms/15000ms

输入

第一行正整数N，其后N行，每行两个整数x,y，输入到N=0结束

限制

4<=N<=50;0<=x,y<=10000

输出

每行一个数，如果存在核，为1，否则为0

样例输入

6
66 13
96 61
76 98
13 94
4 0
45 68
8
27 21
55 14
93 12
56 95
15 48
38 46
51 65
64 31
0

样例输出

1
0

讨论

计算几何，半平面交，和poj 3335几乎一致，除了输入输入的细微差异以及给出顺序相反，都是一样的，不再过多解释
需要说明一点，从一般情况而言，由于半平面交需要求新的交点，而使得点坐标不一定为整数，因而点的结构采用两个整数实际上是不够严谨的行为，应当采用double，但由于是求存在性，不会过于纠结小数部分（否则可能会由于浮点误差而变成坑人数据），因而采用整型也能蒙混过关

题解状态

152K，0MS，C++，2123B

题解代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define INF 0x3f3f3f3f  
#define MAXN 53
#define memset0(a) memset(a,0,sizeof(a))
#define EPS 1e-6

struct Pt//point 点的结构
{
    int x, y;//这里不应当使用整型
}pts[MAXN], pt1[MAXN], pt2[MAXN];//点的原始数据 以及半平面交的两个辅助数组
int hpt1, hpt2;//halfplane_point 配合辅助数组 记录点的数量
int N;//点总数
int xp(int x1, int y1, int x2, int y2, int x3, int y3)
{
    return (x1 - x2)*(y3 - y2) - (y1 - y2)*(x3 - x2);
}
Pt point_of_intersection(double x1, double y1, double x2, double y2, double x3, double y3, double x4, double y4)
{
    Pt poi;
    if (abs(x1 - x2) < EPS) {
        double k3 = (y4 - y3) / (x4 - x3);
        poi.x = x1, poi.y = y3 - k3*x3 + k3*x1;
    }
    else if (abs(x3 - x4) < EPS) {
        double k1 = (y2 - y1) / (x2 - x1);
        poi.x = x3, poi.y = y1 - k1*x1 + k1*x3;
    }
    else {
        double k1 = (y2 - y1) / (x2 - x1);
        double k3 = (y4 - y3) / (x4 - x3);
        poi.x = (y1 - y3 - k1*x1 + k3*x3) / (k3 - k1);
        poi.y = (k1*k3*(x1 - x3) - k3*y1 + k1*y3) / (k1 - k3);
    }
    return poi;
}
bool fun()
{
    for (int p = 0; p < N; p++) {
        scanf("%d%d", &pts[p].x, &pts[p].y);//input
        pt1[p] = pts[p];//顺手初始化点集1
    }
    hpt1 = N;//初始有N个点
    for (int p = 0; p < N; p++) {//枚举相邻的两个点 做判定直线
        hpt2 = 0;//点集2初始没有点
        for (int i = 0; i < hpt1; i++) {//枚举点集1中的点
            if (xp(pt1[i].x, pt1[i].y, pts[p].x, pts[p].y, pts[(p + 1) % N].x, pts[(p + 1) % N].y) <= 0)//在半平面内
                pt2[hpt2++] = pt1[i];
            else {//不在
                if (xp(pt1[(i - 1 + hpt1) % hpt1].x, pt1[(i - 1 + hpt1) % hpt1].y, pts[p].x, pts[p].y, pts[(p + 1) % N].x, pts[(p + 1) % N].y) < 0)//若其与前一点所成直线与判定线交点在半平面内
                    pt2[hpt2++] = point_of_intersection(pt1[i].x, pt1[i].y, pt1[(i - 1 + hpt1) % hpt1].x, pt1[(i - 1 + hpt1) % hpt1].y, pts[p].x, pts[p].y, pts[(p + 1) % N].x, pts[(p + 1) % N].y);
                if (xp(pt1[(i + 1) % hpt1].x, pt1[(i + 1) % hpt1].y, pts[p].x, pts[p].y, pts[(p + 1) % N].x, pts[(p + 1) % N].y) < 0)//若其与后一点所成直线与判定线交点在半平面内
                    pt2[hpt2++] = point_of_intersection(pt1[i].x, pt1[i].y, pt1[(i + 1) % hpt1].x, pt1[(i + 1) % hpt1].y, pts[p].x, pts[p].y, pts[(p + 1) % N].x, pts[(p + 1) % N].y);
            }
        }
        for (int p = 0; p < hpt2; p++)
            pt1[p] = pt2[p];
        hpt1 = hpt2;//将点集2的信息赋给点集1
    }
    return hpt1;//返回点集1中的点数 由于是布尔型 会进行转化
}
int main(void)
{
    //freopen("vs_cin.txt", "r", stdin);
    //freopen("vs_cout.txt", "w", stdout);

    while (~scanf("%d", &N) && N) {//input
        printf("%d\n", fun());//output
    }
}

EOF