题意:给出n个点,坐标都是整数,用一个多边形使所有点严格位于多边形内部,求最小的周长。多边形的边只能是格子的边或者对角线。
这个题可以用凸包解,但我用的并不是凸包,因为题目还是限制条件比较多的,完全自己处理也能处理出来。
首先对于x相同的所有的点,对多边形周长有影响的只有y的最大值和y的最小值。然后考虑y的最大值,对于x不同的情况,如果两个比较大的y中间有一个小的y,中间的y就不会对结果有影响,所以把每一个x对应的y最大的点拿出来,根据y从大到小排列,然后以第一个点为基准,向两边延伸,向x增大和x减小方向是完全相同的。x减小时,如果x减小的比y减小的值大,那么将线按照题目要求的方式连接到当前的点,将基准点更新为当前点,如果x减小值比y减小的值小,那么就不能直接更新到当前点的位置了,只能将基准点的x, y都减小x的减小值。即连一条斜45度的线到当前的x,但是要记录当前点到更新的基准点的距离,将最后一个点的距离加在结果里。对于每一个x对应的y最小的点,将y乘-1加入到点集里就可以按照与另一侧完全相同的方法处理。
然后剩下两端的敞口就比较简单了,因为x都是相同的,所以就是两端的x对应的y的最大与最小的差,当然每一端还要加上两个45度斜线,这样才能严格位于内部。
说的比较多比较乱,看不懂就结合代码看吧
#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <algorithm>
#include <cmath>
#include <stdlib.h>
#include <map>
using namespace std;
struct point{
int x, y;
point(int x, int y):x(x), y(y){}
};
bool cmp(point a, point b) {
if(a.y != b.y) {
return a.y > b.y;
}
return a.x > b.x;
}
map<int, vector<int> > ss;
map<int, vector<int> >::iterator p;
vector<point> t;
void judge(int &u, int &v) {
if(t.size() < 2) return;
int nu = 0, nv = 0;
sort(t.begin(), t.end(), cmp);
int x, y;
int mnx = t[0].x, mny = t[0].y, mxx = t[0].x, mxy = t[0].y;
for(int i = 1; i < t.size(); i++) {
if(t[i].x < mnx) {
x = mnx - t[i].x;
y = mny - t[i].y;
if(x > y) {
u += x - y;
v += y;
mnx = t[i].x;
mny = t[i].y;
nu = 0;
}
else {
mnx -= x;
mny -= x;
v += x;
nu = y - x;
}
}
else if(t[i].x > mxx) {
x = t[i].x - mxx;
y = mxy - t[i].y;
if(x > y) {
u += x - y;
v += y;
mxx = t[i].x;
mxy = t[i].y;
nv = 0;
}
else {
mxx += x;
mxy -= x;
v += x;
nv = y - x;
}
}
}
u += nu + nv;
}
int main() {
int n, x, y;
while(~scanf("%d", &n)) {
ss.clear();
for(int i = 0; i < n; i++) {
scanf("%d%d", &x, &y);
if(ss.count(x) == 0) {
ss[x].push_back(y);
ss[x].push_back(y);
}
else {
if(y < ss[x][0])
ss[x][0] = y;
if(y > ss[x][1])
ss[x][1] = y;
}
}
int ans1 = 0, tmp = 0, ans2 = 4;
t.clear();
for(p = ss.begin(); p != ss.end(); p++) {
if(t.empty()) {
ans1 += p->second[1] - p->second[0];
}
t.push_back(point(p->first, p->second[0] * -1));
tmp = p->second[1] - p->second[0];
}
ans1 += tmp;
judge(ans1, ans2);
t.clear();
for(p = ss.begin(); p != ss.end(); p++) {
t.push_back(point(p->first, p->second[1]));
}
judge(ans1, ans2);
printf("%lf\n", ans1 + sqrt(2.0 * ans2 * ans2));
}
return 0;
}