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

Sample Output

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;
}

BZOJ1295 [SCOI2009]最長距離