BZOJ1003

[BZOJ1003]ZJOI2006物流运输
神奇的dp。。。

目录

思路


  • 首先,若我们知道任意一段时间 [i,j] 的最小花费 cost[i][j] ,很容易得出dp方程:

dp[i]=mini1j=1{dp[j]+cost[j+1][i]+k}

表示到第i天时的最小花费(最开始 dp[i]=cost[1][i] ,就是不该航线的情况)
  • 那么我们就暴力跑最短路好了。。。求出 cost[i][j] 就好了,其实挺水的。

代码

#include <cstdio>
#include <queue>
#include <utility>
#include <vector>
#include <functional>
#include <algorithm>

using std::vector;
using std::greater;
using std::priority_queue;
using std::pair;
using std::min;

const int MAXN=105,MAXM=25,INF=1e9;
int n,m,k,e;
int cost[MAXN][MAXN],dp[MAXN];
bool cant[MAXM][MAXN];

priority_queue<pair<int,int>,vector<pair<int,int> >,greater<pair<int,int> > > hp;

struct node{
    int dis,he;
    bool use,tag;
}d[MAXM];
struct line{
    int to,nex,w;
}ed[MAXM*MAXM<<1];

inline void addE(int u,int v,int w){
    static int cnt=0;
    ed[++cnt]=(line){v,d[u].he,w};
    d[u].he=cnt;
}

inline int Dijsktra(){
    for(int i=1;i<=m;++i){
        d[i].dis=INF;
        d[i].use=false;
    }
    d[1].dis=0;
    hp.push(pair<int,int>(d[1].dis,1));
    pair<int,int> tmp;
    int &u=tmp.second;
    while(hp.size()){
        tmp=hp.top();
        hp.pop();
        if(d[u].use) continue;
        d[u].use=true;
        for(int i=d[u].he,v,w;i;i=ed[i].nex){
            v=ed[i].to,w=ed[i].w;
            if(d[v].dis>d[u].dis+w && !d[v].tag){
                d[v].dis=d[u].dis+w;
                hp.push(pair<int,int>(d[v].dis,v));
            }
        }
    }
    return d[m].dis;
}

inline void deal(){
    int tmpd;
    scanf("%d",&tmpd);
    for(int i=1,p,a,b;i<=tmpd;++i){
        scanf("%d%d%d",&p,&a,&b);
        for(int j=a;j<=b;++j)
            cant[p][j]=true;
    }
    for(int i=1;i<=n;++i){
        for(int j=i;j<=n;++j){
            for(int k=1;k<=m;++k){
                d[k].tag=false;
                for(int l=i;l<=j;++l)
                    d[k].tag|=cant[k][l];
            }
            cost[i][j]=Dijsktra();
            if(cost[i][j]<INF)
                cost[i][j]*=(j-i+1);
        }
    }
}

inline void DP(){
    for(int i=1;i<=n;++i)
        dp[i]=cost[1][i];
    for(int i=2;i<=n;++i){
        for(int j=1;j<i;++j)
            dp[i]=min(dp[i],dp[j]+cost[j+1][i]+k);
    }
}

int main(){
    scanf("%d%d%d%d",&n,&m,&k,&e);
    for(int i=1,u,v,w;i<=e;++i){
        scanf("%d%d%d",&u,&v,&w);
        addE(u,v,w),addE(v,u,w);
    }
    deal();
    DP();
    printf("%d",dp[n]);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值