CH5102 Mobile Service【线性dp】

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 ;
 }