二维凸包
人生第一道计算几何题? 留作纪念
凸包是什么?下图就是一个土包
我们怎么求出凸包呢
我们先把凸包分成上下两个部分
如图,红色部分称为上半部分,黑色部分成为下半部分
然后按照x为第一关键字,y为第二关键字排序
于是第一个就会找到最左边的一个最下面的点,第二个就会找到最右边的一个最上面的点
我们以这两个点为分界线,分出上下两个凸包分别求出周长
我们维护一个栈,栈中的元素即为上(下)凸包所有元素
对于下半部分,我们显然要让斜率尽可能地小
我们先把前两个点(分别为AB)入栈,拓展到第三个点(C),若则把C进栈,否则将B弹出并将C进栈,不断弹出直到
而对于上半部分,我们是要让斜率尽可能大
同理,我们先把后两个点(分别为FE)入栈,拓展到第三个点(D),若则把D进栈,否则将E弹出并将D进栈,不断弹出直到
于是我们就可以得到上下两个凸包的所有点,再用距离公式算就可以了,时间复杂度
代码如下:
#include<bits/stdc++.h>
using namespace std;
#define il inline
#define re register
#define D double
#define inf 123456789
il int read()
{
re int x = 0, f = 1; re char c = getchar();
while(c < '0' || c > '9') { if(c == '-') f = -1; c = getchar();}
while(c >= '0' && c <= '9') x = x * 10 + c - 48, c = getchar();
return x * f;
}
#define rep(i, s, t) for(re int i = s; i <= t; ++ i)
#define drep(i, s, t) for(re int i = t; i >= s; -- i)
#define Next(i, u) for(re int i = head[u]; i; i = e[i].next)
#define mem(k, p) memset(k, p, sizeof(k))
#define maxn 10005
struct node{D x, y;}e[maxn];
int n, st[maxn], top;
D ans;
//il D ff(int x){return x * x;}
#define ff(x) (x) * (x)
il bool cmp(node a, node b){return a.x == b.x ? a.y < b.y : a.x < b.x;}
il D dis(int x, int y){return sqrt(ff(e[x].x - e[y].x) + ff(e[x].y - e[y].y));}
il D getk(int x, int y){return e[x].x == e[y].x ? inf : (e[x].y - e[y].y) / (e[x].x - e[y].x);}//注意这里一定要判断两个x相等的情况
int main()
{
n = read();
rep(i, 1, n) scanf("%lf%lf", &e[i].x, &e[i].y);
sort(e + 1, e + n + 1, cmp);
rep(i, 1, n)
{
st[++ top] = i;
while(top > 2 && getk(st[top], st[top - 2]) < getk(st[top - 1], st[top - 2]))
st[top - 1] = st[top], -- top;
}
rep(i, 1, top - 1) ans += dis(st[i], st[i + 1]);
top = 0;
drep(i, 1, n)
{
st[++ top] = i;
while(top > 2 && getk(st[top], st[top - 2]) < getk(st[top - 1], st[top - 2]))
st[top - 1] = st[top], -- top;
}
rep(i, 1, top - 1) ans += dis(st[i], st[i + 1]);
printf("%.2lf", ans);
return 0;
}
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· PostgreSQL 和 SQL Server 在统计信息维护中的关键差异
· C++代码改造为UTF-8编码问题的总结
· DeepSeek 解答了困扰我五年的技术问题
· 为什么说在企业级应用开发中,后端往往是效率杀手?
· 用 C# 插值字符串处理器写一个 sscanf
· [翻译] 为什么 Tracebit 用 C# 开发
· Deepseek官网太卡,教你白嫖阿里云的Deepseek-R1满血版
· DeepSeek崛起:程序员“饭碗”被抢,还是职业进化新起点?
· 2分钟学会 DeepSeek API,竟然比官方更好用!
· .NET 使用 DeepSeek R1 开发智能 AI 客户端