bzoj 1295 [SCOI2009]最长距离最短路

[SCOI2009]最长距离

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 1715 Solved: 932
[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 #include<cstring>
 2 #include<cmath>
 3 #include<algorithm>
 4 #include<iostream>
 5 #include<cstdio>
 6 
 7 using namespace std;
 8 inline int read()
 9 {
10     int x=0,f=1;char ch=getchar();
11     while(ch>'9'||ch<'0'){if (ch=='-') f=-1;ch=getchar();}
12     while(ch<='9'&&ch>='0'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
13     return x*f;
14 }
15 
16 int n,m,T;
17 double ans=0;
18 bool mp[31][31],inq[31][31],vis[31][31];
19 int dis[31][31];
20 int dx[4]={1,-1,0,0},dy[4]={0,0,-1,1};
21 struct data
22 {
23     int x,y;
24 }q[100007];
25 
26 void getans(int x,int y)
27 {
28     for(int i=x;i<=n;i++)
29         for(int j=1;j<=m;j++)
30             if(dis[i][j]<=T&&(y-j)*(y-j)+(x-i)*(x-i)>ans)
31                 ans=(y-j)*(y-j)+(x-i)*(x-i);
32 }
33 void search(int xx,int yy)
34 {
35     int nowx,nowy,i,t=1,w=1,nx,ny;
36     q[1].x=xx;q[1].y=yy;
37     memset(inq,0,sizeof(inq));
38     memset(dis,127,sizeof(dis));
39     inq[xx][yy]=1;dis[xx][yy]=mp[xx][yy];
40     while(t<=w)
41     {
42         nowx=q[t].x;nowy=q[t].y;
43         t++;
44         for(i=0;i<4;i++)
45         {
46             nx=nowx+dx[i];ny=nowy+dy[i];
47             if(nx>n||nx<xx||ny>m||ny<1)continue;
48             if(dis[nowx][nowy]+mp[nx][ny]<dis[nx][ny])
49             {
50                  dis[nx][ny]=dis[nowx][nowy]+mp[nx][ny];
51                  if(!inq[nx][ny]) q[++w].x=nx,q[w].y=ny,inq[nx][ny]=1;
52             }
53         }
54         inq[nowx][nowy]=0;
55     }
56     getans(xx,yy);
57 }
58 int main()
59 {
60     n=read(),m=read(),T=read();
61     char str[40];
62     for(int i=1;i<=n;i++)
63     {
64         scanf("%s",str);
65         for(int j=0;j<m;j++)
66         mp[i][j+1]=str[j]-'0';
67     }
68     for(int i=1;i<=n;i++)
69         for(int j=1;j<=m;j++)
70             search(i,j);
71     printf("%.6lf",sqrt(ans));
72 }