Description

　　物流公司要把一批貨物從碼頭A運到碼頭B。由於貨物量比較大，需要n天才能運完。貨物運輸過程中一般要轉
停好幾個碼頭。物流公司通常會設計一條固定的運輸路線，以便對整個運輸過程實施嚴格的管理和跟蹤。由於各種
因素的存在，有的時候某個碼頭會無法裝卸貨物。這時候就必須修改運輸路線，讓貨物能夠按時到達目的地。但是
修改路線是一件十分麻煩的事情，會帶來額外的成本。因此物流公司希望能夠訂一個n天的運輸計劃，使得總成本

Input

　　第一行是四個整數n（1<=n<=100）、m（1<=m<=20）、K和e。n表示貨物運輸所需天數，m表示碼頭總數，K表示
每次修改運輸路線所需成本。接下來e行每行是一條航線描述，包括了三個整數，依次表示航線連接的兩個碼頭編
號以及航線長度（>0）。其中碼頭A編號為1，碼頭B編號為m。單位長度的運輸費用為1。航線是雙向的。再接下來
一行是一個整數d，後面的d行每行是三個整數P（ 1 < P < m）、a、b（1< = a < = b < = n）。表示編號為P的碼
頭從第a天到第b天無法裝卸貨物（含頭尾）。同一個碼頭有可能在多個時間段內不可用。但任何時間都存在至少一

Output

　　包括了一個整數表示最小的總成本。總成本=n天運輸路線長度之和+K*改變運輸路線的次數。

Sample Input

Sample Output

HINT

題解：

　　這道題目的數據範圍如此小，處理處cost i，j表示i-j這天的最短路

　　然後取dp即可。

 1 #include<iostream>
 2 #include<cstring>
 3 #include<cstdio>
 4 #include<queue>
 5 using namespace std;
 6 
 7 const int max_n=105;
 8 const int max_m=25;
 9 const int max_e=max_m*max_m*2;
10 const int INF=2e9;
11 
12 int n,m,k,e,x,y,z,d,p,a,b;
13 int tot,point[max_m],nxt[max_e],v[max_e],c[max_e];
14 int dis[max_m]; bool broken[max_m][max_n],vis[max_m],flag[max_m];
15 int f[max_n],cost[max_n][max_n];
16 queue <int> q;
17 
18 inline void addedge(int x,int y,int z)
19 {
20     ++tot; nxt[tot]=point[x]; point[x]=tot; v[tot]=y; c[tot]=z;
21     ++tot; nxt[tot]=point[y]; point[y]=tot; v[tot]=x; c[tot]=z;
22 }
23 inline int spfa()
24 {
25     memset(dis,0x7f,sizeof(dis)); dis[1]=0;
26     memset(vis,0,sizeof(vis)); vis[1]=true;
27     while (!q.empty()) q.pop(); q.push(1);
28 
29     while (!q.empty())
30     {
31         int now=q.front(); q.pop();
32         vis[now]=false;
33         for (int i=point[now];i;i=nxt[i])
34             if (dis[v[i]]>dis[now]+c[i]&&!flag[v[i]])
35             {
36                 dis[v[i]]=dis[now]+c[i];
37                 if (!vis[v[i]])
38                 {
39                     vis[v[i]]=true;
40                     q.push(v[i]);
41                 }
42             }
43     }
44     return dis[m];
45 }
46 int main()
47 {
48     scanf("%d%d%d%d",&n,&m,&k,&e);
49     for (int i=1;i<=e;++i)
50     {
51         scanf("%d%d%d",&x,&y,&z);
52         addedge(x,y,z);
53     }
54     scanf("%d",&d);
55     for (int i=1;i<=d;++i)
56     {
57         scanf("%d%d%d",&p,&a,&b);
58         for (int j=a;j<=b;++j) broken[p][j]=true;
59     }
60     for (int i=1;i<=n;++i)
61         for (int j=i;j<=n;++j)
62         {
63             memset(flag,0,sizeof(flag));
64             for (int k=1;k<=m;++k)
65                 for (int l=i;l<=j;++l)
66                     flag[k]|=broken[k][l];
67             cost[i][j]=spfa();
68         }
69     for (int i=1;i<=n;++i)
70         for (int j=i;j<=n;++j)
71             if (cost[i][j]<INF) cost[i][j]*=(j-i+1);
72     memset(f,0x7f,sizeof(f));
73     for (int i=1;i<=n;++i) f[i]=cost[1][i];
74     for (int i=2;i<=n;++i)
75         for (int j=1;j<i;++j)
76             f[i]=min(f[i],f[j]+cost[j+1][i]+k);
77     printf("%d\n",f[n]);
78 }

bzoj 1003 [ZJOI2006]物流運輸（最短路+dp）