【BZOJ】1295: [SCOI2009]最长距离（spfa+暴力）

http://www.lydsy.com/JudgeOnline/problem.php?id=1295

咳咳。。此题我不会做啊。。一开始认为是多源，可是有移除物品的操作，所以不行。

此题的思想很巧妙！

我们不妨将问题转换一下，对于一个点到另一个点，我们只需算出到达这个点最少需要移除多少个障碍，然后用题目给的障碍判断是否可行，然后暴力算出可行的点之间的欧几里得距离就行了orz。

T_T

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define read(a) a=getnum()
#define print(a) printf("%d", a)
inline int getnum() { int ret=0; char c; for(c=getchar(); c<'0' || c>'9'; c=getchar()); for(; c>='0' && c<='9'; c=getchar()) ret=ret*10+c-'0'; return ret; }

const int N=35;
int d[N][N], a[N][N], vis[N][N], ans, n, m;
int dx[]={-1, 1, 0, 0}, dy[]={0, 0, -1, +1};
struct qu{ int x, y; } q[N*N*2];
int front, tail;

void spfa(const int &s, const int &t) {
	int x, y;
	CC(vis, 0); CC(d, 0x7f);
	d[s][t]=front=tail=0;
	vis[s][t]=1; q[tail].x=s, q[tail++].y=t;
	while(front<tail) {
		x=q[front].x; y=q[front++].y;
		rep(i, 4) if(x+dx[i] && x+dx[i]<=n && y+dy[i] && y+dy[i]<=m && d[x][y]+a[x+dx[i]][y+dy[i]]<d[x+dx[i]][y+dy[i]]) {
			d[x+dx[i]][y+dy[i]]=d[x][y]+a[x+dx[i]][y+dy[i]];
			if(!vis[x+dx[i]][y+dy[i]]) {
				vis[x+dx[i]][y+dy[i]]=1;
				q[tail].x=x+dx[i]; q[tail++].y=y+dy[i];
			}
		}
		vis[x][y]=0;
	}
}

int main() {
	read(n); read(m);
	int t=getnum();
	char c;
	for1(i, 1, n) for1(j, 1, m) {
		for(c=getchar(); c<'0' || c>'9'; c=getchar()); 
		a[i][j]=c-'0';
	}
	int t1;
	for1(i, 1, n) for1(j, 1, m) {
		spfa(i, j);
		t1=t; if(a[i][j]) --t1;
		for1(x, 1, n) for1(y, 1, m)	if(d[x][y]<=t1)
			ans=max(ans, (i-x)*(i-x)+(j-y)*(j-y));
	}
	printf("%.6lf", sqrt((double)ans));
	return 0;
}

---

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

Sample Output

【输出样例一】
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 。

Source

Day2