[BZOJ1295][SCOI2009]最长距离 最短路+枚举
1295: [SCOI2009]最长距离
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 1683 Solved: 912
[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
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 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<cstdlib> 5 #include<cmath> 6 #include<algorithm> 7 #include<queue> 8 using namespace std; 9 const int tx[5]={0,0,1,-1}; 10 const int ty[5]={1,-1,0,0}; 11 struct data { 12 int x,y,val; 13 bool operator <(const data &t)const { 14 return val>t.val; 15 } 16 }; 17 int n,m,T; 18 char s[101][101]; 19 bool vis[35][35]; 20 int dis[35][35]; 21 void dij(int stx,int sty) { 22 priority_queue<data> q; 23 memset(vis,0,sizeof(vis)); 24 memset(dis,97,sizeof(dis)); 25 dis[stx][sty]=s[stx][sty]=='1'?1:0; 26 q.push((data){stx,sty,dis[stx][sty]}); 27 while(!q.empty()) { 28 data now=q.top();q.pop(); 29 if(vis[now.x][now.y]) continue; 30 31 vis[now.x][now.y]=1; 32 for(int i=0;i<4;i++) { 33 int tox=now.x+tx[i],toy=now.y+ty[i]; 34 if(tox>n||toy>m||tox<1||toy<1) continue; 35 if(dis[tox][toy]>now.val+(s[tox][toy]=='1'?1:0)) { 36 dis[tox][toy]=now.val+(s[tox][toy]=='1'?1:0); 37 q.push((data){tox,toy,dis[tox][toy]}); 38 } 39 } 40 } 41 } 42 int main() { 43 int ans=0; 44 scanf("%d%d%d",&n,&m,&T); 45 for(int i=1;i<=n;i++) scanf("%s",s[i]+1); 46 for(int i=1;i<=n;i++) { 47 for(int j=1;j<=m;j++) { 48 dij(i,j); 49 int mx=0; 50 for(int k=1;k<=n;k++) 51 for(int w=1;w<=m;w++) 52 if(dis[k][w]<=T) mx=max(mx,(i-k)*(i-k)+(j-w)*(j-w)); 53 ans=max(ans,mx); 54 } 55 } 56 double tmp=sqrt(ans); 57 printf("%.6lf",tmp); 58 }
O(∩_∩)O~ (*^__^*) 嘻嘻…… O(∩_∩)O哈哈~