BZOJ 1003 dp+最短路

1003: [ZJOI2006]物流运输

题意:m个码头,从码头1到码头m,连续n天都要运送货物。每一天的花费是总路线长度大小,但如果和前一天的路线不一样,要另处加上k元花费。而且有些码头有些天不能用,问这n天的最小费用。

tags:菜鸡一开始真没想到是dp

求n天时最小花费,就要想到以天数为阶段进行规划。dp[i][j]表示第i天到第j天走同一条路线的花费,则f[i]=min( f[i], f[j]+dp[j+1][i]+k )。

#include<bits/stdc++.h>
using namespace std;
#pragma comment(linker, "/STACK:102400000,102400000")
#define FF(i,a,b) for (int i=a;i<=b;i++)
#define F(i,b,a)  for (int i=b;i>=a;i--)
#define mes(a,b)  memset(a,b,sizeof(a))
#define INF 0x3f3f3f3f
typedef long long ll;
const int N = 200;
const ll inf=1e18;

int n, m, k, e1, d1, q[10000], head[N], tot;
ll c[N][N], dis[N], f[N];
bool vis[N], use[N], flag[N][N];
struct Edge{int to, next; ll w;}e[N*N];
void Addedge(int u, int v, ll w) { e[++tot].to=v, e[tot].w=w, e[tot].next=head[u], head[u]=tot; };
ll spfa(int x, int y)
{
    mes(use, 0);  mes(vis, 0);
    FF(i,1,m) dis[i]=inf;   dis[1]=0;
    FF(i,x,y) FF(j,1,m) if(flag[j][i]) use[j]=1;
    int t=0, h=0;
    q[++t]=1, dis[1]=0;
    while(h!=t) {
        int u=q[++h]; vis[u]=1;
        for(int i=head[u]; i; i=e[i].next) {
            int v=e[i].to;  ll w=e[i].w;
            if(use[v]==0 && dis[v]>dis[u]+w) {
                dis[v]=dis[u]+w;
                if(vis[v]==0) q[++t]=v, vis[v]=1;
            }
        }
        vis[u]=0;
    }
    return dis[m]*(dis[m]>=inf ? 1 : (y-x+1));
}
void Init()
{
    scanf("%d %d %d %d", &n, &m, &k, &e1);
    int u, v, p, a, b;
    ll w;
    FF(i,1,e1) {
        scanf("%d %d %lld", &u, &v, &w);
        Addedge(u, v, w);  Addedge(v, u, w);
    }
    scanf("%d", &d1);
    FF(i,1,d1) {
        scanf("%d %d %d", &p, &a, &b);
        FF(j,a,b) flag[p][j]=1;
    }
    FF(i,1,n) FF(j,1,i) c[j][i]=spfa(j, i);
}
void DP()
{
    FF(i,1,n) {
        f[i]=c[1][i];
        FF(j,1,i-1) f[i]=min(f[i], f[j]+c[j+1][i]+k);
    }
    printf("%lld\n", f[n]);
}
int main()
{
    //freopen("in.txt", "r", stdin);
    Init();
    DP();

    return 0;
}
posted @ 2017-02-18 02:03  v9fly  阅读(191)  评论(0编辑  收藏  举报