Description
windy有一块矩形土地,被分为 N*M 块 1*1 的小格子。 有的格子含有障碍物。 如果从格子A可以走到格子B,那么两个格子的距离就为两个格子中心的欧几里德距离。 如果从格子A不可以走到格子B,就没有距离。 如果格子X和格子Y有公共边,并且X和Y均不含有障碍物,就可以从X走到Y。 如果windy可以移走T块障碍物,求所有格子间的最大距离。 保证移走T块障碍物以后,至少有一个格子不含有障碍物。
Input
输入文件maxlength.in第一行包含三个整数,N M T。 接下来有N行,每行一个长度为M的字符串,'0'表示空格子,'1'表示该格子含有障碍物。
Output
输出文件maxlength.out包含一个浮点数,保留6位小数。
Sample Input
【输入样例一】
3 3 0
001
001
110
【输入样例二】
4 3 0
001
001
011
000
【输入样例三】
3 3 1
001
001
001
3 3 0
001
001
110
【输入样例二】
4 3 0
001
001
011
000
【输入样例三】
3 3 1
001
001
001
Sample Output
【输出样例一】
1.414214
【输出样例二】
3.605551
【输出样例三】
2.828427
1.414214
【输出样例二】
3.605551
【输出样例三】
2.828427
HINT
20%的数据,满足 1 <= N,M <= 30 ; 0 <= T <= 0 。 40%的数据,满足 1 <= N,M <= 30 ; 0 <= T <= 2 。 100%的数据,满足 1 <= N,M <= 30 ; 0 <= T <= 30 。
题解
枚举答案中两个格子中的一个,暴力Dijkstra,每经过一个1距离就加1,距离起点距离不超过T的就是移走至多T块障碍物能到的格子。
附代码:
#include <algorithm>
#include <cstdio>
#include <cmath>
#include <queue>
const int N = 50;
int map[N][N];
struct HeapNode{
int i, j, d;
HeapNode(int i = 0, int j = 0, int d = 0) : i(i), j(j), d(d) {}
bool operator<(const HeapNode &x)const{
return d > x.d;
}
};
std::priority_queue<HeapNode> Q;
int dis[N][N], n, m, T, ans;
const int d[][2] = {{0, 1}, {0, -1}, {1, 0}, {-1, 0}};
void Dijkstra(int si, int sj) {
std::fill(dis[0], dis[n], 10000);
while (!Q.empty()) Q.pop();
Q.push(HeapNode(si, sj, dis[si][sj] = map[si][sj]));
while (!Q.empty()) {
HeapNode x = Q.top(); Q.pop();
if (x.d != dis[x.i][x.j]) continue;
if (x.d > T) return;
int ui = x.i, uj = x.j;
ans = std::max(ans, (ui - si) * (ui - si) + (uj - sj) * (uj - sj));
for (int k = 0; k < 4; ++k) {
int vi = ui + d[k][0], vj = uj + d[k][1];
if (~vi && ~vj && vi < n && vj < m && dis[vi][vj] > x.d + map[vi][vj])
Q.push(HeapNode(vi, vj, dis[vi][vj] = x.d + map[vi][vj]));
}
}
}
int main() {
scanf("%d%d%d", &n, &m, &T);
for (int i = 0; i < n; ++i)
for (int j = 0; j < m; ++j)
scanf("%1d", &map[i][j]);
ans = 0;
for (int i = 0; i < n; ++i)
for (int j = 0; j < m; ++j)
Dijkstra(i, j);
printf("%.6f", sqrt((float)ans));
return 0;
}