【ZJOI2006】物流运输trans

Description

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

Input

第一行是四个整数N（1<=N<=100）、M（1<=M<=20）、K和E。 

    N表示货物运输所需天数

    M表示码头总数

    K表示每次修改运输路线所需成本

接下来E行每行是一条航线描述，包括了三个整数，依次表示航线连接的两个码头编号以及航线长度（>0）。其中码头A编号为1，码头B编号为m。单位长度的运输费用为1。航线是双向的。 
再接下来一行是一个整数D，后面的D行每行是三个整数P（ 1 < P < m）、a、b（1 < = a < = b < = n）。 
表示编号为P的码头从第a天到第b天无法装卸货物（含头尾）。同一个码头有可能在多个时间段内不可用。但任何时间都存在至少一条从码头A到码头B的运输路线。

Output

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

Sample Input

5 5 10 8
1 2 1
1 3 3
1 4 2
2 3 2
2 4 4
3 4 1
3 5 2
4 5 2
4
2 2 3
3 1 1             
3 3 3
4 4 5

Sample Output

32

Hint

前三天走1-4-5，后两天走1-3-5，这样总成本为(2+2)*3+(3+2)*2+10=32 

【分析】

DP+SPFA

        仔细分析后我们知道，某些天数的路线可以是一样的，而且一个路线可以持续很多天。这就提示我们可以DP，我们以天数为阶段划分状态。

定义状态数组：f[i]表示1～i天的最小费用。

        状态转移方程 f[i]=min(f[j]+Cost[j+1][i]*(i-j)+K).  （0<=j<i）

        Cost[i][j]表示从第i天到第j天的最小费用（含端点）。表示从j+1天到i天用同一种路线。

        边界f[0]=-K.

注意到数据范围很小，于是我们可以暴力SPFA求出Cost数组，本题便完美解决了。

【代码】

/*
    ID:Ciocio
	LANG:C++
	DATE:2013-11-29
	TASK:Trans
*/
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <queue>

using namespace std;

#define MAXN 105
#define INF 999999999

int N,M,K,E,D;
int G[MAXN][MAXN],Cost[MAXN][MAXN];
bool Fbd[MAXN][MAXN],cntgo[MAXN],vis[MAXN];
int dis[MAXN],f[MAXN];

void _init()
{
	scanf("%d%d%d%d",&N,&M,&K,&E);
	for(int i=0;i<=M;i++)
	{
		for(int j=0;j<=M;j++)
			G[i][j]=INF;
		G[i][i]=0;
	}
	for(int i=1;i<=E;i++)
	{
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		if(G[x][y]>z) G[x][y]=G[y][x]=z;
	}
	scanf("%d",&D);
	for(int i=1;i<=D;i++)
	{
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		for(;y<=z;y++) Fbd[x][y]=true;
	}
}

queue <int> Q;
int _SPFA(int begin,int end)
{
	while(!Q.empty()) Q.pop();
	for(int i=1;i<=M;i++) {dis[i]=INF;vis[i]=false;cntgo[i]=false;}
	for(int i=1;i<=M;i++)
		for(int j=begin;j<=end;j++)
			if(Fbd[i][j])
			{
				cntgo[i]=true;
				break;
			}
	Q.push(1);vis[1]=true;dis[1]=0;
	while(!Q.empty())
	{
		int t=Q.front();Q.pop();
		vis[t]=false;
		for(int i=1;i<=M;i++)
			if((!cntgo[i])&&dis[i]>dis[t]+G[t][i])
			{
				dis[i]=dis[t]+G[t][i];
				if(!vis[i])
				{
					vis[i]=true;
					Q.push(i);
				}
			}
	}
	return dis[M];
}

void _solve()
{
	for(int i=1;i<=N;i++)
		for(int j=i;j<=N;j++)
			Cost[i][j]=_SPFA(i,j);
	for(int i=1;i<=N;i++) f[i]=INF;
	f[0]=-K;
	for(int i=1;i<=N;i++)
		for(int j=0;j<i;j++)
			if(Cost[j+1][i]!=INF)
			    f[i]=min(f[i],f[j]+Cost[j+1][i]*(i-j)+K);
	cout<<f[N]<<endl;
}

int main()
{
	_init();
	_solve();
	return 0;
}