BZOJ P1003 物流运输 （最短路 + 动态规划）

物流运输是我这个蒟蒻做的第一道BZOJ的题…不对！还有A+B problem！！！哇我在BZOJ做了两道题了。

题目链接：物流运输
因为这道题要求最小代价显然是一道最短路，但是因为有的时候需要该路线增加额外代价，所以还要做一个DP，来求最优策略。
我们来构想一下，当我们到第n天的时候，有两种选择，一是继续原先的路，二是换一条新的路，当然，如果在这一天原先的路不能走了就必须要换一条路，然而换成什么路当然是用最短路求啦。
我们设dp[i]为第i天的最优策略。每天用spfa来求一遍最短路，用来状态转移。状态转移方程为：dp[i]=min(dp[i],dp[j-1]+val*(i-j+1)+k) .val为当天最短路，k为换方案所需的代价。

#include<iostream>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<queue>

using namespace std;

inline int read( )
{
	int x = 0; char c = getchar( );
	while(!isdigit(c)) c = getchar( );
	while(isdigit(c)) x = (x<<1) + (x<<3) + (c^48), c = getchar( );
	return x;
}

struct edge{
	int next,to,val;
}z[1000];

int n,m,k,e,head[1000],cnt,dp[200],d,a[21][110],dis[21],vis[21],now[21];
queue<int> q;

void add(int x,int y,int v)
{
	z[++cnt].to = y;
	z[cnt].val = v;
	z[cnt].next = head[x];
	head[x] = cnt;
}

int spfa()
{
	memset(dis,0x3f3f3f3f,sizeof(dis));
	memset(vis,0,sizeof(vis));
	q.push(1);dis[1] = 0;
//	vis[1] = 1;
	while(!q.empty())
	{
		int u = q.front( ); q.pop( );
//		if(vis[u])continue;
		vis[u] = 0;
		for(int j = head[u];j; j = z[j].next)
		{
			if(now[z[j].to]) continue;
			if(dis[z[j].to] > dis[u] + z[j].val)
			{
//				cout<<"@#@";
				dis[z[j].to] = dis[u] + z[j].val;
				if(!vis[z[j].to])
				{
					q.push(z[j].to);
					vis[z[j].to] = 1;
				}
			}
		}
//		vis[u] = 0;
	}
	return dis[m];
}

int main( )
{
//	scanf("%d%d%d%d",&n,&m,&k,&e);
	n = read( );
	m = read( );
	k = read( );
	e = read( );
	for(int i = 1; i <= e; i++)
	{
		int a,b,c;
		a = read( );
		b = read( );
		c = read( );
//		scanf("%d%d%d",&a,&b,&c);
		add(a,b,c);
		add(b,a,c);
	}
	d = read( );
	for(int i = 1; i <= d; i++)
	{
		int qwq;
		qwq = read( );
		int o , k;
		o = read( );
		k = read( );
		for(int f = o; f <= k; f++)
		{
			a[qwq][f] = 1;
		}
	}
	memset(dp, 0x7f7f7f7f, sizeof(dp));
	dp[0] = -k;
	for(int i = 1; i <= n; i++)
	{
		memset(now,0,sizeof(now));
		for(int j = i; j >= 1; j--)
		{		
			for(int l = 1; l <= m; l++)
			{
				now[l]|=a[l][j];
			}
		int vall = spfa( );
		if(dis[m] != 0x3f3f3f3f)
		{
			dp[i]=min(dp[i],dp[j-1]+vall*(i-j+1)+k);
		}
	}
	}
	printf("%d",dp[n]);
	return 0;
}