1. 程式人生 > >BZOJ_1003_[ZJOI2006]物流運輸_最短路+dp

BZOJ_1003_[ZJOI2006]物流運輸_最短路+dp

tar 長度 scanf bsp oid wap ems blank tin

BZOJ_1003_[ZJOI2006]物流運輸_最短路+dp

題意:http://www.lydsy.com/JudgeOnline/problem.php?id=1003

分析:

這種一段一段的顯然要用dp求。

f[i]表示到第i天為止的最小花銷。轉移有f[i]=min{f[j-1]+cost[j][i]*(i-j+1)+k};

其中cost[i][j]表示從i到j天的最短路長度,spfa預處理出來。

代碼:

 1 #include <stdio.h>
 2 #include <string.h>
 3 #include <algorithm>
 4 using namespace
std; 5 int head[30],to[1000],nxt[1000],val[1000],cnt; 6 int n,m,k,t,inb[110][30],f[110],cost[110][110],can[30]; 7 int Q[100],l,r,dis[30],inq[30]; 8 inline void add(int u,int v,int w){ 9 to[++cnt]=v;nxt[cnt]=head[u];head[u]=cnt;val[cnt]=w; 10 } 11 int spfa(int be,int en){ 12 for(int i=1;i<=n;i++)can[i]=0
; 13 for(int i=be;i<=en;i++) 14 for(int j=1;j<=n;j++)if(inb[i][j]) 15 can[j]=1; 16 l=r=0; 17 memset(dis,0x3f,sizeof(dis)); 18 dis[1]=0;Q[r++]=1;inq[1]=1; 19 while(l^r){ 20 int x=Q[l++];if(l==n+10)l=0;inq[x]=0; 21 for(int i=head[x];i;i=nxt[i]){
22 if(can[to[i]])continue; 23 if(dis[to[i]]>dis[x]+val[i]){ 24 dis[to[i]]=dis[x]+val[i]; 25 if(!inq[to[i]]){ 26 inq[to[i]]=1;Q[r++]=to[i];if(r==n+10)r=0; 27 } 28 } 29 } 30 } 31 return dis[n]>100000?-1:dis[n]; 32 } 33 int main(){ 34 scanf("%d%d%d%d",&t,&n,&k,&m); 35 int x,y,z; 36 for(int i=1;i<=m;i++){ 37 scanf("%d%d%d",&x,&y,&z); 38 add(x,y,z);add(y,x,z); 39 } 40 int h; 41 scanf("%d",&h); 42 for(int i=1;i<=h;i++){ 43 scanf("%d%d%d",&x,&y,&z); 44 if(y>z)swap(y,z); 45 for(int j=y;j<=z;j++)inb[j][x]=1; 46 } 47 for(int i=1;i<=t;i++){ 48 for(int j=i;j<=t;j++){ 49 cost[i][j]=spfa(i,j); 50 //printf("%d\n",cost[i][j]); 51 } 52 } 53 f[0]=0; 54 for(int i=1;i<=t;i++){ 55 f[i]=1<<30; 56 for(int j=1;j<=i;j++){ 57 if(cost[j][i]==-1)continue; 58 f[i]=min(f[j-1]+cost[j][i]*(i-j+1)+k,f[i]); 59 } 60 //printf("%d\n",f[i]); 61 } 62 printf("%d",f[t]-k); 63 }

BZOJ_1003_[ZJOI2006]物流運輸_最短路+dp