Acwing 274.移动服务 (DP)

题面

一个公司有三个移动服务员,最初分别在位置1,2,3处。

如果某个位置(用一个整数表示)有一个请求,那么公司必须指派某名员工赶到那个地方去。

某一时刻只有一个员工能移动,且不允许在同样的位置出现两个员工。

从 p 到 q 移动一个员工,需要花费 c(p,q)。

这个函数不一定对称,但保证 c(p,p)=0。

给出N个请求,请求发生的位置分别为 p1~pN。

公司必须按顺序依次满足所有请求,且过程中不能去其他额外的位置,目标是最小化公司花费,请你帮忙计算这个最小花费。

输入格式
第1行有两个整数L,N,其中L是位置数量,N是请求数量,每个位置从1到L编号。

第2至L+1行每行包含L个非负整数,第i+1行的第j个数表示c(i,j) ,并且它小于2000。

最后一行包含N个整数,是请求列表。

一开始三个服务员分别在位置1,2,3。

输出格式
输出一个整数M,表示最小花费。

数据范围
3≤L≤200,
1≤N≤1000
输入样例:
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

思路

首先我们直接考虑四维数组dp,但看起来会炸,那么优化,我们其实可以压缩成为三维的,因为这三个邮递员里面必定有一个是等于当前的要求的地点,所以我们不用单独开一维去存。那么接下来我们考虑状态转移,这道题目我们发现从后往前推并不是很方便,所以我们考虑用当前状态去更新后面的状态,那么具体的转移就是枚举是哪个邮递员去了下一个地点,并且加上邻接矩阵里面的对应距离,就是下一个状态的值了。要记得不能有多个邮递员在同一地点的判断哦。

代码实现

#include<cstdio>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
#define rep(i,f_start,f_end) for (int i=f_start;i<=f_end;++i)
#define per(i,n,a) for (int i=n;i>=a;i--)
#define MT(x,i) memset(x,i,sizeof(x) )
#define rev(i,start,end) for (int i=0;i<end;i++)
#define inf 0x3f3f3f3f
#define mp(x,y) make_pair(x,y)
#define lowbit(x) (x&-x)
#define MOD 1000000007
#define exp 1e-8
#define N 1000005 
#define fi first 
#define se second
#define pb push_back
typedef long long ll;
typedef pair<int ,int> PII;
ll gcd (ll a,ll b) {return b?gcd (b,a%b):a; }
inline int read() {
    char ch=getchar(); int x=0, f=1;
    while(ch<'0'||ch>'9') {
        if(ch=='-') f = -1;
        ch=getchar();
    } 
    while('0'<=ch&&ch<='9') {
        x=x*10+ch-'0';
        ch=getchar();
    }   return x*f;
}

const int maxn=201;
const int M=1010;
int f[M][maxn][maxn];
int maze[maxn][maxn];
int n,m;
int ask[M];

int main () {
    cin>>n>>m;
    rep (i,1,n)
     rep (j,1,n) {
         cin>>maze[i][j];
    }
    MT (ask,0);
    rep (i,1,m) cin>>ask[i];                               
    ask[0]=3;
    int ans=0x3f3f3f3f;
    MT (f,0x3f);
    f[0][1][2]=0;
    rev (i,0,m) 
     rep (x,1,n)
       rep (y,1,n) {
           int z=ask[i],u=ask[i+1],w=f[i][x][y];
           if (x==y||y==z||x==z) continue;
           f[i+1][x][y]=min (f[i+1][x][y],w+maze[z][u]);
           f[i+1][x][z]=min (f[i+1][x][z],w+maze[y][u]);
           f[i+1][z][y]=min (f[i+1][z][y],w+maze[x][u]);
       } 
    rep (i,1,n) 
     rep (j,1,n) {
         int z=ask[m];
         if (i==j||j==z||z==i) continue;
         ans=min (ans,f[m][i][j]);
     }
    cout<<ans<<endl;
    return 0;
}
原文地址:https://www.cnblogs.com/hhlya/p/13425965.html