Mobile Service_纪中_1327

题目大意

    一个公司有三个移动服务员。如果某个地方有一个请求，某个员工必须赶到那个地方去（那个地方没有其他员工），某一时刻只有一个员工能移动。被请求后，他才能移动，不允许在同样的位置出现两个员工。从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.