bzoj1003 [ZJOI2006]物流运输

Description

---

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

1<=n<=100
1<=m<=20
同一个码头有可能在多个时间段内不可用。但任何时间都存在至少一 条从码头A到码头B的运输路线。

Solution

---

我怎么这么弱啊，不如定个小目标写写dp题好了

用w[i,j]记录第i天到第j天的从起点到终点的最短路，f[i]表示前i天的最小花费，那么就很显然了
其实应该没必要开ll，还有注意特判一下w为inf的情况

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))

typedef long long LL;
const int INF=0x3f3f3f3f;
const int N=2005;
const int E=4005;

struct edge {int x,y,w,next;} e[E];

int queue[N],used[N][N];
int dis[N],vis[N];
int ls[N],edCnt;
LL f[N],w[N][N];

void add_edge(int x,int y,int w) {
    e[++edCnt]=(edge) {x,y,w,ls[x]}; ls[x]=edCnt;
    e[++edCnt]=(edge) {y,x,w,ls[y]}; ls[y]=edCnt;
}

int spfa(int st,int ed,int l,int r) {
    fill(dis,21); dis[st]=0;
    fill(vis,0); vis[st]=1;
    int head=1,tail=1; queue[tail]=st;
    while (head<=tail) {
        int now=queue[head++];
        for (int i=ls[now];i;i=e[i].next) {
            if (!(used[r][e[i].y]-used[l-1][e[i].y])&&(dis[now]+e[i].w<dis[e[i].y])) {
                dis[e[i].y]=dis[now]+e[i].w;
                if (!vis[e[i].y]) {
                    queue[++tail]=e[i].y;
                    vis[e[i].y]=1;
                }
            }
        }
        vis[now]=0;
    }
    return dis[ed];
}

int main(void) {
    int n,m,p,q; scanf("%d%d%d%d",&n,&m,&p,&q);
    rep(i,1,q) {
        int x,y,w; scanf("%d%d%d",&x,&y,&w);
        add_edge(x,y,w);
    }
    int d; scanf("%d",&d);
    rep(i,1,d) {
        int x,l,r; scanf("%d%d%d",&x,&l,&r);
        rep(j,l,r) used[j][x]=1;
    }
    rep(i,1,n) rep(j,1,m) used[i][j]+=used[i-1][j];
    fill(w,31);
    rep(i,1,n) rep(j,i,n) {
        w[i][j]=spfa(1,m,i,j);
    }
    rep(i,1,n) {
        if (w[1][i]==w[0][0]) f[i]=INF;
        else f[i]=w[1][i]*(LL)i;
    }
    rep(i,2,n) {
        rep(j,1,i-1) {
            if (w[j+1][i]!=w[0][0]&&f[j]!=INF)
            f[i]=std:: min(f[i],f[j]+w[j+1][i]*(LL)(i-j)+p);
        }
    }
    printf("%lld\n", f[n]);
    return 0;
}