【bzoj1003】【最短路+DP】【ZJOI2006】物流运输trans



物流公司要把一批货物从码头A运到码头B。由于货物量比较大，需要n天才能运完。货物运输过程中一般要转停好几个码头。物流公司通常会设计一条固定的运输路线，以便对整个运输过程实施严格的管理和跟踪。由于各种因素的存在，有的时候某个码头会无法装卸货物。这时候就必须修改运输路线，让货物能够按时到达目的地。但是修改路线是一件十分麻烦的事情，会带来额外的成本。因此物流公司希望能够订一个n天的运输计划，使得总成本尽可能地小。

思路：先预处理出对于任意连续的第L天到第R天，在中间不改变运输路线时的最短路径。然后DP即可，用F[i]表示从第一天到第i天的最小成本，考虑最后一次改变路线是哪一天，就不难写出DP方程了

var
  a:array[1..20,1..20]of longint;
  dis:array[1..20]of longint;
  vis:array[1..20]of boolean;
  use:array[1..20,1..100]of boolean;
 
  f:array[1..100,1..100]of longint;
  dp:array[0..100]of longint;
 
  n,m,k,e,i,j:longint;
  conmax:longint;
 
procedure readdata;
var b,c,d,f:longint;
begin
  readln(n,m,k,e);
  for i:=1 to e do
  begin
    readln(b,c,d);
    a[b,c]:=d;
    a[c,b]:=d;
  end;
 
  fillchar(use,sizeof(use),true);
  readln(d);
  for i:=1 to d do
  begin
    readln(b,c,f);
    for j:=c to f do use[b,j]:=false;
  end;
end;
 
procedure qzdl(x,y:longint);
var i,j,k,min,mini:longint;
begin
  fillchar(dis,sizeof(dis),$7f);
  dis[1]:=0;
  fillchar(vis,sizeof(vis),false);
  for i:=1 to m do
  for j:=x to y do
    if use[i,j]=false then vis[i]:=true;
 
  while vis[m]=false do
  begin
    min:=maxlongint;
    for i:=1 to m do
    if (vis[i]=false)and(dis[i]<min) then
    begin
      min:=dis[i];
      mini:=i;
    end;
    i:=mini;
    vis[i]:=true;
    for j:=1 to m do
    if a[i,j]<>0 then
    if dis[i]+a[i,j]<dis[j] then
      dis[j]:=dis[i]+a[i,j];
  end;
 
  if dis[m]=conmax then f[x,y]:=dis[m]
    else f[x,y]:=dis[m]*(y-x+1);
end;
 
begin
  readdata;
  fillchar(conmax,sizeof(conmax),$7f);
 
  for i:=1 to n do
  for j:=i to n do
    qzdl(i,j);
 
  {for i:=1 to n do
  for j:=i to n do
    writeln(i,' ',j,' ',f[i,j]);}
 
  fillchar(dp,sizeof(dp),$7f);
  dp[0]:=0;
  for i:=1 to n do
  for j:=0 to i-1 do
  if dp[j]+f[j+1,i]+k < dp[i] then
    dp[i]:=dp[j]+f[j+1,i]+k;
 
  writeln(dp[n]-k);
end.