HDU 5046 Airport

HDU 5046 Airport

题目链接

题意：给定一些机场，要求选出K个机场，使得其他机场到其他机场的最大值最小

思路：二分+DLX重复覆盖去判断即可

代码：

1. #include <cstdio>  
2. #include <cstring>  
3.   
4. using namespace std;  
5.   
6. const int MAXNODE = 4005;  
7. const int MAXM = 65;  
8. const int MAXN = 65;  
9. const int INF = 0x3f3f3f3f;  
10.   
11. int K;  
12.   
13. struct DLX {  
14.   
15.     int n, m, size;  
16.   
17.     int U[MAXNODE], D[MAXNODE], R[MAXNODE], L[MAXNODE], row[MAXNODE], col[MAXNODE];  
18.     int H[MAXN], S[MAXM];  
19.     int ansd, ans[MAXN];  
20.   
21.     void init(int n, int m) {  
22.         this->n = n;  
23.         this->m = m;  
24.         ansd = INF;  
25.         for(int i = 0; i <= m; i++) {  
26.             S[i] = 0;  
27.             U[i] = D[i] = i;  
28.             L[i] = i - 1;  
29.             R[i] = i + 1;  
30.         }  
31.         R[m] = 0; L[0] = m;   
32.         size = m;  
33.         for(int i = 1; i <= n; i++)  
34.             H[i] = -1;  
35.     }  
36.   
37.     void Link(int r, int c) {  
38.         ++S[col[++size] = c];  
39.         row[size] = r;  
40.         D[size] = D[c];  
41.         U[D[c]] = size;  
42.         U[size] = c;  
43.         D[c] = size;  
44.         if(H[r] < 0) H[r] = L[size] = R[size] = size;  
45.         else {  
46.             R[size] = R[H[r]];  
47.             L[R[H[r]]] = size;  
48.             L[size] = H[r];  
49.             R[H[r]] = size;  
50.         }  
51.     }  
52.   
53.     void remove(int c) {  
54.         for(int i = D[c]; i != c; i = D[i]) {  
55.             L[R[i]] = L[i];  
56.             R[L[i]] = R[i];  
57.         }  
58.     }  
59.   
60.     void resume(int c) {  
61.         for(int i = U[c]; i != c; i = U[i])  
62.             L[R[i]] = R[L[i]] = i;  
63.     }  
64.   
65.     bool v[MAXNODE];  
66.   
67.     int f() {  
68.         int ret = 0;  
69.         for(int c = R[0]; c != 0; c = R[c]) v[c] = true;  
70.         for(int c = R[0]; c != 0; c = R[c]) {  
71.             if(v[c]) {  
72.                 ret++;  
73.                 v[c] = false;  
74.                 for(int i = D[c]; i != c; i = D[i])  
75.                     for(int j = R[i]; j != i; j = R[j])  
76.                         v[col[j]] = false;  
77.             }  
78.         }  
79.         return ret;  
80.     }  
81.   
82.     bool Dance(int d) {  
83.         if(d + f() > K) return false;  
84.         if(R[0] == 0) {  
85.             return d <= K;  
86.         }  
87.         int c = R[0];  
88.         for(int i = R[0]; i != 0; i = R[i]) {  
89.             if(S[i] < S[c])  
90.                 c = i;  
91.         }  
92.         for(int i = D[c]; i != c; i = D[i]) {  
93.             remove(i);  
94.             for(int j = R[i]; j != i; j = R[j]) remove(j);  
95.             ans[d] = row[i];  
96.             if (Dance(d + 1)) return true;  
97.             for(int j = L[i]; j != i; j = L[j]) resume(j);  
98.             resume(i);  
99.         }  
100.         return false;  
101.     }  
102. } gao;  
103.   
104. typedef long long ll;  
105.   
106. int T, n;  
107.   
108. const int N = 65;  
109.   
110. struct City {  
111.     ll x, y;  
112.     void read() {  
113.         scanf("%I64d%I64d", &x, &y);  
114.     }  
115. } c[N];  
116.   
117. ll dis(City a, City b) {  
118.     ll dx = a.x - b.x; if (dx < 0) dx = -dx;  
119.     ll dy = a.y - b.y; if (dy < 0) dy = -dy;  
120.     return dx + dy;  
121. }  
122.   
123. int main() {  
124.     int cas = 0;  
125.     scanf("%d", &T);  
126.     while (T--) {  
127.         scanf("%d%d", &n, &K);  
128.         for (int i = 1; i <= n; i++) c[i].read();  
129.         ll l = 0, r = 100000000000LL;  
130.         while (l < r) {  
131.             ll mid = (l + r) / 2;  
132.             gao.init(n, n);  
133.             for (int i = 1; i <= n; i++) {  
134.                 for (int j = 1; j <= n; j++) {  
135.                     if (dis(c[i], c[j]) <= mid)  
136.                         gao.Link(i, j);  
137.                 }  
138.             }  
139.             if (gao.Dance(0)) r = mid;  
140.             else l = mid + 1;  
141.         }  
142.         printf("Case #%d: %I64d\n", ++cas, l);  
143.     }  
144.     return 0;  
145. }