穿越矩阵

穿越矩阵

现在有一个 m * n 的整数矩阵，每一格上有一个数字。请你编写一个程序计算出一条从左到右穿过矩阵的路径，并使此路径的费用最小。路径的花费是指这条路径所穿越的所有单元格中的数字之和。

路径从矩阵的左侧的第一列的任意单元格开始，逐步穿过矩阵到达最右侧的一列的任意单元格。每一步是指从某单元格进入它右侧一列的相邻单元格（如下图，可以是横向或斜向）。矩阵的第一行和最后一行实际是相邻的，你可以想象矩阵是包裹在一个横放的圆柱体外面。

穿越两个略有不同的 5 * 6 的矩阵的路径如下图所示，这两个矩阵只有最后一行的数字不同。右侧的路径显示了第一行和最后一行相邻的效果。

输入

输入包括一系列矩阵描述。每个矩阵描述的第一行是 m 和 n，即矩阵的行数和列数；之后的 m 行，每行包括 n 个以空格分开的整数，则是当前矩阵的值，注意矩阵的值未必是正数。

矩阵的行数 m 和列数 n 的范围是：1 <= m <= 10、 1 <= n <= 100；所有路径的费用值都可以用 30bit 的整数表示。

输出

针对每一个矩阵，找出费用最小的路径，并将其输出。每个矩阵的输出包括两行，第一行是路径本身，即输出每一步所在的行，第二行则是该路径的费用。

如果对于同一个矩阵有多条不同的费用最小路径，则输出左端行号较小的一条。

Sample input

5 6

3 4 1 2 8 6

6 1 8 2 7 4

5 9 3 9 9 5

8 4 1 3 2 6

3 7 2 8 6 4

5 6

3 4 1 2 8 6

6 1 8 2 7 4

5 9 3 9 9 5

8 4 1 3 2 6

3 7 2 1 2 3

Sample Output

1 2 3 4 4 5

16

1 2 1 5 4 5

11

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
int a[11][101],dp[11][101];
int num[101];

int main()
{
	int m,n;
	while(scanf("%d%d",&m,&n)!=EOF)
	{
		for(int i=0;i<m;i++)
			for(int j=0;j<n;j++)
				cin>>a[i][j];
		memset(dp,0,sizeof(dp));
		for(int i=0; i<m; i++)
			dp[i][n-1] = a[i][n-1];
		for(int j=n-2; j>=0; j--)
		{
			for(int i=0; i<m; i++)
			{
				dp[i][j] = a[i][j];
				if(i>0 && i<m-1)
					dp[i][j] += min(min(dp[i][j+1],dp[i+1][j+1]),dp[i-1][j+1]);
				else if(i==0)
					dp[i][j] += min(min(dp[i][j+1],dp[i+1][j+1]),dp[m-1][j+1]);
				else if(i==m-1)
					dp[i][j] += min(min(dp[i][j+1],dp[ 0 ][j+1]),dp[i-1][j+1]);
			//	cout<<dp[i][j]<<" ";
			}
		//	cout<<endl;
		}
		int mini, mine = dp[m-1][0];
		for(int i=m-1; i>=0; i--)
		{
			if(mine >= dp[i][0])
			{
				mine = dp[i][0];
				mini = i;
			}
		}
		memset(num,0,sizeof(num));
		num[0] = mini;
		for(int j=1; j<n; j++)
		{
			int flag = 0;
			for(int i=0; i<m; i++)
			{
				if(mini==0 && i!=0 && i!=1 && i!=m-1)
					continue;
				if(mini==m-1 && i!=m-1 && i!=m-2 && i!=0)
					continue;
				if(mini>0 && mini<m-1 && i!=mini && i!=mini+1 && i!=mini-1)
					continue;
				if(dp[mini][j-1] - a[mini][j-1]==dp[i][j])
				{
					mini = i;
					num[j] =mini;
					flag =1;
				}
				if(flag)
					break;
			}
		}
		cout<<num[0]+1;
		for(int j=1;j<n;j++)
			cout<<" "<<num[j]+1;
			cout<<endl;
		cout<<mine<<endl;
	}
	return 0;
}