动态规划Always On the Run poj 1476解题报告（附详细分析）

这道是动态规划，动态方程不难想，只是要实现数据的录入需要一些技巧。

以第i次飞行和第m个城市作为状态dp

动态转移方程：dp[i][m] = min {dp[i][m] ，,dp[i-1][j]+map[j][m][p]}   其中的map[i][k][p]是指第j城市到第k城市的第p个航班的费用,具体看代码和注释。

#include <iostream>
using namespace std;
const int INF = 9999999;
int dp[1001][11], map[11][11][1001];		//三维数组储存每日每次航班的费用，dp状态为第k次飞行，第n个城市
int main()
{
	int n, k, p, counter = 1;
	while (scanf("%d%d", &n, &k))
	{
		if(n + k == 0)
			break;
		memset(map, 0, sizeof(map));
		for(int i = 1; i <= n; i++)
			for(int j = 1; j <= n; j++)
				if(i != j)				//忽略i城市与飞向自己的航班j，即满足条件n(n-1)的航班次数
				{
					scanf("%d", &map[i][j][0]);		//存储当天的航班次数
					for(int m = 1; m <= map[i][j][0]; m++)		
						scanf("%d", &map[i][j][m]);		//存储每次航班的费用
				}
		for(int i = 0; i <= k; i++)		//初始化动态规划，是“装满背包型”问题
			for(int j = 0; j <= n; j++)
				dp[i][j] = INF;
		dp[0][1] = 0;
		for(int i = 1; i <= k; i++)				//动态规划开始
			for(int j = 1; j <= n; j++)
				if(dp[i-1][j] != INF)			
					for(int m = 1; m <= n; m++)
						if(j != m)					//新建一个循环来跳过同一个城市自身的航班比较
						{
							p = i % map[j][m][0];		//得到飞向m城市的航班
							if(p == 0)	p = map[j][m][0];		//一个飞行航班周期
							if(map[j][m][p] == 0)	continue;		//跳过无法到达m城市的航班
							if(dp[i][m] > dp[i-1][j]+map[j][m][p])	//dp[i-1][j]是前一状态的值，map[j][m][p]是指第j城市到第m城市的第p个航班的费用
								dp[i][m] = dp[i-1][j]+map[j][m][p];
						}
		printf("Scenario #%d\n",counter++);
		if(dp[k][n]==INF)				//如果经过所有航班无法到达目的地，那么dp[k][n](经过k次飞行n个城市的费用）的初始值没有改变
			printf("No flight possible.\n\n");
		else
			printf("The best flight costs %d.\n\n",dp[k][n]);
		}
	return 0;
}