半平面交学习笔记

最新推荐文章于 2021-11-29 23:22:09 发布
最新推荐文章于 2021-11-29 23:22:09 发布
阅读量266 收藏 2
点赞数
分类专栏： # 半平面交 计算几何
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ez_lcw/article/details/100024138
版权

半平面交学习笔记

半平面

半平面：一条直线把一个平面分成的两个平面，如图，直线 A B AB AB把平面分成左（上）半平面和右（下）半平面的两个平面。

半平面1

那我们如何判断是在哪一个平面呢？有两种判定方法：

  1. 对于直线 A B AB AB，我们可以通过 A A A、 B B B的坐标算出其解析式，不妨设为 a x + b y + c = 0 ax+by+c=0 ax+by+c=0，那么当 C C C在左平面时，将 C C C的坐标代入直线 A B AB AB的解析式，其值大于0，即 a × C x + b × C y + c &gt; 0 a\times C_x+b\times C_y+c&gt;0 a×Cx​+b×Cy​+c>0；当 D D D在右平面时，将 D D D的坐标也代入直线 A B AB AB的解析式，其值小于0，即 a × D x + b × D y + c &lt; 0 a\times D_x+b\times D_y+c&lt;0 a×Dx​+b×Dy​+c<0。

  2. （前置知识：矢量的叉乘，可见我的矢量&凸包学习笔记）对于直线 A B AB AB顺次的两个端点 A A A、 B B B，如果 c r o s s ( B − A , C − A ) &gt; 0 cross(B-A,C-A)&gt;0 cross(B−A,C−A)>0（这里的相减操作定义为 A − B = ( A x − B x , A y − B y ) A-B=(A_x-B_x,A_y-B_y) A−B=(Ax​−Bx​,Ay​−By​)，所以 B − A B-A B−A就是 B B B以 A A A为原点的坐标， C − A C-A C−A就是 C C C以 A A A为原点的坐标），那么根据右手法则， C C C在左（上）半平面；如果 c r o s s ( B − A , D − A ) &lt; 0 cross(B-A,D-A)&lt;0 cross(B−A,D−A)<0，那么根据右手法则， D D D在右（下）半平面。

半平面交

半平面交，顾名思义，就是许多半平面的交集，一般我们指的是左半平面。

我们先举道例题：

P4196 [CQOI2006]凸多边形

做法：首先我们注意到，某个凸多边形的每条边所在直线的半平面交就是这个凸多边形：

半平面交1

那么我们可以把题目给出的所有凸多边形变成对应边数的直线，然后求这些直线的半平面交。

半平面交怎么求？

首先，对于每一个凸多边形，按逆时针将每一条边建成直线。

然后对于每一条直线，以起点为原点，求出它们各自的极角，并将它们按极角排序。

排序完后，如果有两条直线极角相同，那么我们肯定只留最左侧的直线，我们就要去重。

然后将第 1 1 1、 2 2 2条直线放入双端队列 d e q u e deque deque， d e q u e deque deque里存的是当前半平面交的边，因为半平面交是一个凸多边形

然后每一次加入直线到 d e q u e deque deque之前：

  1. 计算队尾的两条直线（也就是最陡的两条）的交点，如果它在要加入的直线的右边，那么就弹出队尾。

  2. 计算队首的两条直线（也就是最缓的两条）的交点，如果它在要加入的直线的右边，那么就弹出队首。

然后将新直线加入队尾。

加入完所有直线后：

  1. 计算队尾的两条直线（也就是最陡的两条）的交点，如果它在队首的直线的右边，那么就弹出队尾。

  2. 计算队首的两条直线（也就是最缓的两条）的交点，如果它在队尾的直线的右边，那么就弹出队首。

然后算出每两条直线的交点，这就是半平面交——一个凸多边形的顶点。

最后直接算出其面积就好了。（算面积方法详见我的矢量&凸包学习笔记）

最后代码如下：

#include<bits/stdc++.h>

#define eps 1e-6
#define N 510

using namespace std;

struct Point
{
    double x,y;
    Point(){};
    Point(double a,double b)
    {
        x=a,y=b;
    }
}p[N],e[N];

struct Line
{
    Point s,t;
    double d;
    Line(){};
    Line(Point a,Point b)
    {
        s=a,t=b;
    }
}l[N],q[N];

Point operator + (Point a,Point b)
{
    return Point{a.x+b.x,a.y+b.y};
}

Point operator - (Point a,Point b)
{
    return Point{a.x-b.x,a.y-b.y};
}

Point operator * (Point a,double x)
{
    return Point{a.x*x,a.y*x};
}

int n,cnt,cnt1;
double ans;

double cross(Point a,Point b)
{
    return a.x*b.y-b.x*a.y;
}

double dis(Point a,Point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

int compare(double x,double y)
{
    if(fabs(x-y)<eps)return 0;
    return x-y<0?-1:1;
}

bool cmp(Line a,Line b)
{
    if(!compare(a.d,b.d))return cross(a.t-a.s,b.t-a.s)>0;
    return a.d<b.d;
}

bool onright(Point a,Line l)
{
    return compare(cross(l.t-l.s,a-l.s),0)>0?0:1;
}

Point intersection(Line a,Line b)
{
    Point u=a.s-b.s;
    Point v=a.t-a.s;
    Point w=b.t-b.s;
    double t=cross(w,u)/cross(v,w);
    return a.s+v*t;
}

double angle(Point a)
{
    return atan2(a.y,a.x);
}

void half()
{
    sort(l+1,l+cnt+1,cmp);
    cnt1=0;
    for(int i=1;i<cnt;i++)
    {
        if(!compare(l[i+1].d-l[i].d,0))continue;
        l[++cnt1]=l[i];
    }
    l[++cnt1]=l[cnt];
    cnt=cnt1;
    int L=1,R=0;
    q[++R]=l[1],q[++R]=l[2];
    for(int i=3;i<=cnt;i++)
    {
        while(L<R&&onright(intersection(q[R],q[R-1]),l[i]))R--;
        while(L<R&&onright(intersection(q[L],q[L+1]),l[i]))L++;
        q[++R]=l[i];
    }
    while(L<R&&onright(intersection(q[R],q[R-1]),q[L]))R--;
    while(L<R&&onright(intersection(q[L],q[L+1]),q[R]))L++;
    q[++R]=q[L];
    cnt1=0;
    for(int i=L;i<R;i++)e[++cnt1]=intersection(q[i],q[i+1]);
    e[++cnt1]=e[1];
}

int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        int m;
        scanf("%d",&m);
        Point last,st;
        for(int j=1;j<=m;j++)
        {
            double x,y;
            scanf("%lf%lf",&x,&y);
            if(j==1)
            {
                last=st=Point{x,y};
                continue;
            }
            l[++cnt]=Line{last,Point{x,y}};
            last=Point{x,y};
        }
        l[++cnt]=Line{last,st};
    }
    for(int i=1;i<=cnt;i++)l[i].d=angle(l[i].s-l[i].t);
    half();
    if(cnt1<3)
    {
        puts("0.000");
        return 0;
    }
    for(int i=1;i<cnt1;i++)ans+=cross(e[i],e[i+1]);
    printf("%.3lf\n",ans/2.0);
    return 0;
}

练习

poj3130 How I Mathematician Wonder What You Are!/poj3335 Rotating Scoreboard/poj1474 Video Surveillance 判断半平面交是否是平面而不是一条线、一个点或没有。看最后交点数是否大于2不就好了……

代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>

#define eps 1e-6
#define N 510

using namespace std;

struct Point
{
    double x,y;
    Point(){};
    Point(double a,double b)
    {
        x=a,y=b;
    }
}p[N],e[N];

struct Line
{
    Point s,t;
    double d;
    Line(){};
    Line(Point a,Point b)
    {
        s=a,t=b;
    }
}l[N],q[N];

Point operator + (Point a,Point b)
{
    return Point(a.x+b.x,a.y+b.y);
}

Point operator - (Point a,Point b)
{
    return Point(a.x-b.x,a.y-b.y);
}

Point operator * (Point a,double x)
{
    return Point(a.x*x,a.y*x);
}

int n,cnt,cnt1;
double ans;

double cross(Point a,Point b)
{
    return a.x*b.y-b.x*a.y;
}

double dis(Point a,Point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

int compare(double x,double y)
{
    if(fabs(x-y)<eps)return 0;
    return x-y<0?-1:1;
}

bool cmp(Line a,Line b)
{
    if(!compare(a.d,b.d))return cross(a.t-a.s,b.t-a.s)>0;
    return a.d<b.d;
}

bool onright(Point a,Line l)
{
    return compare(cross(l.t-l.s,a-l.s),0)>0?0:1;
}

Point intersection(Line a,Line b)
{
    Point u=a.s-b.s;
    Point v=a.t-a.s;
    Point w=b.t-b.s;
    double t=cross(w,u)/cross(v,w);
    return a.s+v*t;
}

double angle(Point a)
{
    return atan2(a.y,a.x);
}

void half()
{
    sort(l+1,l+cnt+1,cmp);
    cnt1=0;
    for(int i=1;i<cnt;i++)
    {
        if(!compare(l[i+1].d-l[i].d,0))continue;
        l[++cnt1]=l[i];
    }
    l[++cnt1]=l[cnt];
    cnt=cnt1;
    int L=1,R=0;
    q[++R]=l[1],q[++R]=l[2];
    for(int i=3;i<=cnt;i++)
    {
        while(L<R&&onright(intersection(q[R],q[R-1]),l[i]))R--;
        while(L<R&&onright(intersection(q[L],q[L+1]),l[i]))L++;
        q[++R]=l[i];
    }
    while(L<R&&onright(intersection(q[R],q[R-1]),q[L]))R--;
    while(L<R&&onright(intersection(q[L],q[L+1]),q[R]))L++;
    q[++R]=q[L];
    cnt1=0;
    for(int i=L;i<R;i++)e[++cnt1]=intersection(q[i],q[i+1]);
    e[++cnt1]=e[1];
}

int main()
{
	while(~scanf("%d",&n)&&n)
	{
		cnt=0;
		Point last,st;
	    for(int i=1;i<=n;i++)
	    {
	        double x,y;
	        scanf("%lf%lf",&x,&y);
	        if(i==1)
	        {
	        	last=st=Point(x,y);
	        	continue;
	        }
	        l[++cnt]=Line(last,Point(x,y));
	        last=Point(x,y);
	    }
	    l[++cnt]=Line(last,st);
	    for(int i=1;i<=cnt;i++)l[i].d=angle(l[i].s-l[i].t);
	    half();
	    ans=0;
	    if(cnt1<3)ans=0;
	    else for(int i=1;i<cnt1;i++)ans+=cross(e[i],e[i+1]);
	    if(!ans)puts("0");
	    else puts("1");
	}
    return 0;
}

poj1279 Art Gallery 模板……代码不放了。

poj3384 Feng Shui 将凸多边形的每条边内移 r r r，那么圆心就在这些直线的半平面交的端点上，最后暴力枚举或旋转卡壳即可。

代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>

#define eps 1e-6
#define N 110
#define INF 0x7fffffff

using namespace std;

struct Point
{
    double x,y;
    Point(){};
    Point(double a,double b)
    {
        x=a,y=b;
    }
}e[N],p[N],p1,p2;

Point operator + (Point a,Point b)
{
    return Point(a.x+b.x,a.y+b.y);
}

Point operator - (Point a,Point b)
{
    return Point(a.x-b.x,a.y-b.y);
}

Point operator * (Point a,double x)
{
    return Point(a.x*x,a.y*x);
}

double angle(Point a)
{
    return atan2(a.y,a.x);
}

struct Line
{
    Point s,t;
    double d;
    Line(){};
    Line(Point a,Point b)
    {
        s=a,t=b,d=angle(b-a);
    }
}l[N],q[N];

int n,cnt,cnt1;
double r,ans=-INF;

double cross(Point a,Point b)
{
    return a.x*b.y-b.x*a.y;
}

double dis(Point a,Point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

int compare(double x,double y)
{
    if(fabs(x-y)<eps)return 0;
    return x-y<0?-1:1;
}

bool cmp(Line a,Line b)
{
    if(!compare(a.d,b.d))return compare(cross(a.t-a.s,b.t-a.s),0)>0;
    return a.d<b.d;
}

bool onright(Point a,Line l)
{
    return compare(cross(l.t-l.s,a-l.s),0)>=0?0:1;
}

Point intersection(Line a,Line b)
{
    Point u=a.s-b.s;
    Point v=a.t-a.s;
    Point w=b.t-b.s;
    double t=cross(w,u)/cross(v,w);
    return a.s+v*t;
}

void half()
{
    sort(l+1,l+cnt+1,cmp);
    cnt1=0;
    for(int i=1;i<cnt;i++)
    {
        if(!compare(l[i+1].d-l[i].d,0))continue;
        l[++cnt1]=l[i];
    }
    l[++cnt1]=l[cnt];
    cnt=cnt1;
    int L=1,R=0;
    q[++R]=l[1],q[++R]=l[2];
    for(int i=3;i<=cnt;i++)
    {
        while(L<R&&onright(intersection(q[R],q[R-1]),l[i]))R--;
        while(L<R&&onright(intersection(q[L],q[L+1]),l[i]))L++;
        q[++R]=l[i];
    }
    while(L<R&&onright(intersection(q[R],q[R-1]),q[L]))R--;
    while(L<R&&onright(intersection(q[L],q[L+1]),q[R]))L++;
    q[++R]=q[L];
    cnt1=0;
    for(int i=L;i<R;i++)e[++cnt1]=intersection(q[i],q[i+1]);
}

int main()
{
    scanf("%d%lf",&n,&r);
    for(int i=1;i<=n;i++)scanf("%lf%lf",&p[i].x,&p[i].y);
    p[0]=p[n];
    for(int i=n;i>=1;i--)
    {
        double k=r/dis(p[i-1],p[i]);
        double x1=k*(p[i-1].y-p[i].y);
        double y1=k*(p[i].x-p[i-1].x);
        l[++cnt]=Line(Point(p[i].x-x1,p[i].y-y1),Point(p[i-1].x-x1,p[i-1].y-y1));
    }
    half();
    for(int i=1;i<=cnt1;i++)
    {
        for(int j=i;j<=cnt1;j++)
        {
            double d=dis(e[i],e[j]);
            if(d>ans)
            {
                ans=d;
                p1=e[i];
                p2=e[j];
            }
        }
    }
    printf("%.4lf %.4lf %.4lf %.4lf\n",p1.x,p1.y,p2.x,p2.y);
    return 0;
}
ez_lcw
关注
  • 0
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
半平面交 学习笔记
永爱米蒂
07-04 351
前置知识： 文章目录定义 定义
[计算几何]求半平面交的面积
lb2003
05-07 488
题面描述 传送门 思路 先看一下SCY的视频吧，讲得真的很清(luo)楚(suo)。 忽略掉一些奇奇怪怪的语言 首先我们先来了解一下SCY开头说的$$
开关电源中的“右半平面零点问题（RHPZ）”
qq_32419593的博客
11-29 1万+
title: 开关电源中的“右半平面零点问题（RHPZ） date: 2021-11-29 23:18:54 about: categories: 开关电源 tags: 开关电源 学习笔记 拓扑类 推挽拓扑 mathjax: true 开关电源中的“右半平面零点问题（RHPZ） 在小信号频率补偿中，极点和零点通常位于复数s平面的左半部分左半平面极点会使增益下降、相位滞后，而零点则相反，会使增益上升、相位超前。RHPZ的问题在于它的效果是让增益增加（类似于一个传统的零点），但是相位滞后。这个特性即使给.
求多边形相交面积（暴力 / 半平面交）
Khasehemwy的博客
04-28 1249
暴力： 貌似只适用于求两个多边形的相交面积，但是可以求任意多边形（凹凸）。 1.设两个多边形点集为P1，P2，相交区域点集为P。 2.遍历P1的点是否在P2内，若在内部，则加入点集P；对P2也做相同操作。 3.遍历P1的边和P2的边，求出两个多边形的交点，将交点加入P。 4.求点集P的面积。 半平面交： 这个只适用凸包，但是可以求任意个凸包最终相交的面积。 1.把每个凸包的点按照逆时针排序。（当然也看自己的算法是要什么顺序，大多数需要逆时针） 2.将每个凸包的边放入边集L。 3.求边集L的半平面交，再求半平
求半平面交的面积模板
OZY的博客
03-06 313
题意 在一个有限大（-10 0000&lt;=x,y&lt;=10 0000)的平面坐标系上有n个半平面（注意有限的），每个半平面给出一条有向线段（x1,y1)——&gt;(x2,y2）。 每个半平面的有效区域都是左侧。求这n个半平面的交的面积。 题解 半平面交的模板 以前偷懒，学了个(n2)(n2)(n^2)的，但是已经忘了啊 于是今天趁没网学了个的nlogn的nlogn的nlo...
csu 1812（半平面交求面积）
Infinity_Izayoi的博客
05-06 608
半平面交
计算几何半平面交学习笔记
最新发布
12-25
计算几何学习笔记
机器学习笔记
08-10
《机器学习笔记》是一份深度探讨机器学习领域的宝贵资料，其详尽的讲解与实例解析使得读者能够系统地掌握这一关键技术。在这个数字化时代，大数据的崛起使得机器学习成为理解和预测复杂现象的重要工具，也是推动各行...
notes-machine-learning-master_机器学习笔记_机器学习_
10-03
【机器学习笔记_机器学习】概述 机器学习是人工智能的一个重要分支，它研究计算机如何从经验中学习，以便在没有明确编程的情况下改进其性能。这个领域由统计学、计算机科学和神经科学等多个学科交叉而成，旨在使...
机器学习与数据挖掘学习笔记.zip
04-27
《机器学习与数据挖掘学习笔记》是一份综合性的学习资料，涵盖了这两个领域的重要概念、算法和技术。这份笔记的目的是为了帮助读者深入理解机器学习和数据挖掘的基础知识，并提供实际操作的指导。 首先，我们来探讨...
斯坦福大学机器学习课程个人学习笔记（上）
09-10
根据是否需要人类标注数据，机器学习可以分为监督学习、无监督学习和半监督学习。监督学习是最常见的形式，包括回归和分类问题。在回归中，目标是预测连续变量，如房价；而在分类中，目标是将数据分配到离散的类别，...
【详解】半平面交算法入门详解(计算几何)
热门推荐
Winsoul Blog
10-30 3万+
半平面交 简介 博客背景 笔者在学习半平面交时，网上找入门博客资源甚少，且大部分难以理解，故在稍稍入门了半平面交后，写此博客，希望能对大家有所帮助。若有错误，麻烦指出。 半平面交是什么？ 我们知道一条直线可以把平面分为两部分，其中一半的平面就叫半平面。 那半平面交，就是多个半平面的相交部分。我们在学习线性规划时就有用过。 半平面交有什么用？ 1.求解一个区域，可以看到给定图形的各个角落。(多边形的...
【计算几何】求半平面交的面积
curry的博客
05-07 271
【题意】 在一个有限大（-10 0000<=x,y<=10 0000)的平面坐标系上有n个半平面（注意有限的），每个半平面给出一条有向线段（x1,y1)——>(x2,y2）。 每个半平面的有效区域都是左侧。求这n个半平面的交的面积。 【输入文件】 第一行一个整数n（1 <= n <= 2 0000） 下来n个半平面。每个半平面一行四个整数x1,y1,x2,y2。（-1...
计算几何之半平面交、多边形的核
ACMer_quack&cheer
04-01 2961
半平面指一个平面的一半。比如在二维空间中，一条直线可以把整个平面分为两部分，左边是一个半平面右边是一个半平面。而半平面交就是一堆半平面的交集，可以想象成线性规划方程组A0x+B0y+C0>=0,A1x+B1y+C1<=0......A_0x+B_0y+C_0>=0,A_1x+B_1y+C_1<=0......的解集 如果一根直线是有向的，我们可以定义它的左边是它对应的半平面（也可以是右边），那么半
半平面交初级
绝尘
04-15 2004
半平面：顾名思义，半平面就是平面的一半，那么半平面应该如何表示呢？我们知道一条直线可以将平面分成两半，如果直线的方程为ax+by+c=0ax+by+c = 0,那么半平面就可以表示为ax+by+c<0或者ax+by+c>0ax+by+c< 0 或者 ax+by+c>0,如果我们这样定义半平面，那么半平面交也就很好解释了。 半平面交：每一个半平面就是一个不等方程，那么半平面交就是满足这些半平面不等方程
计算几何学习之半平面交
04-22 1546
首先解决问题：什么是半平面？ 顾名思义，半平面就是指平面的一半，我们知道，一条直线可以将平面分为两个部分，那么这两个部分就叫做两个半平面。 然后，半平面怎么表示呢？ 二维坐标系下，直线可以表示为ax + by + c = 0，那么两个半平面则可以表示为ax + by + c >= 0 和ax + by + c 还有，半平面的交是神马玩意？ 其实就是一个方程组，让你画出满足若干个式子的
poj 1279 求半平面交的 面积(推荐)
weixin_34074740的博客
04-29 81
poj 1279 求半平面交的 面积 题目来源： http://poj.org/problem?id=1279 分析： 求半平面交的 面积 代码如下： const double EPS = 1e-8 ; const int Max_N = 1505 ; double add(double a, double b){ return (fabs(a + b) <...
半平面交总结and模板
GLSilence的专栏
11-06 4706
这两天刷了POJ上几道半平面交，对半平面交有了初步的体会，感觉半平面交还是个挺实用的知识点。 半平面交主要是看的ZZY的国家队论文，他提出的是一种O(n×log(n))的排序增量法。 附论文地址： 算法合集之《半平面交的新算法及其实用价值》。 POJ 3335 Rotating Scoreboard 题目大意： World finals 要开始了，比赛场
poj1279——求半平面交后凸多边形的面积
mumei的博客
09-17 191
题目链接：http://poj.org/problem?id=1279 The art galleries of the new and very futuristic building of the Center for Balkan Cooperation have the form of polygons (not necessarily convex). When a big exhib...
半平面交 python
09-08
半平面交（Half-plane intersection）是计算凸多边形交集的一种常见算法。在 Python 中，可以使用库如 `shapely` 来实现半平面交算法。 首先，确保已经安装了 `shapely` 库。可以使用以下命令进行安装： ``` pip ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值