链接:https://ac.nowcoder.com/acm/contest/984/J
来源:牛客网
时间限制:C/C++ 1秒,其他语言2秒
空间限制:C/C++ 32768K,其他语言65536K
64bit IO Format: %lld
题目描述
为了防止口渴的食蚁兽进入他的农场,Farmer John决定在他的农场周围挖一条护城河。农场里一共有N(8<=N<=5,000)股泉水,并且,护城河总是笔直地连接在河道上的相邻的两股泉水。护城河必须能保护所有的泉水,也就是说,能包围所有的泉水。泉水一定在护城河的内部,或者恰好在河道上。当然,护城河构成一个封闭的环。
挖护城河是一项昂贵的工程,于是,节约的FJ希望护城河的总长度尽量小。请你写个程序计算一下,在满足需求的条件下,护城河的总长最小是多少。
所有泉水的坐标都在范围为(1..10,000,000,1..10,000,000)的整点上,一股泉水对应着一个唯一确定的坐标。并且,任意三股泉水都不在一条直线上。
以下是一幅包含20股泉水的地图,泉水用"*"表示:
...*-----------------*......
../..........*........\.....
./.....................\....
*......................*\...
|..........*........*....\..
|*........................\.
|..........................*
*..........................|
.\*........................|
..\.....................*..|
...\........*............*.|
....\..................*...*
.....\..*..........*....../.
......\................../..
.......*----------------*...
(请复制到记事本中用等宽字体查看)
图中的直线,为护城河的最优挖掘方案,即能围住所有泉水的最短路线。
路线从左上角起,经过泉水的坐标依次是:(18,0),(6,-6),(0,-5),(-3,-3),(-17,0),(-7,7),(0,4),(3,3)。绕行一周的路径总长为70.8700576850888(...)。答案只需要保留两位小数,于是输出是70.87。
输入描述:
第1行: 一个整数,N 第2..N+1行: 每行包含2个用空格隔开的整数,x[i]和y[i],即第i股泉水的位置坐标
输出描述:
第1行: 输出一个数字,表示满足条件的护城河的最短长度。保留两位小数
示例1
输入
复制
20 2 10 3 7 22 15 12 11 20 3 28 9 1 12 9 3 14 14 25 6 8 1 25 1 28 4 24 12 4 15 13 5 26 5 21 11 24 4 1 8
输出
复制
70.87
凸包模板题
#include <bits/stdc++.h>
#define fuck(x) cout<<'['<<#x<<' '<<(x)<<"]\n"
using namespace std;
typedef long long ll;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7;
const int MX = 1e6 + 5;
struct lp {
ll x, y;
bool operator<(const lp&b) const {
if(x == b.x) return y < b.y;
return x < b.x;
}
} cw[MX],R[MX],op[5],ok;
ll cross(lp a, lp b, lp c) {
return ((b.x - a.x) * (c.y - a.y) - (c.x - a.x) * (b.y - a.y));
}
int convex(int n) {
int tn = 0, k;
sort(cw, cw + n);
for(int i = 0; i < n; i++) {
while(tn > 1 && cross(R[tn - 1], cw[i], R[tn - 2]) <= 0) tn--;
R[tn++] = cw[i];
}
k = tn;
for(int i = n - 2; i >= 0; i--) {
while(tn > k && cross(R[tn - 1], cw[i], R[tn - 2]) <= 0) tn--;
R[tn++] = cw[i];
}
if(n > 1) tn--;
return tn;
}
double dis(lp a,lp b){
return sqrt((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("%lld%lld",&cw[i].x,&cw[i].y);
}
n = convex(n);
double ans=0;
R[n+1]=R[0];
for(int i=0;i<=n;i++){
ans+=dis(R[i+1],R[i]);
}
printf("%.2lf\n",ans);
return 0;
}