【bzoj 1295】[SCOI2009]最长距离（spfa）

1295: [SCOI2009]最长距离

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 1380  Solved: 745
[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 。

Source

[Submit][Status][Discuss]

【题解】【spfa】

【无法想象这竟然是省选题。。。】

【裸spfa，每次判断到每个点的最短路，是否在经过<=T个障碍物的情况下，如果满足，最大距离就是两点之间的欧几里得距离】

【调了辣么久，因为我sb的打错了spfa】

#include<cmath>
#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
struct node{
	int hang,line;
};
int d[4][2]={{-1,0},{1,0},{0,-1},{0,1}};
int dis[40][40],a[40][40];
int n,m,T;
double ans=0;
bool p[40][40];
inline double len(double x,double y)
{
	return x*x+y*y;
}
inline void spfa(int sx,int sy)
{
	memset(dis,127/3,sizeof(dis));
	memset(p,0,sizeof(p));
	queue<node>que;
	node x; 
	x.hang=sx; x.line=sy;
	que.push(x);  p[sx][sy]=1;
	if(a[sx][sy]) dis[sx][sy]=1;
	 else dis[sx][sy]=0;
	while(!que.empty())
     {
     	node u=que.front(); que.pop();
     	int x=u.hang,y=u.line; 
     	for(int i=0;i<4;++i)
     	 {
     	 	int xx=x+d[i][0],yy=y+d[i][1];
     	 	if(xx>0&&xx<=n&&yy>0&&y<=m)
     	 	 if(dis[xx][yy]>dis[x][y]+a[xx][yy])
     	 	  {
     	 	  	dis[xx][yy]=dis[x][y]+a[xx][yy];
     	 	  	if(!p[xx][yy])
     	 	  	 {
     	 	  	 	p[xx][yy]=1;
     	 	  	 	node v;
     	 	  	 	v.hang=xx; v.line=yy;
     	 	  	 	que.push(v);
			    }
			}
	    }
	    p[x][y]=0;
	}	
}
int main()
{
//	freopen("int.txt","r",stdin);
//	freopen("my.txt","w",stdout);
	int i,j;	
	scanf("%d%d%d\n",&n,&m,&T);
	for(i=1;i<=n;++i)
	 {
	 	char s[40];
	 	scanf("%s",s+1);
	 	for(j=1;j<=m;++j) a[i][j]=s[j]-'0';
	 }
	for(i=1;i<=n;++i)
	 for(j=1;j<=m;++j)
	  {
	  	spfa(i,j);
	  	for(int h=1;h<=n;++h)
	  	 for(int l=1;l<=m;++l)
	  	  if(dis[h][l]<=T) ans=max(ans,sqrt(len((double)i-h,(double)j-l)));
	  }
	printf("%.6lf\n",ans);
	return 0;
}