曼哈顿最小生成树

曼哈顿最小生成树

poj 3241

学习博客
注意\(k<n\), 因此可以可以通过减少计算的边来降低复杂度。

查找\(y_j-x_j>y_i-x_i\),且\(y_j+x_j\)最小的j点
具体为什么可以这样,可以看上面博客的证明

关于为什么只要四种坐标变化

只用看x的正半轴就可以了,因为是对称的。

ac code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>
#include <vector>
using namespace std;
const int maxn = 1e6+10;
int INF = 0x3f3f3f3f;
int n, k;
struct Point{
int x, y;
int id;
}p[maxn];
struct BIT{
int min_val, pos;
void init(){
min_val = INF;
pos = -1;
}
}bit[maxn];
struct Edge{
int u, v, d;
}edge[maxn<<2];
bool cmp_edge(Edge a,Edge b){
return a.d<b.d;
}
int tot;
int fa[maxn];
int find_fa(int x){
int a = x;
while(x!=fa[x]) x = fa[x];
while(a!=fa[a]){
int z = a;
a = fa[a];
fa[z] = x;
}
return x;
}
void add_edge(int u, int v, int w){
edge[tot].v = v;
edge[tot].u = u;
edge[tot++].d = w;
}
int low_bit(int x){
return x&(-x);
}
//i表示从[1, i]区间开始的
void update(int i, int val, int pos){
while(i>0){
if(val<bit[i].min_val){
bit[i].min_val = val;
bit[i].pos = pos;
}
i -= low_bit(i);
}
}
int ask(int i, int m){
int min_val = INF, pos = -1;
while(i<=m){
if(bit[i].min_val<min_val){
min_val = bit[i].min_val;
pos = bit[i].pos;
}
i += low_bit(i);
}
return pos;
}
bool cmp(Point a, Point b){
if(a.x!=b.x) return a.x<b.x;
else return a.y<b.y;
}
bool cmpedge(Edge a, Edge b){
return a.d<b.d;
}
int dis(Point a, Point b){
return abs(a.x-b.x)+abs(a.y-b.y);
}
//数组开在里面会爆栈
int a[maxn], b[maxn];
void Manha_minimum_spanning_tree(int n, int k){
int tot = 0;
for(int dir=0; dir<4; dir++){
if(dir == 1|| dir == 3){
for(int i=0; i<n; i++){
swap(p[i].x, p[i].y);
}
}
else if(dir == 2){
for(int i=0; i<n; i++){
p[i].x = -p[i].x;
}
}
sort(p, p+n, cmp);
for(int i=0; i<n; i++){
a[i] = b[i] = p[i].y-p[i].x;
}
sort(b, b+n);
int m = unique(b, b+n)-b;
for(int i=1; i<=m; i++){
bit[i].init();
}
for(int i=n-1; i>=0; i--){
int pos = lower_bound(b, b+m, a[i]) - b+1;//坐标从1开始
int ans = ask(pos, m);
if(ans != -1){
add_edge(p[i].id, p[ans].id, dis(p[i], p[ans]));
}
update(pos, p[i].x+p[i].y, i);
}
}
}
int solve(int k){
Manha_minimum_spanning_tree(n, k);
for(int i=0; i<=n; i++)fa[i] = i;
sort(edge, edge+tot, cmpedge);
for(int i=0; i<tot; i++){
int u = edge[i].u;
int v = edge[i].v;
int d = edge[i].d;
int fau = find_fa(u);
int fav = find_fa(v);
if(fau!=fav){
fa[fau] = fav;
k--;
if(k ==0 ) return edge[i].d;
}
}
}
int main(){
ios::sync_with_stdio(false);
while(~scanf("%d%d", &n, &k)&&n){
for(int i=0; i<n; i++){
scanf("%d%d", &p[i].x, &p[i].y);
p[i].id = i;
}
printf("%d\n", solve(n-k));
}
return 0;
}

未解决的问题

文章目录
  1. 1. poj 3241
    1. 1.0.1. 关于为什么只要四种坐标变化
  2. 1.1. ac code
  • 2. 未解决的问题
  • {{ live2d() }}