Poj 1113 Wall [凸包]

题目大意:

给定多边形城堡的n个顶点,绕城堡外面建一个围墙,围住所有点,并且墙与所有点的距离至少为L,求这个墙最小的长度

思路:

结果等于凸包周长 + 一个完整圆的周长。

当绕城堡外面走一圈,在拐角处所形成的弧度相加为一个长度为L的圆。

  1 /*
  2     poj 1113
  3     凸包 
  4     大意:
  5     给定多边形城堡的n个顶点,绕城堡外面建一个围墙,围住所有点,
  6     并且墙与所有点的距离至少为L,求这个墙最小的长度 
  7 */
  8 
  9 #include <iostream>
 10 #include <math.h>
 11 #include <cstring> 
 12 #include <algorithm>
 13 #include <stdlib.h>
 14 
 15 using namespace std;
 16 #define pi acos(-1.0)
 17 #define eps 1e-8
 18 #define zero(x) (((x)>0?(x):-(x))<eps)
 19 
 20 struct point
 21 {
 22     double x, y;
 23 }p[1005], convex[1005];
 24 
 25 //两点间距离 
 26 double dis(point a, point b)
 27 {
 28     return sqrt((a.x - b.x)*(a.x - b.x) + (a.y - b.y) * (a.y - b.y));
 29 }
 30 
 31 //计算cross product (P1-P0)x(P2-P0)
 32 //叉积 
 33 double xmult(point p1, point p2, point p0)
 34 {
 35     return (p1.x - p0.x)*(p2.y - p0.y) - (p2.x - p0.x)*(p1.y - p0.y);
 36 }
 37 //graham算法顺时针构造包含所有共线点的凸包,O(nlogn)
 38 point p1, p2;
 39 
 40 int graham_cp(const void* a, const void* b)
 41 {
 42     double ret = xmult(*((point*) a), *((point*) b), p1);
 43     double ret2 =  xmult(*((point*) a), *((point*) b), p2);
 44     return zero(ret) ? (ret2> 0 ? 1 : -1) : (ret > 0 ? 1 : -1);
 45 }
 46 
 47 void _graham(int n, point* p, int& s, point* ch)
 48 {
 49     int i, k = 0;
 50     for (p1 = p2 = p[0], i = 1; i<n; p2.x += p[i].x, p2.y += p[i].y, i++)
 51         if (p1.y - p[i].y>eps || (zero(p1.y - p[i].y) && p1.x > p[i].x))
 52             p1 = p[k = i];
 53     p2.x /= n, p2.y /= n;
 54     p[k] = p[0], p[0] = p1;
 55     qsort(p + 1, n - 1, sizeof(point), graham_cp);
 56     for (ch[0] = p[0], ch[1] = p[1], ch[2] = p[2], s = i = 3; i < n; ch[s++] = p[i++])
 57         for (; s>2 && xmult(ch[s - 2], p[i], ch[s - 1]) < -eps; s--);
 58 }
 59 
 60 int wipesame_cp(const void *a, const void *b)
 61 {
 62     if ((*(point *) a).y < (*(point *) b).y - eps) return -1;
 63     else if ((*(point *) a).y > (*(point *) b).y + eps) return 1;
 64     else if ((*(point *) a).x < (*(point *) b).x - eps) return -1;
 65     else if ((*(point *) a).x > (*(point *) b).x + eps) return 1;
 66     else return 0;
 67 }
 68 
 69 int _wipesame(point * p, int n)
 70 {
 71     int i, k;
 72     qsort(p, n, sizeof(point), wipesame_cp);
 73     for (k = i = 1; i < n; i++)
 74         if (wipesame_cp(p + i, p + i - 1) != 0) p[k++] = p[i];
 75     return k;
 76 }
 77 
 78 //构造凸包接口函数,传入原始点集大小n,点集p(p原有顺序被打乱!)
 79 //返回凸包大小,凸包的点在convex中
 80 //参数maxsize为1包含共线点,为0不包含共线点,缺省为1
 81 //参数dir为1顺时针构造,为0逆时针构造,缺省为1
 82 //在输入仅有若干共线点时算法不稳定,可能有此类情况请另行处理!
 83 int graham(int n, point* p, point* convex, int maxsize = 1, int dir = 1)
 84 {
 85     point* temp = new point[n];
 86     int s, i;
 87     n = _wipesame(p, n);
 88     _graham(n, p, s, temp);
 89     for (convex[0] = temp[0], n = 1, i = (dir ? 1 : (s - 1)); dir ? (i < s) : i; i += (dir ? 1 : -1))
 90         if (maxsize || !zero(xmult(temp[i - 1], temp[i], temp[(i + 1)%s])))
 91             convex[n++] = temp[i];
 92     delete []temp;
 93     return n;
 94 }
 95 
 96 int main()
 97 {
 98     int n;
 99     double d;
100     while (cin >> n >> d)
101     {
102         memset(p,0,sizeof(p));
103         memset(convex, 0, sizeof(convex));
104         double dist = 0.0;
105         for (int i = 0; i < n; i++)
106         {
107             cin >> p[i].x >> p[i].y;
108         }
109         //cout << n << endl;
110         int size = graham(n, p, convex,0);
111 //        cout << size << endl;
112         for (int i = 0; i < size - 1; i++)
113         {
114             dist += dis(convex[i], convex[i + 1]);
115         }
116         dist += dis(convex[size - 1], convex[0]);//别忘了起始到结束的长度
117         cout << (int) (dist + 2 * pi * d + 0.5) << endl;
118     }
119 }


 

 

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
SQLAlchemy 是一个 SQL 工具包和对象关系映射(ORM)库,用于 Python 编程语言。它提供了一个高级的 SQL 工具和对象关系映射工具,允许开发者以 Python 类和对象的形式操作数据库,而无需编写大量的 SQL 语句。SQLAlchemy 建立在 DBAPI 之上,支持多种数据库后端,如 SQLite, MySQL, PostgreSQL 等。 SQLAlchemy 的核心功能: 对象关系映射(ORM): SQLAlchemy 允许开发者使用 Python 类来表示数据库表,使用类的实例表示表中的行。 开发者可以定义类之间的关系(如一对多、多对多),SQLAlchemy 会自动处理这些关系在数据库中的映射。 通过 ORM,开发者可以像操作 Python 对象一样操作数据库,这大大简化了数据库操作的复杂性。 表达式语言: SQLAlchemy 提供了一个丰富的 SQL 表达式语言,允许开发者以 Python 表达式的方式编写复杂的 SQL 查询。 表达式语言提供了对 SQL 语句的灵活控制,同时保持了代码的可读性和可维护性。 数据库引擎和连接池: SQLAlchemy 支持多种数据库后端,并且为每种后端提供了对应的数据库引擎。 它还提供了连接池管理功能,以优化数据库连接的创建、使用和释放。 会话管理: SQLAlchemy 使用会话(Session)来管理对象的持久化状态。 会话提供了一个工作单元(unit of work)和身份映射(identity map)的概念,使得对象的状态管理和查询更加高效。 事件系统: SQLAlchemy 提供了一个事件系统,允许开发者在 ORM 的各个生命周期阶段插入自定义的钩子函数。 这使得开发者可以在对象加载、修改、删除等操作时执行额外的逻辑。
SQLAlchemy 是一个 SQL 工具包和对象关系映射(ORM)库,用于 Python 编程语言。它提供了一个高级的 SQL 工具和对象关系映射工具,允许开发者以 Python 类和对象的形式操作数据库,而无需编写大量的 SQL 语句。SQLAlchemy 建立在 DBAPI 之上,支持多种数据库后端,如 SQLite, MySQL, PostgreSQL 等。 SQLAlchemy 的核心功能: 对象关系映射(ORM): SQLAlchemy 允许开发者使用 Python 类来表示数据库表,使用类的实例表示表中的行。 开发者可以定义类之间的关系(如一对多、多对多),SQLAlchemy 会自动处理这些关系在数据库中的映射。 通过 ORM,开发者可以像操作 Python 对象一样操作数据库,这大大简化了数据库操作的复杂性。 表达式语言: SQLAlchemy 提供了一个丰富的 SQL 表达式语言,允许开发者以 Python 表达式的方式编写复杂的 SQL 查询。 表达式语言提供了对 SQL 语句的灵活控制,同时保持了代码的可读性和可维护性。 数据库引擎和连接池: SQLAlchemy 支持多种数据库后端,并且为每种后端提供了对应的数据库引擎。 它还提供了连接池管理功能,以优化数据库连接的创建、使用和释放。 会话管理: SQLAlchemy 使用会话(Session)来管理对象的持久化状态。 会话提供了一个工作单元(unit of work)和身份映射(identity map)的概念,使得对象的状态管理和查询更加高效。 事件系统: SQLAlchemy 提供了一个事件系统,允许开发者在 ORM 的各个生命周期阶段插入自定义的钩子函数。 这使得开发者可以在对象加载、修改、删除等操作时执行额外的逻辑。
GeoPandas是一个开源的Python库,旨在简化地理空间数据的处理和分析。它结合了Pandas和Shapely的能力,为Python用户提供了一个强大而灵活的工具来处理地理空间数据。以下是关于GeoPandas的详细介绍: 一、GeoPandas的基本概念 1. 定义 GeoPandas是建立在Pandas和Shapely之上的一个Python库,用于处理和分析地理空间数据。 它扩展了Pandas的DataFrame和Series数据结构,允许在其中存储和操作地理空间几何图形。 2. 核心数据结构 GeoDataFrame:GeoPandas的核心数据结构,是Pandas DataFrame的扩展。它包含一个或多个列,其中至少一列是几何列(geometry column),用于存储地理空间几何图形(如点、线、多边形等)。 GeoSeries:GeoPandas中的另一个重要数据结构,类似于Pandas的Series,但用于存储几何图形序列。 二、GeoPandas的功能特性 1. 读取和写入多种地理空间数据格式 GeoPandas支持读取和写入多种常见的地理空间数据格式,包括Shapefile、GeoJSON、PostGIS、KML等。这使得用户可以轻松地从各种数据源中加载地理空间数据,并将处理后的数据保存为所需的格式。 2. 地理空间几何图形的创建、编辑和分析 GeoPandas允许用户创建、编辑和分析地理空间几何图形,包括点、线、多边形等。它提供了丰富的空间操作函数,如缓冲区分析、交集、并集、差集等,使得用户可以方便地进行地理空间数据分析。 3. 数据可视化 GeoPandas内置了数据可视化功能,可以绘制地理空间数据的地图。用户可以使用matplotlib等库来进一步定制地图的样式和布局。 4. 空间连接和空间索引 GeoPandas支持空间连接操作,可以将两个GeoDataFrame按照空间关系(如相交、包含等)进行连接。此外,它还支持空间索引,可以提高地理空间数据查询的效率。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值