穿越矩阵

问题描述

      现在有一个 m * n 的整数矩阵，请你编写一个程序计算出一条从左到右穿过矩阵的路径，并使此路径的费用最小。路径从矩阵的左侧的第一列的任意单元格开始，逐步穿过矩阵到达最右侧的一列的任意单元格。每一步是指从某单元格进入它一列的相邻单元格（如下图，可以是横向或斜向）。矩阵的第一行和最后一行实际是相邻的，你可以想象矩阵是包裹在一个横放的圆柱体外面。

路径的花费是指这条路径所穿越的所有单元格中的数字之和。

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

输入

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

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

输出

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

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

来源

UVa: 116

	测试输入	期待的输出	时间限制	内存限制	额外进程
测试用例 1	以文本方式显示 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↵	以文本方式显示 1 2 3 4 4 5↵ 16↵ 1 2 1 5 4 5↵ 11↵	1秒	1024KB	0

题解思路

大致思路：

每次找最小的走，这样最后就可以找到最短的路径了，然后保存路径最后动归输出就好了。

具体实现：

可以从倒数第二列向前递推，也可以从前往后递归。

递推的方式：可以将这个矩阵竖着放置，这样就是从最下面向上找最小的路径和，(即从后向前推)那么即是：

    运用两重循环，倒数第二列开始向前找。

dp[i][j]=map[i][j]+min{dp[i][j+1],dp[(i+1)%m][j+1],dp[(i-1+m)][j+1]}

递归的方式：即从前向后递归找出。

运用一次循环，对于每一次计算调用DP(i,0).

dp[i][j]=map[i][j]+min{dp[i][j+1],dp[(i+1)%m][j+1],dp[(i-1+m)][j+1]}

注意：

（1）考虑只有一行或是只有一列的情况。

因为数据是多组输入的，所以每一次需要对路径的多少初始化。

实现代码

<span style="font-family:Microsoft YaHei;font-size:14px;">#include<stdio.h>
#include<string.h>
#define M 1500000000
int main()
{
	int i,j,hangshu;
	long dp[11][105];
	long temp1,temp2,temp3,temp,min;
	int lujing[11][105],map[11][105];
	int m,n,k;
	while(scanf("%d%d",&m,&n)!=EOF)
	{
		memset(dp,0,sizeof(dp));
		memset(lujing,0,sizeof(lujing));
		memset(map,0,sizeof(map));
		for(i=0;i<m;i++)
		{
			for(j=0;j<n;j++)
			{
				scanf("%d",&map[i][j]);
				if(j==n-1)
				{
					dp[i][j]=map[i][j];
				}
				else
				{
					dp[i][j]=M;
				}
			}
		}
		for(i=n-2;i>=0;i--)
		{
			for(j=0;j<m;j++)
			{
				k=M;
			 	temp1=dp[j][i+1];
				temp2=dp[(j+1)%m][i+1];
				temp3=dp[(j-1+m)%m][i+1];
				temp=temp1;
				if(temp>temp2)
				{
					temp=temp2;
				}
				if(temp>temp3)
				{
					temp=temp3;
				}
				dp[j][i]=map[j][i]+temp;
				if(temp==temp1)
				{
					k=j;
				}
				if(temp==temp2)
				{
					if(k>(j+1)%m)
					{
						k=(j+1)%m;
					}
				}
				if(temp==temp3)
				{
					if(k>(j-1+m)%m)
					{
						k=(j-1+m)%m;
					}
				}
				lujing[j][i]=k;
				
			}
		}
		min=dp[0][0];
		hangshu=0;
		for(i=1;i<m;i++)
		{
			if(min>dp[i][0])
			{
				min=dp[i][0];
				hangshu=i;
			}	
		}
		printf("%d",hangshu+1);
		hangshu=lujing[hangshu][0];
		for(i=1;i<n;i++)
		{
			printf(" %d",hangshu+1);
			hangshu=lujing[hangshu][i];
		}
		printf("\n%d\n",min);
	}
	return 0;
}</span>