poj 1279

最新推荐文章于 2020-01-20 11:50:42 发布

liurui39660

最新推荐文章于 2020-01-20 11:50:42 发布

阅读量372

点赞数

分类专栏：计算几何.半平面交 poj

本文链接：https://blog.csdn.net/liurui39660/article/details/52067395

版权

poj 同时被 2 个专栏收录

114 篇文章 0 订阅

订阅专栏

计算几何.半平面交

9 篇文章 0 订阅

订阅专栏

题目概述

按顺时针顺序给定一个简单N边形的所有顶点坐标x,y，求其核的可行域的面积

时限

1000ms/3000ms

输入

第一行正整数times，其后有times组数据，每组数据第一行正整数N，其后N行，每行两个整数x,y

限制

5<=N<=1500;x,y在16位整型范围内

输出

每行一个保留两位小数的浮点数，为所求可行域的面积

样例输入

1
7
0 0
4 4
4 7
9 7
13 -1
8 -6
4 -4

样例输出

80.00

讨论

计算几何，半平面交，直到这里WA了额才意识到点不该用整型存，至于讨论版所说的什么要加EPS之类的，并没有这种感觉，基本思路没有太大变化，和前面三个题poj 1474,poj 3130,poj 3335都是差不多的，只是多了个求面积的过程，这个也没有必要再展开说，很熟练了已经
题目没说默认给的是顺时针，如果是乱序的，这个题恐怕就复杂了，写代码的时候是猜的顺时针，所幸没猜错
1500组数据0ms就跑完，数据规模有点夸大了吧，再怎么说这也是平方级算法，给个0ms真的好？

题解状态

216K，0MS，C++，2346B

题解代码

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

struct Pt//点的结构
{
    double x, y;//最后才意识到应该用浮点数
}pts[MAXN], pt1[MAXN], pt2[MAXN];
int hpt1, hpt2;
int N;//点总数
double xp(double x1, double y1, double x2, double y2, double x3, double 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;
}
double fun()
{
    for (int p = 0; p < N; p++) {
        scanf("%lf%lf", &pts[p].x, &pts[p].y);//input
        pt1[p] = pts[p];
    }
    hpt1 = N;
    for (int p = 0; p < N; p++) {
        hpt2 = 0;
        for (int i = 0; i < hpt1; i++) {
            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;
    }
    double A = 0;//area 初始化面积为0
    for (int p = 0; p < hpt1 - 1; p++)
        A += xp(pt1[(p + 1) % hpt1].x, pt1[(p + 1) % hpt1].y, pt1[0].x, pt1[0].y, pt1[p].x, pt1[p].y);
    return A / 2.0;//最后统一除以2返回
}
int main(void)
{
    freopen("vs_cin.txt", "r", stdin);
    //freopen("vs_cout.txt", "w", stdout);

    int times;//数据组数
    scanf("%d", &times);//input
    while (times--) {
        scanf("%d", &N);//input
        printf("%.2lf\n", fun());//output
    }
}