BZOJ 1295 Scoi2009 最长距离

Posted 2020-10-15 zhangenming

1295: [SCOI2009]最长距离

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

很水的一道最短路问题

由于n和m很小，我们可以暴力对于每一个点求出它到其他点所需要最小移除的障碍

然后如果到这个点所需要的移除障碍小于T的化，就用他们之间的距离来更新答案

 1 #include <bits/stdc++.h>
 2 #define ll long long
 3 using namespace std;
 4 inline int read(){
 5     int x=0;int f=1;char ch=getchar();
 6     while(!isdigit(ch)) {if(ch==‘-‘) f=-1;ch=getchar();}
 7     while(isdigit(ch)) {x=x*10+ch-‘0‘;ch=getchar();}
 8     return x*f;
 9 }
10 const int MAXN=1e6+10;
11 const int dx[4]={1,0,-1,0};
12 const int dy[4]={0,-1,0,1};
13 namespace zhangenming{
14     struct node{
15         int x,y;
16     }q[MAXN];
17     int head,tail,n,m,mapp[40][40],T,f[40][40],vis[40][40];
18     double ans;
19     inline void spfa(int xx,int yy){
20         head=tail=0;memset(vis,0,sizeof(vis));
21         memset(f,10,sizeof(f));
22         vis[xx][yy]=1;
23         f[xx][yy]=mapp[xx][yy];
24         q[++tail].x=xx;q[tail].y=yy;
25         while(head<tail){
26             node tn=q[++head];
27             for(int i=0;i<4;i++){
28                 int x=tn.x+dx[i];int y=tn.y+dy[i];
29                 if(x==0||y==0||x==n+1||y==m+1) continue;
30                 if(f[tn.x][tn.y]+mapp[x][y]<f[x][y]&&f[tn.x][tn.y]<=T){
31                     f[x][y]=f[tn.x][tn.y]+mapp[x][y];
32                     if(!vis[x][y]){
33                         q[++tail].x=x;q[tail].y=y;
34                         vis[x][y]=1;
35                     }
36                 }
37             }
38             vis[tn.x][tn.y]=0;
39         }
40         for(int i=1;i<=n;i++){
41             for(int j=1;j<=m;j++){
42                 if(f[i][j]<=T){
43                     ans=max(ans,sqrt(1.0*((i-xx)*(i-xx)+(j-yy)*(j-yy))));
44                 }
45             }
46         }
47     }
48     void init(){
49         n=read();m=read();T=read();
50         for(int i=1;i<=n;i++){
51             char ch[40];
52             scanf("%s",ch);
53             for(int j=0;j<m;j++){
54                 if(ch[j]==‘0‘) mapp[i][j+1]=0;
55                 else mapp[i][j+1]=1;
56             }
57         }
58     }
59     void solve(){
60         for(int i=1;i<=n;i++){
61             for(int j=1;j<=m;j++){
62                 spfa(i,j);
63             }
64         }
65         cout<<setiosflags(ios::fixed)<<setprecision(6);
66         cout<<ans<<endl;
67     }
68 }
69 int main(){
70     //freopen("All.in","r",stdin);
71     //freopen("zhang.out","w",stdout);
72     using namespace zhangenming;
73     init();
74     solve();
75     return 0;
76 }

代码

 1 #include <bits/stdc++.h>
 2 #define ll long long
 3 using namespace std;
 4 inline int read(){
 5     int x=0;int f=1;char ch=getchar();
 6     while(!isdigit(ch)) {if(ch==‘-‘) f=-1;ch=getchar();}
 7     while(isdigit(ch)) {x=x*10+ch-‘0‘;ch=getchar();}
 8     return x*f;
 9 }
10 int main(){
11     freopen("All.in","w",stdout);
12     srand(int(time(NULL)));
13     int n=rand()%29+1;int m=rand()%29+1;int T=rand()%107;
14     cout<<n<<‘ ‘<<m<<‘ ‘<<T<<endl;
15     for(int i=1;i<=n;i++){
16         for(int j=1;j<=m;j++){
17             printf("%d",rand()%2);
18         }
19         printf("\n");
20     }
21     return 0;
22 }

对拍代码