17.8.19 校内赛 解题报告【求线段交点+凸包+求多边形面积】【判定点是否在多边形内】【二分答案+半平面交】

最新推荐文章于 2020-01-28 01:23:15 发布
最新推荐文章于 2020-01-28 01:23:15 发布
阅读量348 收藏
点赞数 1
分类专栏： 考试总结 文章标签： 计算几何
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/theljt/article/details/77416150
版权

Problem 1. land

解题报告
显然，这道题就是线段求交点+凸包+求多边形面积
问题就是我们只有直线求交点的模板和判线段相交的模板，没有直接求线段相交的。我是这样解决的：先判断两线段是否相交，然后再把线段存储为直线，求出两者的交点。
代码如下：

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#define Point Vector
#define NAME "land"
using namespace std;
const double eps=1e-8;
const int N=1000;
int sign(double a){return a>eps?1:(a<-eps?-1:0);}
struct Vector
{
    double x,y;
    double len(){return sqrt(x*x+y*y);}
    double ang(){return atan2(y,x);}
    Vector(double a=.0,double b=.0):x(a),y(b){}
    Vector operator+(const Vector &s){return Vector(x+s.x,y+s.y);}
    Vector operator-(const Vector &s){return Vector(x-s.x,y-s.y);}
    Vector operator*(double s){return Vector(x*s,y*s);}
    Vector operator/(double s){return Vector(x/s,y/s);}
};
struct Seg
{
    Point A,B;
    Seg(){}
    Seg(Point a,Point b):A(a),B(b){}
}Segs[N+5];
struct Line
{
    Point p;
    Vector u;
    double ang;
    Line(){}
    Line(Point p,Vector v):p(p),u(v),ang(u.ang()){}
    bool operator<(const Line &s){return ang<s.ang;}
    bool operator==(const Line &s){return sign(ang-s.ang)==0;}
};
bool cmp(Point r,Point s)
{
    if(r.x==s.x)return r.y<s.y;
    return r.x<s.x;
}
double dot(const Vector &a,const Vector &b){return a.x*b.x+a.y*b.y;}//向量点积 
double cross(const Vector &a,const Vector &b){return a.x*b.y-a.y*b.x;}//向量叉积
bool onleft(Point a,Point b,Point p)
{
    return sign(cross(b-a,p-a))>0;
}
bool onleft(Line l,Point p)
{
    return sign(cross(l.u,p-l.p))>0;
}
vector<Point>convex(vector<Point>pts)//求凸包 
{
    int n=(int)pts.size(),m=0;
    vector<Point>cvx;
    sort(pts.begin(),pts.end(),cmp);
    for(int i=0;i<n;i++)
    {
        while(m>1&&onleft(cvx[m-2],cvx[m-1],pts[i])){cvx.pop_back();m--;}
        cvx.push_back(pts[i]);m++;
    }
    int k=m;
    for(int i=n-2;i>=0;i--)
    {
        while(m>k&&onleft(cvx[m-2],cvx[m-1],pts[i])){cvx.pop_back();m--;}
        cvx.push_back(pts[i]);m++;
    }
    if(n>1){m--;cvx.pop_back();}
    return cvx;
}
double area(vector<Point>poly)//求多边形的面积 
{
    double rt=0.0;
    int n=(int)poly.size();
    for(int i=0;i<n;i++)rt+=cross(poly[i],poly[(i+1)%n]);
    return fabs(rt/2.0);
}
Point line_intersect(Point P,Vector u,Point Q,Vector v) 
{
    double t=cross(Q-P,v)/cross(u,v);
    return P+u*t;
}
Point line_intersect(const Line &a,const Line &b)//求出两直线交点 
{
    return line_intersect(a.p,a.u,b.p,b.u);
}
bool seg_intersect(Point A,Point B,Point C,Point D)//判定两线段（AB,CD）是否相交 
{
    double c1=cross(B-A,C-A),c2=cross(B-A,D-A);
    double c3=cross(D-C,A-C),c4=cross(D-C,B-C);
    return sign(c1)*sign(c2)<0&&sign(c3)*sign(c4)<0;
}
int n;
vector<Point>pts,poly; 
int main()
{
    freopen(NAME".in","r",stdin);
    freopen(NAME".out","w",stdout);
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        double x1,y1,x2,y2;
        scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
        Segs[i]=Seg(Point(x1,y1),Point(x2,y2));
    }
    for(int i=1;i<=n-1;i++)
    for(int j=i+1;j<=n;j++)
    {
        if(seg_intersect(Segs[i].A,Segs[i].B,Segs[j].A,Segs[j].B))
        pts.push_back(line_intersect(Line((Segs[i].A),Segs[i].B-Segs[i].A),Line((Segs[j].A),Segs[j].B-Segs[j].A)));
    }
    poly=convex(pts);
    double ans=area(poly);
    printf("%.2lf",ans);
    return 0;
}

Problem 2. escape

解题报告
直接套模板上：

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#define Point Vector
#define NAME "escape"
using namespace std;
const double eps=1e-8;
int sign(double a){return a>eps?1:(a<-eps?-1:0);}
struct Vector
{
    double x,y;
    double len(){return sqrt(x*x+y*y);}
    double ang(){return atan2(y,x);}
    Vector(double a=.0,double b=.0):x(a),y(b){}
    Vector operator+(const Vector &s){return Vector(x+s.x,y+s.y);}
    Vector operator-(const Vector &s){return Vector(x-s.x,y-s.y);}
    Vector operator*(double s){return Vector(x*s,y*s);}
    Vector operator/(double s){return Vector(x/s,y/s);}
};
double dot(const Vector &a,const Vector &b){return a.x*b.x+a.y*b.y;}//向量点积 
double cross(const Vector &a,const Vector &b){return a.x*b.y-a.y*b.x;}//向量叉积
bool point_on_line(Point p,Point A,Point B)//判定点是否在线段上 
{
    return sign(cross(B-A,p-A))==0&&sign(dot(A-p,B-p))<=0;
}
bool in_polygon(Point p,vector<Point>poly)//判定点是否在多边形内 
{
    int n=(int)poly.size();
    int counter=0;
    for(int i=0;i<n;++i)
    {
        Point a=poly[i],b=poly[(i+1)%n];
        if(point_on_line(p,a,b))return true; // bounded included
        int x=sign(cross(p-a,b-a));
        int y=sign(a.y-p.y);
        int z=sign(b.y-p.y);
        if(x>0&&y<=0&&z>0)counter++;
        if(x<0&&z<=0&&y>0)counter--;
    }
    return counter!=0;
}
int n,q;
vector<Point>poly;
int main()
{
    freopen(NAME".in","r",stdin);
    freopen(NAME".out","w",stdout);
    scanf("%d%d",&n,&q);
    for(int i=1;i<=n;i++)
    {
        double x,y;
        scanf("%lf%lf",&x,&y);
        poly.push_back(Point(x,y));
    }
    while(q--)
    {
        double x,y;
        scanf("%lf%lf",&x,&y);
        if(in_polygon(Point(x,y),poly))printf("YES\n");
        else printf("NO\n");
    }
    return 0;
}

Problem 3. island

解题报告
这是蓝书上的例题，可以看我的这篇博客
代码如下：

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#define NAME "island"
#define Point Vector
using namespace std;
const double eps=1e-8;
const int N=1e5;
int sign(double a){return a>eps?1:(a<-eps?-1:0);}
struct Vector
{
    double x,y;
    double len(){return sqrt(x*x+y*y);}
    double ang(){return atan2(y,x);}
    Vector(double a=.0,double b=.0):x(a),y(b){}
    Vector operator+(const Vector &s){return Vector(x+s.x,y+s.y);}
    Vector operator-(const Vector &s){return Vector(x-s.x,y-s.y);}
    Vector operator*(double s){return Vector(x*s,y*s);}
    Vector operator/(double s){return Vector(x/s,y/s);}
};
Vector Normal(Vector a){double L=a.len();return Vector(-a.y/L,a.x/L);}
struct Line
{
    Point p;
    Vector u;
    double ang;
    Line(){}
    Line(Point p,Vector v):p(p),u(v),ang(u.ang()){}
    bool operator<(const Line &s){return ang<s.ang;}
    bool operator==(const Line &s){return sign(ang-s.ang)==0;}
    bool operator<(const Line&L)const{return ang<L.ang;}
};
double dot(const Vector &a,const Vector &b){return a.x*b.x+a.y*b.y;}//向量点积 
double cross(const Vector &a,const Vector &b){return a.x*b.y-a.y*b.x;}//向量叉积
Point line_intersect(Point P,Vector u,Point Q,Vector v) 
{
    double t=cross(Q-P,v)/cross(u,v);
    return P+u*t;
}
Point line_intersect(const Line &a,const Line &b)//求出两直线交点 
{
    return line_intersect(a.p,a.u,b.p,b.u);
}
bool seg_intersect(Point A,Point B,Point C,Point D)//判定两线段（AB,CD）是否相交 
{
    double c1=cross(B-A,C-A),c2=cross(B-A,D-A);
    double c3=cross(D-C,A-C),c4=cross(D-C,B-C);
    return sign(c1)*sign(c2)<0&&sign(c3)*sign(c4)<0;
}
bool cmp1(Line r,Line s)
{
    return r<s;
}
bool onleft(Point a,Point b,Point p)
{
    return sign(cross(b-a,p-a))>0;
}
bool onleft(Line l,Point p)
{
    return sign(cross(l.u,p-l.p))>0;
}
vector<Point>half_plane_intersect(vector<Line>L)
{
    int n=L.size();
    sort(L.begin(),L.end());
    int first,last;
    vector<Point> p(n);
    vector<Line> q(n);
    vector<Point> ans;
    q[first=last=0]=L[0];
    for(int i=1;i<n;i++)
    {
        while(first<last&&!onleft(L[i],p[last-1]))last--;
        while(first<last&&!onleft(L[i],p[first]))first++;
        q[++last]=L[i];
        if(fabs(cross(q[last].u,q[last-1].u))<eps)
        {
            last--;
            if(onleft(q[last],L[i].p))q[last]=L[i];
        }
        if(first<last)p[last-1]=line_intersect(q[last-1],q[last]);
    }
    while(first<last&&!onleft(q[first],p[last-1]))last--;
    if(last-first<=1)return ans;
    p[last]=line_intersect(q[last],q[first]);
    for(int i=first;i<=last;i++)ans.push_back(p[i]);
    return ans;
}
int n;
int main()
{
    freopen(NAME".in","r",stdin);
    freopen(NAME".out","w",stdout);
    scanf("%d",&n);
    vector<Point>p,v,normal;
    for(int i=1;i<=n;i++)
    {   
        double x,y;
        scanf("%lf%lf",&x,&y);
        p.push_back(Point(x,y));

    }
    for(int i=0;i<=n-1;i++)
    {
        v.push_back(p[(i+1)%n]-p[i]);
        normal.push_back(Normal(v[i]));
    }
    double lf=0,rg=20000;
    while(rg-lf>eps)
    {
        vector<Line>h;
        double mid=lf+(rg-lf)/2.0;
        for(int i=0;i<=n-1;i++)
        h.push_back(Line(p[i]+normal[i]*mid,v[i]));
        vector<Point>poly=half_plane_intersect(h);
        if(poly.empty())rg=mid;
        else lf=mid;
    }
    printf("%.2lf",lf);
    return 0;
}
theljt
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
凸包问题的五种解法
baningtao1470的博客
05-29 3011
前言： 首先，什么是凸包？ 假设平面上有p0~p12共13个点，过某些点作一个多边形，使这个多边形能把所有点都“包”起来。当这个多边形是凸多边形的时候，我们就叫它“凸包”。如下图： 然后，什么是凸包问题？ 我们把这些点放在二维坐标系里面，那么每个点都能用 (x,y) 来表示。 现给出点的数目13，和各个点的坐标。求构成凸包的点？ 解一：穷举法（...
线段交点计算与凸包的简单实现
月无痕的博客
02-21 634
本文主要介绍了平面凸包计算。采用递增策略构建算法，并进行简单实现。同时介绍了一种求解线段交点的方法。
叉积、线段相交判断、凸包
albon arith
09-04 9321
一、叉积 叉积的计算是线段方法的核心。对于向来p1和p2，叉积是由点(0,0)、p1、p2和p1+p2构成的平行四边形的有向面积。另一种与之等价但更有效的的叉积定义方式是将其看做矩阵行列式： p1×p2 = x1y2 - x2y1 = - p2×p1 若p1×p2为正，则相对于原点(0,0)来说，p1位于p2顺时针方向；若p1×p2为负，p1位于p2逆时针方向；若为0则方向相同，
POJ - 1584 A Round Peg in a Ground Hole(综合几何)
Falcon的博客
01-28 220
题目链接：点击查看 题目大意：给出n个点，以及一个圆心和半径，首先判断这n个点能否构成凸包，若能，继续判断圆是否在凸包内 题目分析：这个题目确实非常综合，考察了判断凸包问题，判断圆是否在凸包内，只要保证圆心在凸包内，以及圆心到每条边的距离小于等于半径即可，说到底还是考察了模板是否靠谱，不得不佩服kuangbin大神的模板真的牛，一套用下来没有一点bug，直接贴代码吧，全是模板，懒得删减，所以有...
凸包
一个弱鸡的博客
08-04 448
一篇很不错的博客https://blog.csdn.net/bone_ace/article/details/46239187 贴个大佬博客找的板子 #include &lt;algorithm&gt; #include &lt;cstdio&gt; #include &lt;cmath&gt; const int N=1e5+5; const double eps=1e-6; int n,...
Python求凸包及多边形面积教程
09-17
在Python编程中，处理几何问题时，有时我们需要计算多边形的凸包和面积。凸包是一组点中所有点连接形成的最小多边形，该多边形包含了所有原始点。本教程将介绍如何使用Python实现这一功能。 首先，凸包的计算通常...
JTS.rar_JTS将面转为线_jts_jts 求凸多边形_jts 计算中心点_jts计算点中心
最新发布
09-21
接着，"jts_求凸多边形"指的是计算几何对象的凸包。JTS提供了一个`ConvexHull`类，可以通过传入一个几何对象来计算其凸包。凸包是包含几何对象所有点的最小凸多边形，这对于理解和简化复杂形状或者快速检测碰撞很有...
ff.rar_grahamscan_n个点求其凸包_凸包测试数据_求其凸包。
09-19
grahamScan扫描法求凸包 输入有若干组测试数据。每一组测试数据的第一行上有整数n，表示该组测试数据有n个点组成的。接下来有n行，其每一行上有二个正整数，之间用一个或几个空格隔开。当输入行上只有一个数0时，...
tubao.zip_visual c_凸多边形matlab_多边形凸包
09-20
在计算机图形学和几何处理领域，凸包（Convex Hull）是一个重要的概念，它是指一个几何对象所有点的最小凸集合。在这个“tubao.zip”压缩包中，包含的是关于使用C++编程语言和MATLAB软件实现多边形凸包算法的相关...
凸包 —— 五种解法
lxt_Lucia的博客
10-17 7055
欢迎访问https://blog.csdn.net/lxt_Lucia～～ 宇宙第一小仙女\(^o^)/～～萌量爆表求带飞=≡Σ((( つ^o^)つ~ dalao们点个关注呗～～   关于凸包，之前一直漏掉了这一个，上周末组队赛，听说有一题就是凸包的模板，可惜不会，听说不算很难，正好这周空余的课很少，就顺便搞了8...    目录 前言 解一：穷举法（蛮力法） 解二：分治法 解三...
判断线段相交，凸包
(ง •̀_•́)ง
03-11 824
判断线段相交  #include #include using namespace std; typedef struct dotnode{ double x; double y; }dot; double Direction(dot pi,dot pj,dot pk){ return (pk.x-pi.x)*(pj.y-pi.y)-(pj.x-pi.x)*(pk.y-p
判断点是否在凸包内
weixin_43829465的博客
08-20 4584
周测5 描述 二维平面上，给定n个点{ai}和m个点{bi}，且保证这n+m个点中，任意两个点的x坐标或y坐标均不相同。 对于每个bi，判断是否存在由3个ai,aj,ak(1≤i,j,k≤n,i≠j≠k)点组成的三角形包含bi（在三角形边上也算包含；允许三点共线的三角形，此时只有bi在三点中任意两点的线段上才算包含）。 输入 第一行为一个整数n。接下来n行，其中第i行有两个整数，表示ai的横纵坐标...
凸包问题
zy010101博客
08-19 1804
定义1：平面上的点集，如果以该集合中的任意两点P和Q为端点构成的线段属于该集合，就称该集合是凸的。 定义2：一个点集S的凸包是包含S的最小凸集合。 定理：任意包含n &gt; 2个点的集合S的凸包是以S中的某些点为顶点的凸多边形。（如果所有点是共线的，多边形退化为线段） 因此，直观看来，任意的凸多边形都是凸集合。 凸包问题是为一个包含n个点的集合构造一个凸包。 根据上面的定理设计了一个基...
NOIP2009解题报告(C/C++)（潜伏者）（Hankson的趣味题）（最优贸易）（靶形数独）
中国计算机学会成都七中（高新校区）委员会委员活动室
03-08 2299
NOIP2009解题报告
NOIP2010提高组复赛 解题报告（C/C++）（机械翻译）（乌龟棋）（关押罪犯）（引水入城）
中国计算机学会成都七中（高新校区）委员会委员活动室
02-18 1185
NOIP2010解题报告
17.8.27 校内赛 解题报告【卢卡斯定理】【dfs+hash】【线段树】
中国计算机学会成都七中（高新校区）委员会委员活动室
08-27 581
17.8.27 校内赛
17.4.4 数据结构 校内赛（setmod）（area）（intkth）
中国计算机学会成都七中（高新校区）委员会委员活动室
04-12 524
数据结构 校内赛
17.4.5 图论 校内赛
中国计算机学会成都七中（高新校区）委员会委员活动室
04-21 423
遥祝jyb同志万寿无疆！

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值