ZOJ 1430 / POJ 2679 The Erythea Campaign (bfs+dijkstra)

繁，但是不难。

 

首先，建图，题中给出的是m*n的矩阵，里面包含了所有stronghold的分布，但是你需要把他转换成一个图。但是路径不能穿过堡垒，多以说两个堡垒之间的两个点不能相互到达。还有就是每个点应该保存它可以到达的点的信息，当使用dijkstra求最短路径的时候，这些点可以当邻接表使用，很方便的。

 

接下来是每个的危险值，求法是m+n-d，不会出现负值。d的意思是它到达离他最近的堡垒的距离。对于每个堡垒上的点做bfs求出所有点到他的最短距离，如果求出的距离比之前更新的距离小的话，就存那个小的。

 

#include<iostream> #include<cstring> #include<queue> const int INF = 99999999; const int MAXLIST = 7000; const int MAXMAP = 82; using namespace std; struct node { int value; //危险值 bool isedge; //是否是边界点 int order; //点的序号 int reach[4]; //可以到达的点的序号 } map[MAXMAP][MAXMAP],list[MAXLIST];//list是类似于临接表的东西 int m,n;//n行，m列 char c; int cases; int startx,starty,endx,endy; //开始和结束的坐标 int start,end; //开始和结束的点在list里面的位置 int stronghold[MAXMAP][MAXMAP]; // 堡垒的分布 int flag[MAXLIST],dist[MAXLIST]; //dijkstra时候用的 void bfs(int n) //求出到所有点到n的距离 { queue<int> q; int inq[MAXLIST],t,dist; for(int i= 0;i <MAXLIST;i++) inq[i] = 0; inq[n] = 1;q.push(n);list[n].value=0; while( !q.empty() ) { n = q.front();q.pop(); dist = list[n].value+1; for(int j=0; j < 4; j++) { t = list[n].reach[j]; if( t==0 ) continue; if( list[t].value > dist ) { list[t].value = dist; q.push(t);inq[t] = 1;} } } } void dijkstra(int a) { int now,min,v = (n+1)*(m+1); for(int i=1;i <= v;i++) {dist[i]=INF;flag[i]=0;} dist[now=a] = list[a].value;flag[a] = 1; for(int i = 1; i < v; i++) { for(int j = 0; j < 4; j++) { int to = list[now].reach[j]; if(to == 0) continue; if( dist[to] > dist[now] + list[to].value ) dist[to] = dist[now] + list[to].value; } min = INF; for(int j = 1; j <= v; j++) if( !flag[j] && dist[j] < min ) min = dist[now=j]; flag[now] = 1; } } int main(void) { cin >> cases; while( cases-- ) { cin >> n >> m; cin >> startx >> starty >> endx >> endy; starty++;endy++; start = (m+1)*startx + starty; end = (m+1)*endx + endy; memset(stronghold,0,sizeof(stronghold)); for( int i = 0; i < n; i++ ) //输入stronghold { getchar(); for( int j = 0; j < m; j++ ) stronghold[i][j] = (c=getchar()) - '0'; } for( int i = 0,order = 1; i <= n; i++ ) //将图中的order和isedge初始化 for( int j = 0; j <= m; j++ ) { map[i][j].order = order++; map[i][j].isedge = false; } for( int i = 0; i < n; i++ ) //判断某个点是不是堡垒的边界点 for( int j = 0; j < m; j++ ) if( stronghold[i][j] == 1 ) map[i][j].isedge = map[i+1][j].isedge = map[i][j+1].isedge = map[i+1][j+1].isedge = 1; for( int i = 0; i <= n; i++ ) //点map【i】【j】可以到达的点在list中的位置 for( int j = 0; j <= m; j++ ) //不可以到的为零 { map[i][j].reach[0] = i>0?map[i-1][j].order:0; map[i][j].reach[1] = j>0?map[i][j-1].order:0; map[i][j].reach[2] = i<n?map[i+1][j].order:0; map[i][j].reach[3] = j<m?map[i][j+1].order:0; } for( int i = 0; i < n; i++ ) //两个堡垒相邻的话他们中间那个路不能走 for( int j = 0; j < m; j++ ) { if( stronghold[i][j] == 1 && stronghold[i][j+1] == 1 ) map[i][j+1].reach[2] = map[i+1][j+1].reach[0] = 0; if( stronghold[i][j] == 1 && stronghold[i+1][j] == 1 ) map[i+1][j].reach[3] = map[i+1][j+1].reach[1] = 0; } for( int i = 0,order = 1; i <= n; i++ ) //把map中的东西转移到list中 for( int j = 0; j <= m; j++ ) { map[i][j].value = INF; list[order++] = map[i][j]; } for( int i = 0; i <= n; i++ ) // 每个点的到所有堡垒边界点的最短距离 for( int j= 0; j <= m; j++ ) if( map[i][j].isedge ) bfs(map[i][j].order); for(int i = 1;i <= (n+1)*(m+1); i++ ) list[i].value = m+n-list[i].value; //危险值 dijkstra(start); if (dist[end] == INF ) cout << "no solution/n"; else cout << dist[end] <<endl; } return 0; }