SZTUOJ 1010.最优对称路径

Description

给一个n行n列的网格，每个格子里有一个1到9的数字。你需要从左上角走到右下角，其中每一步只能往上、下、左、右四个方向之一走到相邻格子，不能斜着走，也不能走出网格，但可以重复经过一个格子。为了美观，你经过的路径还必须关于“左下-右上”这条对角线对称。

你的任务是统计所有合法路径中，数字之和最小的路径有多少条。

Input

输入最多包含25组测试数据。每组数据第一行为一个整数n（2<=n<=200）。以下n行每行包含n个1到9的数字，表示输入网格。输入结束标志为n=0。

Output

对于每组数据，输出合法路径中，数字之和最小的路径条数除以1,000,000,009的余数。

Sample Input

2
1 1
1 1
3
1 1 1
1 1 1
2 1 1
0

Sample Output

2
3

Hint

Source

湖南省第七届大学生计算机程序设计竞赛

题目解析

乍一看可能像搜索题。但是确实不好锁定状态。分析题目。首先发现对称性可以把图折叠成只剩左上半边。把数字叠一下就行。然后发现要求最短路。对于这种自带边权的网格图。我们直接bfs求出出发点到所有点的最短路。然后选择对称轴上最小的点做起点。向前利用最短路的性质做一遍记忆化搜索。dp设计状态为dp[x] [y]代表的是代表x，y到1，1上的最短路有dp[x] [y]条，且dp[0] [0]=1.非常经典的一个最短路结合dp的题目。spfa也可以求出最短路。这里为了简单就用bfs直接求了。具体实现细节可以参考代码

代码实现

#include<bits/stdc++.h>
using namespace std;
const int N=210,mod=1000000009;
int a[N][N],dp[N][N],dis[N][N],n,minn;
int dir[4][2]={-1,0,0,1,1,0,0,-1};
struct node{
	int x,y;
};
void bfs()
{
	queue<node>q;
	q.push({1,1});
	int x,y,x1,y1;
	dis[1][1]=a[1][1];
	while(!q.empty())
	{
		node tmp=q.front();
		q.pop();
		x=tmp.x;
		y=tmp.y;
		for(int i=0;i<4;i++)
		{
			x1=x+dir[i][0];
			y1=y+dir[i][1];
			if(a[x1][y1]&&dis[x1][y1]>dis[x][y]+a[x1][y1])//用a[x1][y1]判断边界。 
			{
				dis[x1][y1]=dis[x][y]+a[x1][y1];
				q.push({x1,y1});
			}
		}
	}
	minn=dis[1][n];
	for(int i=2;i<n;i++)
	minn=min(minn,dis[i][n-i+1]);
}
int dfs(int x,int y)
{
	int i,dx,dy;
	if(x==1&&y==1)return 1;
	if(dp[x][y])return dp[x][y];
	for(int i=0;i<4;i++)
	{
		dx=x+dir[i][0];
		dy=y+dir[i][1];
		if(a[dx][dy]&&dis[dx][dy]==dis[x][y]-a[x][y])
		dp[x][y]=(dp[x][y]+dfs(dx,dy))%mod;
	}
	return dp[x][y];
} 
int main()
{
	while(cin>>n)
	{
		if(n==0)break;
		memset(a,0,sizeof(a));
		memset(dp,0,sizeof(dp));
		for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		{
			cin>>a[i][j];
			dis[i][j] = mod+ 2;
		}
		for(int i=1;i<=n;i++)
		for(int j=1;j<=n-i;j++)//对称性重新建图 
		{
			a[i][j]+=a[n-j+1][n-i+1];
			a[n-j+1][n-i+1]=0;
		}
		bfs();
		int ans=0;
		for(int i=1;i<=n;i++)
		if(dis[i][n-i+1]==minn)
		ans=(ans+dfs(i,n-i+1))%mod;
		cout<<ans<<endl; 
	}
	return 0;
}