题目大意
一个公司有三个移动服务员。如果某个地方有一个请求,某个员工必须赶到那个地方去(那个地方没有其他员工),某一时刻只有一个员工能移动。被请求后,他才能移动,不允许在同样的位置出现两个员工。从p到q移动一个员工,需要花费c(p,q)。这个函数没有必要对称,但是c(p,p)=0。公司必须满足所有的请求。目标是最小化公司花费。
分析
这是一个dp,这一眼就可以看出。
但是暴力的dp是会爆空间和时间的。
所以要改一改。
首先,要用滚动数组,因为我们的空间不足。
具体方程见程序。
还要手动把第一个请求设为初值。
代码
var
f:array[0..1,1..300,1..300] of longint;
c:array[1..300,1..300] of longint;
a:array[0..2000] of longint;
i,j,k,x,y:longint;
n,m:longint;
ans:longint;
procedure init;
begin
readln(n,m);
for i:=1 to n do
begin
for j:=1 to n do
read(c[i,j]);
readln;
end;
for i:=1 to m do
read(a[i]);
end;
function max(x,y:longint):longint;
begin
if x<y then max:=x
else max:=y;
end;
begin
init;
fillchar(f,sizeof(f),$7f);
y:=f[0,1,2];
f[0,1,2]:=c[3,a[1]];
f[0,1,3]:=c[2,a[1]];
f[0,2,3]:=c[1,a[1]];
f[0,2,1]:=f[0,1,2];
f[0,3,1]:=f[0,1,3];
f[0,3,2]:=f[0,2,3];
x:=0;
for i:=2 to m do
begin
x:=x xor 1;
for j:=1 to n do for k:=1 to n do f[x,j,k]:=y;
for j:=1 to n do
for k:=1 to n do
if (j<>k) and (j<>a[i-1]) and (k<>a[i-1])
then
begin
f[x,j,k]:=max(f[x,j,k],f[x xor 1,j,k]+c[a[i-1],a[i]]);
f[x,k,j]:=f[x,j,k];
f[x,a[i-1],k]:=max(f[x,a[i-1],k],f[x xor 1,j,k]+c[j,a[i]]);
f[x,k,a[i-1]]:=f[x,a[i-1],k];
f[x,a[i-1],j]:=max(f[x,a[i-1],j],f[x xor 1,j,k]+c[k,a[i]]);
f[x,j,a[i-1]]:=f[x,a[i-1],j];
end;
end;
ans:=maxlongint;
for i:=1 to n do
for j:=1 to n do
if ans>f[x,i,j] then ans:=f[x,i,j];
write(ans);
end.