[bzoj1295]: [SCOI2009]最长距离

1295: [SCOI2009]最长距离

Time Limit: 10 Sec   Memory Limit: 162 MB
Submit: 1576   Solved: 853
[ Submit][ Status][ Discuss]

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 。

题解

用spfa 因为n,m最大30 所以可以枚举每个起点 求出移走T个障碍的最大距离

spfa中的做法就是用dis[i][j]表示从枚举的起点到达(i,j)最少移走的障碍物 

两种情况讨论 i,j是否为障碍物 是就+1

最后对于小于T的点 算出他和起点的距离找最大的即可

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
const int N=35;
int n,m,T;char s[N];
struct node{
	int x,y;
	node(int a,int b){
		x=a;y=b;
	}
};
queue<node>q;
int flag[N][N];
int dis[N][N],a[N][N];
int dx[]={0,1,-1,0,0};
int dy[]={0,0,0,1,-1};
double ans;
void getans(int x,int y){
	for(int i=1;i<=n;i++)
	   for(int j=1;j<=m;j++)
	      if(dis[i][j]<=T&&sqrt((y-j)*(y-j)+(x-i)*(x-i))>ans)
	         ans=sqrt((y-j)*(y-j)+(x-i)*(x-i));//求最大
}
void spfa(int x,int y){//dis[i][j]=k表示从(x,y)出发移走k个障碍才可以到达(i,j)
	q.push(node(x,y));
	memset(dis,127,sizeof(dis));
	memset(flag,0,sizeof(flag));
	dis[x][y]=a[x][y];
	flag[x][y]=1;
	while(!q.empty()){
		int kx=q.front().x;int ky=q.front().y;q.pop();flag[kx][ky]=0;
		for(int i=1;i<=4;i++)
		{
			int tx=kx+dx[i];int ty=ky+dy[i];
			if(tx>n||tx<=0||ty>m||ty<=0) continue;
			if(a[tx][ty]&&dis[kx][ky]+1<dis[tx][ty]){//如果是障碍
				dis[tx][ty]=dis[kx][ky]+1;//如果从tx,ty这条路上移走这个障碍带来的d[tx][ty]的值可以松弛 就松弛它
				if(!flag[tx][ty]){
					flag[tx][ty]=1;
					q.push(node(tx,ty));
				}
			}
			else if(!a[tx][ty]&&dis[kx][ky]<dis[tx][ty]){//不是障碍
				dis[tx][ty]=dis[kx][ky];
				if(!flag[tx][ty]){
					flag[tx][ty]=1;
					q.push(node(tx,ty));
				}
			}
		}
	}
	getans(x,y);
}
int main(){
	scanf("%d%d%d",&n,&m,&T);
	for(int i=1;i<=n;i++)
		{
			scanf("%s",s+1);
			for(int j=1;j<=m;j++)
				a[i][j]=s[j]-'0';
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			spfa(i,j);
	printf("%.6lf",ans);
	return 0;
}