【计算几何】 凸包训练题解 （待更新）

最新推荐文章于 2021-05-12 20:50:00 发布
最新推荐文章于 2021-05-12 20:50:00 发布
阅读量889 收藏 1
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Qiuker_jl/article/details/110094498
版权

涉及知识点：
1.向量叉积及求凸包

文章目录

A - Beauty Contest

题目链接：POJ2187

题目大意：给定平面中N个点，求最远点对

分析

朴素思路：O( N 2 N^2 N2)枚举

优化：
①最远点对必定在凸包上
②暴力枚举凸包上的点，找最远距离
因为凸包上的点个数小于等于N，所以减少了时间。

上面的思路已经可以过题了（因为数据弱），但是细想一下，如果它给出的数据极端一点，所有点都在凸包上，那仍然会T；也就是说这种做法并不能百分百缩小枚举的范围

更优解法：
最远点对必定是对踵点，而对踵点对数不超过 3 2 N \frac{3}{2}N 23​N（网上看到的，并没有证明过）
所以如果可以枚举所有对踵点，就可以 O ( N ) O(N) O(N)解出
对踵点可以利用旋转卡壳O（N）求得

旋转卡壳的思想可以参考这篇某个巨佬写的文章：旋转卡壳，非常通俗易懂。

但是貌似有一点点小瑕疵，我做了一些修改。

凸包+暴力代码

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;

const double eps=1e-8;

typedef struct{
 int x,y;
}Point;

typedef Point Vector;

Point P[50005];
Point Ch[50005];

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

int Cross(Vector a,Vector b){
  return a.x*b.y-a.y*b.x;
}

bool cmp(Point a,Point b){
   if(a.x!=b.x) return a.x<b.x;
   return a.y<b.y;
}

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

int main()
{
    int N;
    scanf("%d",&N);
    for(int i=0;i<N;i++)
    scanf("%d%d",&P[i].x,&P[i].y);
    sort(P,P+N,cmp);

    int cnt=0;
    for(int i=0;i<N;i++){
      while(cnt>1&&Cross(Ch[cnt-2]-Ch[cnt-1],Ch[cnt-2]-P[i])<=0)
          cnt--;
      Ch[cnt++]=P[i];
    }
    int vis=cnt;
    for(int i=N-2;i>=0;i--){
       while(cnt>vis&&Cross(Ch[cnt-2]-Ch[cnt-1],Ch[cnt-2]-P[i])<=0)
          cnt--;
       Ch[cnt++]=P[i];
    }
    int ans=0;
    for(int i=0;i<cnt;i++)
        for(int j=i+1;j<cnt;j++)
            ans=max(ans,dis(Ch[i],Ch[j]));
    printf("%d\n",ans);
    return 0;
}

旋转卡壳代码

仔细观察Andrew求凸包的算法就会发现，初始点是会被重复计算两次的，所以凸包实际上点的个数应该是cnt-1（多算了一个）

给出一组数据

3
0 2
2 0
0 0

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

const double eps=1e-8;
const double PI=acos(-1);

typedef struct{
 int x,y;
}Point;

typedef Point Vector;

Point P[50005];
Point Ch[50005];

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

int Cross(Vector a,Vector b){
  return a.x*b.y-a.y*b.x;
}

bool cmp(Point a,Point b){
   if(a.x!=b.x) return a.x<b.x;
   return a.y<b.y;
}

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

int main()
{
    int N;
    scanf("%d",&N);
    for(int i=0;i<N;i++)
    scanf("%d%d",&P[i].x,&P[i].y);
    sort(P,P+N,cmp);

    int cnt=0;
    for(int i=0;i<N;i++){
      while(cnt>1&&Cross(Ch[cnt-2]-Ch[cnt-1],Ch[cnt-2]-P[i])<=0)
          cnt--;
      Ch[cnt++]=P[i];
    }
    int vis=cnt;
    for(int i=N-2;i>=0;i--){
       while(cnt>vis&&Cross(Ch[cnt-2]-Ch[cnt-1],Ch[cnt-2]-P[i])<=0)
          cnt--;
       Ch[cnt++]=P[i];
    }
    cnt--;
    //旋转卡壳：
    int ans=0;
    int j=2;
    for(int i=0;i<cnt;i++){
       while(abs(Cross(Ch[i]-Ch[i+1],Ch[i+1]-Ch[j]))<abs(Cross(Ch[i]-Ch[i+1],Ch[i+1]-Ch[j+1])))
         j=(j+1)%cnt;
       ans=max(ans,dis(Ch[i],Ch[j]));
    }
    printf("%d\n",ans);
    return 0;
}

B - Wall

题目链接：POJ1113

题目大意:

给出一个n边形的n个顶点，现要修建一堵墙，将这个n边形围起来，要求墙上任意一点距离n边形至少距离为L，求墙的最小周长

分析

第一想法很简单，就是让这个n边形外扩L
①题目没有说这个n边形是凸多边形，画图即可发现，如果是凹多边形，转化为凸包可以满足条件，且墙周长更短

②在凸多边形的角的地方可以磨圆，即以凸多边形的顶点为圆心，L为半径画圆即可，所得圆与外扩L的n边形的边相切，得到一个扇形。不难证明扇形的圆心角的和必然为360°：

①扇形的圆心角和其所对角和为180°，其所对角为凸多边形的内角
②总共n个角，所以圆心角之和加凸多边形内角之和为180n
③n凸边形内角和180（n-2）
④相减得圆心角之和为360

因此答案为凸包的周长加上以L为半径的圆的周长

代码

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

const double eps=1e-8;
const double PI=acos(-1);

typedef struct{
 int x,y;
}Point;

typedef Point Vector;

Point P[50005];
Point Ch[50005];

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

int Cross(Vector a,Vector b){
  return a.x*b.y-a.y*b.x;
}

bool cmp(Point a,Point b){
   if(a.x!=b.x) return a.x<b.x;
   return a.y<b.y;
}

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


int main()
{
    int N,L;
    scanf("%d%d",&N,&L);
    for(int i=0;i<N;i++)
    scanf("%d%d",&P[i].x,&P[i].y);
    sort(P,P+N,cmp);

    int cnt=0;
    for(int i=0;i<N;i++){
      while(cnt>1&&Cross(Ch[cnt-2]-Ch[cnt-1],Ch[cnt-2]-P[i])<=0)
          cnt--;
      Ch[cnt++]=P[i];
    }
    int vis=cnt;
    for(int i=N-2;i>=0;i--){
       while(cnt>vis&&Cross(Ch[cnt-2]-Ch[cnt-1],Ch[cnt-2]-P[i])<=0)
          cnt--;
       Ch[cnt++]=P[i];
    }

     double ans=0;
     for(int i=0;i<cnt;i++)
         ans+=dis(Ch[i],Ch[(i+cnt-1)%cnt]);
     ans+=2.0*PI*L;
     printf("%.0f\n",ans);
    return 0;
}

E - Triangle

题目大意：
给定平面坐标中N个点，从这N个点中任取三个，求组成的最大的三角形的面积是多少

分析

易得面积最大的三角形肯定由凸包上的三个点组成

朴素做法：
暴力枚举三个点
O ( N 3 ) O(N^3) O(N3)

优化：
枚举第一个顶点，旋转卡壳枚举第二个和第三个顶点，因为在凸包中两点固定时，第三个点按逆时针旋转时，所得三角形面积为单峰函数，所以复杂度大致降到了 O ( N 2 ) O(N^2) O(N2)

代码

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

const double eps=1e-8;
const double PI=acos(-1);

typedef struct{
 int x,y;
}Point;

typedef Point Vector;

Point P[50005];
Point Ch[50005];

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

int Cross(Vector a,Vector b){
  return a.x*b.y-a.y*b.x;
}

bool cmp(Point a,Point b){
   if(a.x!=b.x) return a.x<b.x;
   return a.y<b.y;
}

int main()
{
    int N;
    while(~scanf("%d",&N)&&N!=-1){
       for(int i=0;i<N;i++)
           scanf("%d%d",&P[i].x,&P[i].y);
       sort(P,P+N,cmp);

       int cnt=0;
       for(int i=0;i<N;i++){
          while(cnt>1&&Cross(Ch[cnt-2]-Ch[cnt-1],Ch[cnt-2]-P[i])<=0)
            cnt--;
          Ch[cnt++]=P[i];
       }
       int vis=cnt;
       for(int i=N-2;i>=0;i--){
          while(cnt>vis&&Cross(Ch[cnt-2]-Ch[cnt-1],Ch[cnt-2]-P[i])<=0)
            cnt--;
          Ch[cnt++]=P[i];
       }
       cnt--;
       int ans=0;
       int j=1;
       int k=2;
       for(int i=0;i<cnt;i++){
         while(abs(Cross(Ch[i]-Ch[j],Ch[j]-Ch[k]))<abs(Cross(Ch[i]-Ch[j],Ch[j]-Ch[k+1])))
           k=(k+1)%cnt;
         while(abs(Cross(Ch[i]-Ch[j],Ch[j]-Ch[k]))<abs(Cross(Ch[i]-Ch[j+1],Ch[j+1]-Ch[k])))
           j=(j+1)%cnt;
         ans=max(ans,abs(Cross(Ch[i]-Ch[j],Ch[j]-Ch[k])));
       }
       printf("%.2f\n",(double)ans/2.0);
    }
    return 0;
}

K - Romeo Meets Juliet

题目链接：POJ-2189

这题暂时没想到和凸包有啥关系， 1 0 3 10^3 103的数据量也完全可以 N 2 N^2 N2暴力

代码

#include <cstdio>
#include <iostream>
#include <queue>
#include <cstring>
using namespace std;

int cow[1005];

int  main()
{
    int N,P,C;
    scanf("%d%d%d",&N,&P,&C);
    for(int i=1;i<=N;i++){
        int temp;
        scanf("%d",&temp);
        cow[temp]++;
    }

    int ans=0;
    for(int i=1;i<P;i++){
        int cnt=0;
        for(int j=i;j<P;j++){
              cnt+=cow[j];
              if(cnt>C) break;
              ans=max(ans,j-i+1);
        }
    }
    printf("%d\n",ans);
    return 0;
}

L - Scrambled Polygon

题目链接：POJ-2007

题目大意:

给定一个n边形，n个顶点的坐标，保证其中一个顶点在原点，求这些顶点逆时针排序后的结果

分析

利用叉积判断顺逆时针关系

代码

#include<cstdio>
#include <algorithm>
#include <cmath>
using namespace std;

typedef struct{
  int x;
  int y;
}Point;

typedef Point Vector;

int Cross(Vector a,Vector b){
    return a.x*b.y-a.y*b.x;
}

bool cmp(Vector a,Vector b){
    return a.x*b.y-a.y*b.x>0;
}

Vector V[1000];

int main() {
  int x,y;
  int cnt=0;
  while(~scanf("%d%d",&x,&y))
     V[cnt++]=Point{x-0,y-0};
  sort(V+1,V+cnt,cmp);
  for(int i=0;i<cnt;i++)
     printf("(%d,%d)\n",V[i].x,V[i].y);
	return 0;
}
Nefu_qky
关注
计算几何 | 凸包入门及例题
Encore47的博客
11-11 528
学习笔记 Andrew算法 算法做两次扫描，先从最左边的点沿下凸包扫描到最右边，再从最右边的点沿上凸包扫描到最左边，上下合在一起是完整的凸包。 例题 模板题：求凸包周长 以下convex_hull函数为核心求解凸包的函数，凸包顶点放在ch数组中，返回值为顶点数。 //计算几何 int sgn(db x) { //判断浮点数是否为0，为0时返回0 if (fabs(x) < eps) return 0; return x < 0 ? -1 : 1; } struct poi
计算几何之凸包问题
m0_49959202的博客
07-17 157
文章目录计算几何之凸包问题POJ 1113 WallPOJ 2007 Scrambled PolygonPOJ 1873 The Fortified ForestPOJ 1228 Grandpa's EstatePOJ 3348 Cows 计算几何之凸包问题 POJ 1113 Wall 大意： 给出一个n边形的城堡，国王要求修建一个城墙，城墙到城堡的距离最少是L，问最少的城堡边长是多少，答案要求四舍五入到整数 [外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-sZDZ
学习凸包(三）:凸包练习 POJ 1113
04-14
NULL 博文链接：https://128kj.iteye.com/blog/1748635
【总结】Graham求凸包 及 习题
秋叶红于二月花
10-19 219
计算几何模板之一： 凸包 大致意思就是先极角排序，再用单调栈维护凸包 有参考计算几何之凸包模板 实现中，令s[cnt+1]=s[1]，方便循环。 struct Pt{ double x,y; }p[M], s[M]; //调用graham(n)，将p[1..n]的点求凸包放在s[1..cnt]中，返回cnt namespace Convex //求凸包 { double cross(Pt ...
几道凸包题
@火柴
08-07 718
算法
凸包（讲解+模板+例题）
September_的博客
07-31 931
凸包讲解 例题POJ 2187 模板 #include<cstdio> #include<iostream> #include<algorithm> #include<cstring> #include<cmath> using namespace std; struct node { int x,y; }; node...
算法参考资料国际大学生程序设计竞赛例题解1数论,计算几何,搜索算法专集
最新发布
04-15
此资料集题为《算法参考资料国际大学生程序设计竞赛例题解1数论、计算几何、搜索算法专集》。从标题来看，这份资料是针对国际大学生程序设计竞赛（ICPC）、区域性的算法竞赛（如ACM-ICPC）或是其他算法比赛中的问题...
题解 计算几何刷题记录
Love_Jacques的博客
05-12 328
题单源自洛谷：16408 计算几何 - 从入门到跳楼： 入门： P1652 圆： 题意：给n个不相交的圆，求从Point1到Point2一条曲线最少穿越多少个圆。 思路：n个圆都不相交，且连接两点的曲线可以任意弯曲延展，则除了只包括point1或者只包括point2的圆不能避开，其他都可以避开（同时包括point1和point2不会被影响）。 代码： #include<bits/stdc++.h> #define eps 1e-6 using namespace std; struct p
HDU 4273(计算几何+凸包重心)
wust_cyl的博客
09-21 529
问题描述: I work at NASA outer space rescue team which needs much courage and patient. In daily life, I always receive a lot of mission, and I must complete it right now.  Today, team leader announced m
计算几何 ( 求凸包，计算三角形面积 )——最大三角形 ( HDU 2202 )
FeBr2的博客
07-25 2693
题目链接： http://acm.hdu.edu.cn/showproblem.php?pid=2202 分析： 给出平面上的N个点，在其中找出三个点，使得其构成的三角形的面积最大。 题解： 先用找出凸包上的所有点，然后一次枚举遍历即可。 1.求凸包：int cmp(point a, point b) //水平排序 { if(a.x==b.x)return a.y<b.y;
凸包详细代码+解说+题目
08-22
zoj上凸包题目，可用作模板，相信您用了之后一定不会后悔的
chebycenter(A,b,r0):计算多边形的切比雪夫中心和切比雪夫半径。-matlab开发
05-30
此函数返回由线性不等式 Ax<=b 定义的多胞体的切比雪夫中心和切比雪夫半径。 多胞体的切比雪夫中心是由多胞体包围的最大超球面的中心。 类似地，切比雪夫半径是由多面体包围的最大超球面的半径。
[hdu5251]矩形面积 凸包+旋转卡壳 计算几何
LvYanchang的博客
05-23 257
Problem Description小度熊有一个桌面，小度熊剪了很多矩形放在桌面上，小度熊想知道能把这些矩形包围起来的面积最小的矩形的面积是多少。 Input第一行一个正整数 T，代表测试数据组数(1≤T≤20)，接下来 T 组测试数据。每组测试数据占若干行，第一行一个正整数 N(1≤N&lt;≤1000)，代表矩形的数量。接下来 N 行，每行 8 个整数x1,y1,x2,y2,x3,y3,x4...
hdu 1348 Wall（凸包模板题 Graham）
08-15 1821
hdu 1348 Wall（凸包模板题）
凸包例题-UVa-10652
qq_40905284的博客
02-20 296
题意：有n块矩形木板，你的任务是用一个面积尽量小的凸多边形把它们包起来，并计算出木板占整个包装面积的百分比。 第一行输入的是数据的组数，（T<50)。每组数据第一行为木板个...
POJ 1113 Wall 凸包(模板题)
Masamiiiiii
04-18 1314
用不严谨的话来讲，给定二维平面上的点集，凸包就是将最外层的点连接起来构成的凸多边型，它能包含点集中所有的点。 题意:给出n个点 n=L? 若不考虑L 则答案为凸包的周长,考虑条件L,由如下图形 答案为:凸包周长+以L为半径的圆的周长   求凸包用的是 Graham's Scan法 可以根据归纳法来证明,栈中保存的是前i-1个点集的凸包的极点 因为按极角排序后,若(p[i],栈顶
小蒟蒻的题解（代码）系列-----凸包入门题目集合
子衿君的博客
03-04 688
学习过斜率优化的童鞋们应该对凸包都不陌生，这里就不做详细的介绍了 手动上代码 P2742 【模板】二维凸包 #include &lt;algorithm&gt; #include &lt;iostream&gt; #include &lt;cstdio&gt; #include &lt;cmath&gt; using namespace std; const int N=100010; const...
大白 计算几何专题 凸包、半平面交、平面区域 部分例题练习题总结
lchi1997的博客
04-26 584
LA 2797 Monster Trap 1. 姿势，学习学习， 做法就是将每条线段都延长一些些， 然后判断点到点是否可以直接到达再bfs。 2. 最重要的就是避免穿墙， 除了书上讲的情况外，还要排除 0 0 4 0 1 0 2 0 这一组数据 或者 0 0 2 0 2 0 4 0 这一组数据 当然，如果判断点在线段上的条件不同对这组数据的鲁棒性就不同。 代码：
习题5-1(凸壳、拓扑排序)
Markland_l的博客
10-31 476
习题5-1(凸壳、拓扑排序) 凸包 卷包裹 求凸壳 解法1 求上下凸壳，类似卷包裹，但是不是极角排序(任取一个点，做极角排序) 按照先x轴后y轴的方式升序排序，这是为了找到一个一定在凸包上的点 然后求上凸壳和下凸壳，可以想象在y轴负无穷处有一个点，然后求凸包得到上凸壳，y轴正无穷有一个点，然后求凸包得到下凸壳 单调栈 代码解析 convex函数，对a数组求凸包，求的结果存储在b数组中 可以想象在y轴负无穷处有一个点，然后求凸包得到上凸壳，上凸壳加上y轴无穷远处的点，可以把内部
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值