关闭

[POJ1279]Art Gallery(半平面交)

373人阅读 评论(0) 收藏 举报
分类:

题目描述

传送门
题意:求多边形的核可行区域的面积

题解

用半平面交求出凸壳了之后求面积
都是基本操作

代码

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define N 10005

const double inf=1e9;
const double eps=1e-9;
int dcmp(double x)
{
    if (x<=eps&&x>=-eps) return 0;
    return (x>0)?1:-1;
}
struct Vector
{
    double x,y;
    Vector(double X=0,double Y=0)
    {
        x=X,y=Y;
    }
};
typedef Vector Point;
Vector operator + (Vector a,Vector b) {return Vector(a.x+b.x,a.y+b.y);}
Vector operator - (Vector a,Vector b) {return Vector(a.x-b.x,a.y-b.y);}
Vector operator * (Vector a,double p) {return Vector(a.x*p,a.y*p);}

int T,n,cnt,ncnt;
double x,y,ans;
Point p[N],poly[N],npoly[N];

double Cross(Vector A,Vector B)
{
    return A.x*B.y-A.y*B.x;
}
bool insLS(Point A,Point B,Point C,Point D)
{
    Vector v=B-A,w=C-A,u=D-A;
    return dcmp(Cross(v,w))!=dcmp(Cross(v,u));
}
Point GLI(Point P,Vector v,Point Q,Vector w)
{
    Vector u=P-Q;
    double t=Cross(w,u)/Cross(v,w);
    return P+v*t;
}
void init()
{
    cnt=0;
    poly[++cnt]=Point(inf,inf);
    poly[++cnt]=Point(inf,-inf);
    poly[++cnt]=Point(-inf,-inf);
    poly[++cnt]=Point(-inf,inf);
}
void halfp(Point A,Point B)
{
    ncnt=0;
    Point C,D;
    for (int i=1;i<=cnt;++i)
    {
        C=poly[i%cnt+1];
        D=poly[(i+1)%cnt+1];
        if (dcmp(Cross(B-A,C-A))<=0)
            npoly[++ncnt]=C;
        if (insLS(A,B,C,D))
            npoly[++ncnt]=GLI(A,B-A,C,D-C);
    }
    cnt=ncnt;
    for (int i=1;i<=cnt;++i)
        poly[i]=npoly[i];
}
double Area()
{
    double area=0;
    for (int i=2;i<cnt;++i)
        area+=Cross(poly[i]-poly[1],poly[i+1]-poly[1]);
    return fabs(area/2.0);
}
int main()
{
    scanf("%d",&T);
    while (T--)
    {
        scanf("%d",&n);
        for (int i=1;i<=n;++i)
        {
            scanf("%lf%lf",&x,&y);
            p[i]=Point(x,y);
        }
        init();
        for (int i=1;i<=n;++i)
            halfp(p[i%n+1],p[(i+1)%n+1]);
        ans=Area();
        printf("%.2lf\n",ans);
    }
}
0
0
查看评论
发表评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场

poj1279Art Gallery【半平面交求内核面积】

Language: Default Art Gallery Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 6027   Accep...
  • R1986799047
  • R1986799047
  • 2015-09-14 23:45
  • 618

半平面交初级

半平面:顾名思义,半平面就是平面的一半,那么半平面应该如何表示呢?我们知道一条直线可以将平面分成两半,如果直线的方程为ax+by+c=0ax+by+c = 0,那么半平面就可以表示为ax+by+c<0或者ax+by+c>0ax+by+c< 0 或者 ax+by+c>0,如果我...
  • huayunhualuo
  • huayunhualuo
  • 2016-04-15 20:15
  • 801

半平面交总结and模板

这两天刷了POJ上几道半平面交,对半平面交有了初步的体会,感觉半平面交还是个挺实用的知识点。 半平面交主要是看的ZZY的国家队论文,他提出的是一种O(n×log(n))的排序增量法。 附论文地址: 算法合集之《半平面交的新算法及其实用价值》。 POJ 3335 Ro...
  • xuechelingxiao
  • xuechelingxiao
  • 2014-11-06 16:44
  • 2597

计算几何学习之半平面交

首先解决问题:什么是半平面? 顾名思义,半平面就是指平面的一半,我们知道,一条直线可以将平面分为两个部分,那么这两个部分就叫做两个半平面。 然后,半平面怎么表示呢? 二维坐标系下,直线可以表示为ax + by + c = 0,那么两个半平面则可以表示为ax + by + c...
  • kep159
  • kep159
  • 2014-04-22 15:51
  • 988

POJ 2451 nlog(n)半平面交裸题。

Uyuw's Concert Time Limit: 6000MS   Memory Limit: 65536K Total Submissions: 6587   Accepted: 2611 Description...
  • u012358934
  • u012358934
  • 2014-05-04 18:44
  • 974

【半平面交】[ZJOI2008][HYSBZ\BZOJ1038]瞭望塔

题目链接分析看穿题目考虑组成村庄的每一条线段,显然,我们要在这条线段所在直线上方的半平面内才能看见它。所以,瞭望塔必须要在所有组成村庄的线段的所在直线的上方的半平面的交集内,才能 从瞭望塔的顶端可以看到H村的任意位置 所以,这道题就是求村庄的地面到这个半平面交的最短距离。实现方式很多同学看到半平...
  • outer_form
  • outer_form
  • 2016-02-05 17:37
  • 416

HDU 2297 半平面交

Run Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 640  ...
  • u012358934
  • u012358934
  • 2014-05-15 17:28
  • 1206

LA 2512 —— Art Gallery(半平面交)

题目:Art Gallery 题意:在一个建筑物li'mi
  • hongrock
  • hongrock
  • 2014-07-02 15:41
  • 517

[半平面交对偶凸包] BZOJ 1007 [HNOI2008]水平可见直线

半平面交对偶转凸包问题 神奇的新姿势 还有待研究#include<cstdio> #include<cstdlib> #include<algorithm> using namespace std; typedef long long ll;inline cha...
  • u014609452
  • u014609452
  • 2017-02-14 10:28
  • 192

POJ 1279 Art Gallery 半平面交求面积

和CQOI 2006 凸多边形就多组数据的区别。。 #include #include #include using namespace std; const int N = 1502; struct Point { double x, y; Point(){} Po...
  • huanghongxun
  • huanghongxun
  • 2016-01-10 15:38
  • 268
    个人资料
    • 访问:599408次
    • 积分:17464
    • 等级:
    • 排名:第647名
    • 原创:1151篇
    • 转载:7篇
    • 译文:0篇
    • 评论:168条
    联系方式
    QQ:1209628538 加好友请备注省份和ID
    E-mail:fiona_2000@126.com
    最新评论