bzoj1295 最长距离

1295: [SCOI2009]最长距离

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

题解：

因为数据范围较小，我们直接对每个点暴力迪杰斯特拉，两个相邻点之间连边，如果边的终点有障碍物，那么权值为1，否则为0。需要注意，如果一开始的位置有障碍物，那么自己到自己的最短路长度为1，而不是为0。

两个点之间如果去掉T个点可以到达，那么它们之间最短路的长度一定小于等于T（显然，因为当长度大于T时，你就算去掉的每个障碍物都在两点之间的最短路上，两点之间还是会有障碍物把它们堵住）。

代码：

#include<bits/stdc++.h>
using namespace std;
int n,m,e[10000],head[10000],next[10000],d[1000],g[1000],q[100000],jia[10000],b[10],c[10],t1,ans,tot,a[100][100],ii;
char s[1000];
void build(int t,int k,int s){
	tot++;
	e[tot]=k;jia[tot]=s;
	next[tot]=head[t];head[t]=tot;
}
void dj(int x,int k,int ll,int rr){
	int l=0,r=1,v,i,lll,rrr,u;
	q[1]=x;
	memset(d,63,sizeof(d));
	//printf("%d\n",d[1]);
	memset(g,0,sizeof(g));
	d[x]=k;g[x]=1;
	while(l<r){
		//printf("%d %d\n",l,r);
	    l++;
	    u=q[l];
	    v=0;
	    for(i=head[u];i;i=next[i])
	     if(d[e[i]]>d[u]+jia[i]){
	     	d[e[i]]=d[u]+jia[i];
	     	if(!g[e[i]]){
	     		r++;
	     		q[r]=e[i];
	     		g[e[i]]=1;
			 }
		 }
		g[u]=0; 
	}
	for(i=1;i<=n*m;i++)
	 if(d[i]<=t1){
	 	lll=i/m;
	 	if(i%m){
	 		rrr=i%m;
		 lll++;}
		 else{
		 	rrr=m;
		 }
		 ans=max(ans,(ll-lll)*(ll-lll)+(rr-rrr)*(rr-rrr));
	 }
}
int main(){
	b[1]=1;b[2]=-1;b[3]=0;b[4]=0;
	c[1]=0;c[2]=0;c[3]=1;c[4]=-1;
	int i,j,x,y;
	scanf("%d%d%d",&n,&m,&t1);
	for(i=1;i<=n;i++){
		scanf("%s",&s);
		for(j=0;j<m;j++)a[i][j+1]=s[j]-'0';
	}
	for(i=1;i<=n;i++)
	 for(j=1;j<=m;j++){
	 	for(ii=1;ii<=4;ii++){
	 		x=i+b[ii];y=j+c[ii];
	 		if(x<1||y<1||x>n||y>m)continue;
	 		if(a[x][y])build((i-1)*m+j,(x-1)*m+y,1);
	 		 else build((i-1)*m+j,(x-1)*m+y,0);
		 }
	 }
	for(i=1;i<=n;i++)
	 for(j=1;j<=m;j++){
	 	dj((i-1)*m+j,a[i][j],i,j);
	 } 
	 printf("%.6lf",sqrt((double)ans));
}