博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
HDU 1392 Surround the Trees(凸包)题解
阅读量:5837 次
发布时间:2019-06-18

本文共 2745 字,大约阅读时间需要 9 分钟。

题意:给一堆二维的点,问你最少用多少距离能把这些点都围起来

思路:

凸包:

我们先找到所有点中最左下角的点p1,这个点绝对在凸包上。接下来对剩余点按照相对p1的角度升序排序,角度一样按距离升序排序。因为凸包有一个特点,从最左下逆时针走,所有线都在当前这条向量的左边,根据这个特点我们进行判断。我们从栈顶拿出两个点s[top-1],s[top],所以如果s[top-1] -> p[i] 在 s[top-1] -> s[top] 右边,那么s[top]就不是凸包上一点,就这样一直判断下去。判断左右可以用叉乘。

参考:

模板(考虑n <= 2):

struct node{    double x, y;}p[maxn], s[maxn];int n, top;double dis(node a, node b){    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));}bool cmp(node a, node b){    double A = atan2((a.y - p[1].y), (a.x - p[1].x));    double B = atan2((b.y - p[1].y), (b.x - p[1].x));    if(A != B) return A < B;    else{        return dis(a, p[1]) < dis(b, p[1]);    }}double cross(node a, node b, node c){   //(a->b)X(a->c)    return (b.x - a.x) * (c.y - a.y) - (c.x - a.x) * (b.y - a.y);}void solve(){    int pos = 1;    for(int i = 2; i <= n; i++){        if(p[i].y < p[pos].y || (p[i].y == p[pos].y && p[i].x < p[pos].x)){            pos = i;        }    }    swap(p[1], p[pos]);    sort(p + 2, p + n + 1, cmp);    s[0] = p[1], s[1] = p[2];    top = 1;    for(int i = 3; i <= n; i++){        while(top >= 1 && cross(s[top - 1], p[i], s[top]) >= 0){ //向右转出栈            top--;        }        s[++top] = p[i];    }}

 

代码:

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
typedef long long ll;using namespace std;const int maxn = 150 + 10;const int MOD = 1e9 + 7;const int INF = 0x3f3f3f3f;struct node{ double x, y;}p[maxn], s[maxn];int n, top;double dis(node a, node b){ return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));}bool cmp(node a, node b){ double A = atan2((a.y - p[1].y), (a.x - p[1].x)); double B = atan2((b.y - p[1].y), (b.x - p[1].x)); if(A != B) return A < B; else{ return dis(a, p[1]) < dis(b, p[1]); }}double cross(node a, node b, node c){ //(a->b)X(a->c) return (b.x - a.x) * (c.y - a.y) - (c.x - a.x) * (b.y - a.y);}void solve(){ int pos = 1; for(int i = 2; i <= n; i++){ if(p[i].y < p[pos].y || (p[i].y == p[pos].y && p[i].x < p[pos].x)){ pos = i; } } swap(p[1], p[pos]); sort(p + 2, p + n + 1, cmp); s[0] = p[1], s[1] = p[2]; top = 1; for(int i = 3; i <= n; i++){ while(top >= 1 && cross(s[top - 1], p[i], s[top]) >= 0){ //向左转出栈 top--; } s[++top] = p[i]; }}int main(){ while(~scanf("%d", &n) && n){ for(int i = 1; i <= n; i++){ scanf("%lf%lf", &p[i].x, &p[i].y); } if(n == 1){ printf("0.00\n"); continue; } if(n == 2){ printf("%.2lf\n", dis(p[1], p[2])); continue; } solve(); double ans = 0; for(int i = 0; i < top; i++){ ans += dis(s[i], s[i + 1]); } ans += dis(s[top], s[0]); printf("%.2lf\n", ans); } return 0;}

 

转载于:https://www.cnblogs.com/KirinSB/p/10414378.html

你可能感兴趣的文章
【C#】protected 变量类型
查看>>
Ubuntu解压
查看>>
爬虫_房多多(设置随机数反爬)
查看>>
藏地密码
查看>>
爬虫去重(只是讲了去重的策略,没有具体讲实现过程,反正就是云里雾里)...
查看>>
react中将px转化为rem或者vw
查看>>
8816
查看>>
avcodec_open2()分析
查看>>
何如获取单选框中某一个选中的值
查看>>
paip.输入法编程----删除双字词简拼
查看>>
tcp状态
查看>>
QQ悬浮返回顶部
查看>>
MySQL建表语句的一些特殊字段
查看>>
《Unix环境高级编程》读书笔记 第8章-进程控制
查看>>
腾讯前端二面题目详解
查看>>
mascara-1
查看>>
Jquery Form表单取值
查看>>
Python version 2.7 required, which was not found in the registry
查看>>
Android API level 与version对应关系
查看>>
Team Name
查看>>