295: [SCOI2009]最长距离
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 1455 Solved: 789
[ 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
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
1.414214
【输出样例二】
3.605551
【输出样例三】
2.828427
把每个节点是否为1,看做从一个节点到另一个节点的代价,这样对每一个点做一次最短路,在暴力更新答案就好了。
#include<cstdio>
#include<cstring>
#include<cmath>
#define max(x,y) ((x)>(y)?(x):(y))
using namespace std;
int n,m,t;
char mp[75][75];
int d[75][75],l,r;
double ans;
bool vis[75][75];
struct p{
int x,y;
p(){}
p(int a,int b):x(a),y(b){}
}q[20000];
void spfa(int x,int y)
{
memset(d,0x3f,sizeof(d));
memset(vis,0,sizeof(vis));
vis[x][y]=1;l=r=0;
d[x][y]=(mp[x][y]=='1');
q[r++]=p(x,y);
while(l<r)
{
p t=q[l++];
for(int i=-1;i<=1;i++)
for(int j=-1;j<=1;j++)
if((i==0||j==0)&&t.x+i>=0&&t.x+i<n&&t.y+j>=0&&t.y+j<m&&d[t.x][t.y]+(mp[t.x+i][t.y+j]=='1')<d[t.x+i][t.y+j])
{
d[t.x+i][t.y+j]=d[t.x][t.y]+(mp[t.x+i][t.y+j]=='1');
if(vis[t.x+i][t.y+j]==false){
vis[t.x+i][t.y+j]=true;
q[r++]=p(t.x+i,t.y+j);
}
}
vis[t.x][t.y]=0;
}
}
int main(){
scanf("%d%d%d",&n,&m,&t);
for(int i=0;i<n;i++)
scanf("%s",mp[i]);
for(int i=0;i<n;i++)
for(int j=0;j<m;j++){
spfa(i,j);
for(int ii=0;ii<n;ii++)
for(int jj=0;jj<m;jj++)
if(d[ii][jj]<=t){
ans = max(ans,sqrt( (ii-i)*(ii-i) + (jj-j)*(jj-j) ) );
}
}
printf("%.6f",ans);
return 0;
}