Mobile Service TYVJ1061（线性DP）

题意：一个公司有三个移动服务员，开始三个服务员分别在位置1，2，3。
如果某个位置（用一个整数表示）有一个请求，某个员工必须赶到那个地方去（那个地方没有其他员工），某一时刻只有一个员工能移动。被请求后，他才能移动，不允许在同样的位置出现两个员工。从p到q移动一个员工，需要花费c(p,q)，c(p,p)=0。公司必须满足所有的请求。目标是最小化公司花费。

输入格式：
第一行有两个整数L,N(3<=L<=200, 1<=N<=1000)。L是位置数；N是请求数。
每个位置从1到L编号。下L行每行包含L个非负整数。第i+1行的第j个数表示c(i,j) ，并且它小于2000。最后一行包含N个数，是请求列表。

输出格式：
一个数M，表示最小服务花费。

输入样例 #1
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

输出样例 #1
5

思路：
使用动态规划解决本题，容易发现DP的“阶段”就是“已完成的请求数量”，通过指派一个服务员，可以把一个“完成i-1个请求”状态转移到“完成i个请求”的状态。

为了计算指派服务员的花费，就必须知道状态转移时每个服务员的位置，最直接的想法就是把三个服务的位置放在DP的“状态”中。

设F[i，x，y，z]表示完成i个请求，三个服务员位于x，y，z时，公司花费最小。

考虑一下F[i，x，y，z]能够更新i+1阶段的那些状态。就是派三个服务员之一去第i+1个请求处。

F[i+1,Pi+1,y,z]=min(F[i+1,Pi+1,y,z],F[i，x，y，z]+c(x,Pi+1))
F[i+1,x,Pi+1,z]=min(F[i+1,x,Pi+1,z],F[i，x，y，z]+c(y,Pi+1))
F[i+1,x,y,Pi+1]=min(F[i+1,x,y,Pi+1],F[i，x，y，z]+c(z,Pi+1))
该算法的规模大约在1000*200³，不能承受。仔细观察发现，在第i个请求完成是，一定有某个员工位于位置Pi，只需知道阶段i和另两个员工的位置即可描述一个状态。
因此，可用F[i，x，y]表示完成了前i个请求，其中一个员工位于Pi，另两个员工位于x和y时，公司当前的最小花费。

三种转移则分别表示让位于Pi，x和y的员工前往Pi+1处理请求
F[i+1,x,y]=min(F[i+1,x,y],F[i,x,y]+c(Pi，Pi+1))
F[i+1,Pi,y]=min(F[i+1,Pi,y],F[i,x,y]+c(x,Pi+1))
F[i+1,x,Pi]=min(F[i+1,x,Pi],F[i,x,y]+c(y,Pi+1))
不妨设P0=3，于是初值可设为F[0,1,2]=0,目标F[N，？，？]。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int f[2][210][210],p[1010],c[210][210];
int n,m;
int main()
{
	cin>>m>>n;
	for(int i=1;i<=m;i++)
		for(int j=1;j<=m;j++)
			cin>>c[i][j];
	for(int i=1;i<=n;i++) cin>>p[i];
	p[0]=3;
	memset(f,0x3f,sizeof(f));
	f[0][1][2]=0;
	for(int i=1;i<=n;i++)
		for(int x=1;x<=m;x++)
			for(int y=1;y<=m;y++)
				if(f[i-1&1][x][y]!=0x3f3f3f3f)
				{
					int z=p[i-1];
					// the problem doesn't allow two people at the same place
					if(y!=p[i]&&z!=p[i]) f[i&1][y][z]=min(f[i&1][y][z],f[i-1&1][x][y]+c[x][p[i]]);
					if(x!=p[i]&&z!=p[i]) f[i&1][x][z]=min(f[i&1][x][z],f[i-1&1][x][y]+c[y][p[i]]);
					if(x!=p[i]&&y!=p[i]) f[i&1][x][y]=min(f[i&1][x][y],f[i-1&1][x][y]+c[z][p[i]]);
					f[i-1&1][x][y]=0x3f3f3f3f;
				}
	int ans=1<<30;
	for(int x=1;x<=m;x++)
		for(int y=1;y<=m;y++)
			ans=min(ans,f[n&1][x][y]);
	cout<<ans<<endl;
}

这道题给我们两点启发：
1、求解线性DP问题，一般先确定“阶段”。若“阶段”不足以表示一个状态，则可以把所需的附加信息也作为状态的维度。
在转移时，若总是从一个阶段转移到下一个阶段（本题i到i+1）,则没有必要关心附加信息维度大小变化情况（本题x，y，z在转移前后大小不定），因为无后效性已经由“阶段”保证。
2、在确定DP状态是，要选择最小的能够覆盖整个状态空间的“维度集合”。