5102 Mobile Service 0x50「动态规划」例题
描述
一个公司有三个移动服务员，最初分别在位置1，2，3处。
如果某个位置（用一个整数表示）有一个请求，那么公司必须指派某名员工赶到那个地方去。某一时刻只有一个员工能移动，且不允许在同样的位置出现两个员工。从 p 到 q 移动一个员工，需要花费 c(p,q)。这个函数不一定对称，但保证 c(p,p)=0。
给出N个请求，请求发生的位置分别为 p_1~p_N。公司必须按顺序依次满足所有请求，目标是最小化公司花费，请你帮忙计算这个最小花费。N≤1000，位置是1~200的整数。
输入格式
第一行有两个整数L,N(3<=L<=200, 1<=N<=1000)。L是位置数；N是请求数。每个位置从1到L编号。下L行每行包含L个非负整数。第i+1行的第j个数表示c(i,j) ，并且它小于2000。最后一行包含N个数，是请求列表。一开始三个服务员分别在位置1，2，3。
输出格式
一个数M，表示最小服务花费。
样例输入
5 9
0 1 1 1 1
1 0 2 3 2
1 1 0 4 1
2 1 5 0 1
4 2 3 4 0
4 2 4 1 5 4 3 2 1
样例输出
5
题意：
有n个点标号为1-n，现在有三个人站在1,2,3处
给出一个矩阵 表示从i到j需要的花费
给出m个询问，每次给一个pos表示这一次pos要有人
问 m次询问总的cost 是多少
思路：
dp[i, x, y, z]表示对于第i次查询，服务员分别在xyz时的总花费
那么第i+1次时，只有三种可能。
并且可以发现，其实只需要三维即可，因为三个服务员中一定有一个位置是确定的
他在i+1时一定是在pi，并且服务员之间交换位置是不影响答案的
所以i+1时只需要维护dp[i+1, x, y], dp[i+1, pi, y], dp[i+1, x, pi]
dp[i+1, x, y] = min(dp[i+1, x, y], dp[i, x, y] + c(pi, pi+1), z从pi到pi+1
dp[i+1, pi, y] = min(dp[i+1,pi, y], dp[i, x, y] + c(x, pi+1)，x到pi+1，z与x交换位置
dp[i+1, x, pi] = min(dp[i+1,y, pi], dp[i, x, y] + c(y, pi+1)，y到pi+1， z与y交换位置
需要判断一下 xyz中是否有与pi+1相同的点，若有则只有一种走法了
虐狗宝典笔记：
1.求解线性DP问题，一般先确定“阶段”。若“阶段”不足以表示一个状态，则可以把所需的附加信息也作为状态的维度。
转移时，若总是从一个阶段转移到下一个阶段，则没有必要关心附加信息维度的大小变化情况，因为“无后效性”已经由“阶段”保证
2.在确定dp状态时，要选择最小的能够覆盖整个状态空间的“维度集合”
若DP状态由多个维度构成，则应检查这些维度之间能否相互导出，用尽量少的维度覆盖整个状态空间，排除冗余维度。
//#include <bits/stdc++.h>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<stdio.h>
#include<cstring>
#include<map> #define inf 0x3f3f3f3f
using namespace std;
typedef long long int LL; const int maxn = ;
int c[maxn][maxn], dp[][maxn][maxn], p[];
int l, n; int main()
{
scanf("%d%d", &l, &n);
for(int i = ; i <= l; i++){
for(int j = ; j <= l; j++){
scanf("%d", &c[i][j]);
}
}
for(int i = ; i <= n; i++){
scanf("%d", &p[i]);
}
memset(dp, 0x3f, sizeof(dp)); dp[][][] = ;
p[] = ;
for(int i = ; i <= n; i++){
for(int x = ; x <= l; x++){
for(int y = ; y <= l; y++){
if(dp[i - ][x][y] != inf){
int z = p[i - ];
if(y != p[i] && x != p[i]){
dp[i][x][y] = min(dp[i - ][x][y] + c[z][p[i]], dp[i][x][y]);
}
if(y != p[i] && z != p[i]){
dp[i][y][z] = min(dp[i - ][x][y] + c[x][p[i]], dp[i][y][z]);
}
if(x != p[i] && z != p[i]){
dp[i][x][z] = min(dp[i - ][x][y] + c[y][p[i]], dp[i][x][z]);
}
dp[i - ][x][y] = 0x3f3f3f3f;
}
}
}
} int ans = inf;
for(int x = ; x <= l; x++){
for(int y = ; y <= l; y++){
ans = min(ans, dp[n][x][y]);
}
}
printf("%d\n", ans);
return ;
}
