[ZJOI2006]BZOJ 1003 物流运输 - 最短路 - dp

121 篇文章 0 订阅
40 篇文章 0 订阅

传送门

这个题想明白了还是还是非常的,就是注意到对T天分为若干段,每段的费用一定选择可行最短路。

然后就是dp,没了。

//ZJOI 2006
//BZOJ 1003
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<climits>
#include<algorithm>
#define INF (LLONG_MAX/3)
#define MAXN 30
#define MAXT 110
#define MAXM 1000
#define lint long long
#define open oppppen
using namespace std;
int t,n,k,m,etop,h[MAXN];
lint g[MAXN][MAXN],f[MAXT];
lint d[MAXN];bool vis[MAXN];
bool open[MAXN][MAXT],nc[MAXN];
struct edges{
    int to,wgt,pre;
}e[MAXM];
inline int add_edge(int u,int v,int w)
{
    etop++;
    e[etop].to=v;
    e[etop].wgt=w;
    e[etop].pre=h[u];
    h[u]=etop;
    return 0;
}
inline int init_edges(int L,int R)
{
    memset(h,0,sizeof(h));etop=0;
    memset(nc,true,sizeof(nc));
    for(int i=1;i<=n;i++)
        for(int j=L;j<=R;j++)
            if(!open[i][j])
            {
                nc[i]=false;
                break;
            }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if((i^j)&&nc[i]&&nc[j]&&g[i][j]<INF)
                add_edge(i,j,g[i][j]);
    return 0;
}
struct node{
    int id;lint dist;
    node(int _i,lint _d)
    {
        id=_i,dist=_d;
    }
    bool operator<(const node &n)const
    {
        return this->dist>n.dist;
    }
};
priority_queue<node> q;
lint dijkstra(int s,int t,int L,int R)
{
    memset(vis,false,sizeof(vis));
    while(!q.empty()) q.pop();
    for(int i=1;i<=n;i++) d[i]=INF;
    q.push(node(s,d[s]=0));
    init_edges(L,R);
    while(!q.empty())
    {
        int x=q.top().id;q.pop();
        if(vis[x]) continue;vis[x]=true;
        for(int i=h[x];i;i=e[i].pre)
            if(!vis[e[i].to]&&d[e[i].to]>d[x]+e[i].wgt)
            {
                d[e[i].to]=d[x]+e[i].wgt;
                q.push(node(e[i].to,d[e[i].to]));
            }
    }
    return d[t];
}
int main()
{
    scanf("%d%d%d%d",&t,&n,&k,&m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if(i^j) g[i][j]=INF;
            else g[i][j]=0;
    for(int i=1;i<=m;i++)
    {
        int u,v,w;scanf("%d%d%d",&u,&v,&w);
        g[u][v]=g[v][u]=min(g[u][v],(lint)w);
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=t;j++)
            open[i][j]=true;
    int T;scanf("%d",&T);
    while(T--)
    {
        int x,u,v;scanf("%d%d%d",&x,&u,&v);
        for(int i=u;i<=v;i++) open[x][i]=false;
    }
    for(int i=1;i<=t;i++) f[i]=INF;
    for(int i=1;i<=t;i++)
        for(int j=0;j<i;j++)
        {
            lint v=dijkstra(1,n,j+1,i);
            if(v>=INF) continue;
            f[i]=min(f[i],f[j]+v*(i-j)+((int)(j>0))*k);
        }
    printf("%lld\n",f[t]);return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值