BZOJ_1295_[SCOI2009]最長距離_dij

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

Sample Output

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很小，我們跑一遍最短路處理出任意兩點之間最少經過多少障礙。

然後O(n^4)枚舉兩個點，判斷兩點距離是否小於等於T，更新答案即可。

代碼：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <ext/pb_ds/priority_queue.hpp>
using namespace std;
using namespace __gnu_pbds;
#define p(i,j) ((i-1)*m+j)
__gnu_pbds::priority_queue<pair<int,int> >q;
int head[1050],to[200050],nxt[200050],val[200050],cnt,ans;
int dis[1050][1050],n,m,T,map[35][35],vis[1050];
inline void add(int u,int v,int w) {
	to[++cnt]=v; nxt[cnt]=head[u]; head[u]=cnt; val[cnt]=w;
}
int ds(int x,int y,int z,int w) {
	return (x-z)*(x-z)+(y-w)*(y-w);
}
void dij(int s) {
	dis[s][s]=0; q.push(make_pair(0,s));
	memset(vis,0,sizeof(vis));
	while(!q.empty()) {
		int x=q.top().second,i; q.pop();
		if(vis[x]) continue;
		vis[x]=1;
		for(i=head[x];i;i=nxt[i]) {
			if(dis[s][to[i]]>dis[s][x]+val[i]) {
				dis[s][to[i]]=dis[s][x]+val[i];
				q.push(make_pair(-dis[s][to[i]],to[i]));
			}
		}
	}
}
int main() {
	// freopen("distance.in","r",stdin);
	// freopen("distance.out","w",stdout);
	scanf("%d%d%d",&n,&m,&T);
	int i,j,k,l;
	for(i=1;i<=n;i++) {
		for(j=1;j<=m;j++) {
			scanf("%1d",&map[i][j]);
			if(i>1) add(p(i-1,j),p(i,j),map[i][j]),add(p(i,j),p(i-1,j),map[i-1][j]);
			if(j>1) add(p(i,j-1),p(i,j),map[i][j]),add(p(i,j),p(i,j-1),map[i][j-1]);
		}
	}
	memset(dis,0x3f,sizeof(dis));
	for(i=1;i<=n;i++) {
		for(j=1;j<=m;j++) {
			dij(p(i,j));
		}
	}
	for(i=1;i<=n;i++) {
		for(j=1;j<=m;j++) {
			for(k=1;k<=n;k++) {
				for(l=1;l<=m;l++) {
					if(dis[p(i,j)][p(k,l)]+map[i][j]<=T) ans=max(ans,ds(i,j,k,l));
				}
			}
		}
	}
	printf("%.6f\n",1.0*sqrt(1.0*ans));
}

BZOJ_1295_[SCOI2009]最長距離_dij