洛谷 P1772 [ZJOI2006]物流运输 spfa+dp

思路

最短路+dp（比较暴力吧）

题目描述

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

输入输出格式

输入格式：

第一行是四个整数n(l≤n≤100)、m(l≤m≤20)、K和e。n表示货物运输所需天数，m表示码头总数，K表示每次修改运输路线所需成本，e表示航线条数。接下来e行每行是一条航线描述，包括了三个整数，依次表示航线连接的两个码头编号以及航线长度(>0)。其中码头A编号为1，码头B编号为m。单位长度的运输费用为1。航线是双向的。再接下来一行是一个整数d，后面的d行每行是三个整数P(1

输出格式：

包括了一个整数表示最小的总成本。总成本=n天运输路线长度之和+K*改变运输路线的次数

题解

spfa预处理costi,jcosti,j表示第i天到第j天都跑这条路的最短路。也就是说保证路径上的所有点在i到j天都可以跑的前提下的最短路。
然后dp。fifi表示前i天的最小费用，那么

$$f_j=min(f_i+cost_{i+1,j}+k , f_j)$$

初始值

$$f_i=cost_{1,i}$$

值得注意

每一个车站的关闭并不是单日，对于dp算法时间区间【i,j】内若需要用同一路线，那么只要第x天($x\in[i,j]$)不可通行，那么整个时间段就不可以，所以对于spfa中判断点的可行性，需要给每个站一个持久化的禁止通行标记，在每一次 i 循环 初始化这些标记，利用无后效性，在内层spfa之前加一次车站 i （$i\in[1,n]$）的“ | ”运算即可。

代码

#include<cstdio>
#include<iostream>
#include<queue>
#include<cstring>
using namespace std;
#define ll long long
#define mem(a,b) memset(a,b,sizeof(a))
int day,n,k,m;
struct edge
{
    int to,co,next;
}arr[500];
int head[50],cnt;
void adds(int u,int v,int w)
{
    arr[++cnt].to=v;
    arr[cnt].co=w;
    arr[cnt].next=head[u];
    head[u]=cnt;
}
bool banl[50],ban[50][200],vis[50];

int dis[50],cost[200][200],dp[200];

void spfa()
{
    queue<int>q;
    mem(vis,0);mem(dis,0x3f);
    vis[1]=1;dis[1]=0;
    q.push(1);
    while(!q.empty())
    {
        int u=q.front();q.pop();vis[u]=0;
        for(int e=head[u];~e;e=arr[e].next)
        {
            int v=arr[e].to;
            if(banl[v]==1)continue;
            if(dis[v]>dis[u]+arr[e].co)
            {
                dis[v]=dis[u]+arr[e].co;
                if(vis[v]==0)
                {
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
}

int main()
{
    scanf("%d%d%d%d",&day,&n,&k,&m);
    mem(head,-1);mem(ban,0);
    for(int i=1,u,v,w;i<=m;i++)
    {
        scanf("%d%d%d",&u,&v,&w);
        adds(u,v,w);adds(v,u,w);
    }
    int qday;scanf("%d",&qday);
    for(int i=1,id,l,r;i<=qday;i++)
    {
        scanf("%d%d%d",&id,&l,&r);
        for(int j=l;j<=r;j++)ban[id][j]=1;
    }

    mem(cost,0x3f);
    for(int i=1;i<=day;i++)
    {
        mem(banl,0);
        for(int j=i;j<=day;j++)
        {
            for(int i=1;i<=n;i++)banl[i]|=ban[i][j];
            spfa();cost[i][j]=dis[n];
        }
    }

    for(int i=1;i<=day;i++)
    for(int j=i;j<=day;j++)
    if(cost[i][j]<0x3f3f3f3f)cost[i][j]*=(j-i+1);

    mem(dp,0x3f);
    for(int i=1;i<=day;i++)dp[i]=cost[1][i];
    for(int j=2;j<=day;j++)
    for(int i=1;i<j;i++)
    dp[j]=min(dp[j],dp[i]+cost[i+1][j]+k);
    printf("%d",dp[day]);
//system("pause");
return 0;
}