bzoj 1003 【zjoi2006】物流运输

好吧。我坦白。因为题意有点很多地方不了解查了tijie
Description

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

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

分析：
所以看完题解后从新分析。。。

一开始没有想明白的是改变运输路线的次数是什么玩意。。总感觉有可能经过最短路的时候并没有经过不可以经过的点所以好麻烦。。就比如样例三天可以经过一条路线。。这怎么办？
但是有一点很清楚。就是要求最短路。。。

所以几天走同一条路径既ok又花费最少呢。。
这个很像dp模型啊。类似由上一个有用的状态转移到这个状态。
所以先把dp方程搞出来。
f[i]表示1-第i所用的成本。

f[i]=min(f[i],f[j]+temp+K); 0<=j<=i-1;

这个tmp代表着从第j+1天-第i天中间可以走的话的代价。这个用spfa搞定就ok 结果dis[m]*（i-j+1）。不然就是无穷大。

又有问题出现了。怎么判断能不能走呢？
利用前缀和就好。记录path[i][j]为码头i在第j天是否可以走，可以走为0，不能走为1 ， path[i][j]+=path[i][j-1]
如果在dp时的path[v][i]-path[v][j-1]!=0就是中间有天不能走。

其实感觉这道题的精华在于dp方程。
其实很多题是分层剖析可以套出来2+个模型的。
看代码喽~

#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;
//by mars_ch
int n,m,K,E,d;
struct data
{
    int u,v,nxt,w;
}e[10005*2];
int first[1005],tot;
int dis[1005],f[1005],path[25][105],inq[105];
void add(int a,int b,int c)
{
    e[tot].u=a;
    e[tot].v=b;
    e[tot].w=c;
    e[tot].nxt=first[a];
    first[a]=tot++;
}
int spfa(int a,int b)
{
    queue<int> q;
    memset(dis,0x3f,sizeof(dis));
    memset(inq,0,sizeof(inq));
    q.push(1);
    dis[1]=0;
    inq[1]=1;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        inq[u]=0;

        for(int i=first[u];i!=-1;i=e[i].nxt)
        {
            int v=e[i].v;
            if(path[v][b]-path[v][a-1]!=0) continue;
            if(dis[v]>dis[u]+e[i].w)
            {
                dis[v]=dis[u]+e[i].w;
                if(!inq[v])
                {
                    inq[v]=1;
                    q.push(v);
                }
            }
        }
    }
    return dis[m]>=1061109567?dis[m]:dis[m]*(b-a+1);   //乘错。。。。。写的是（b-a) 
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&K,&E);
    memset(first,-1,sizeof(first));
    for(int i=1;i<=E;i++)
    {
        int a,b,c;
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
        add(b,a,c);
    }
    for(int i=1;i<=m;i++)
    {
        for(int j=1;j<=n;j++)
        {
            path[i][j]=0;
        }
    }
    scanf("%d",&d);
    for(int i=1;i<=d;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        for(int j=y;j<=z;j++)
        {
            path[x][j]=1;
        } 
    } 
    for(int i=1;i<=m;i++)
    {
        for(int j=1;j<=n;j++)
        {
            path[i][j]+=path[i][j-1];
        }
    }
    memset(f,0x3f,sizeof(f));
    f[0]=0;
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<=i-1;j++)
        {
            int tmp=spfa(j+1,i);
            if(tmp>=1000000000) continue;
            f[i]=min(f[i],f[j]+tmp+K);
        }
    }
    printf("%d\n",f[n]-K);
    return 0;
}