[ZJOI2006] bzoj 1003 物流运输 最短路 动态规划

转移方程很好想 f[i]=f[j]+cost[j+1][i]*(i-j)+k,其中f[i]表示前i天的最小花费，cost[i][j]表示从第i天到第j天可行的最短路，可预处理出来。

注意我们在选择第一条路线时，本不用加k，但我们的方程是看做从0转化了一次路线，已经加上了k，所以最后答案应为f[n]-k。

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<queue>
#define inf 10000000
using namespace std;
int n,m,x[100000],y[100000],z[100000],k,e,d,jishu=0;
int xuan[500],ying[500],num[500],f[200],cost[200][200],head[50],dis[50];
int xx,yy,zz;
bool flag[500];
queue<int> q;
struct bian
{
  int qi,zhong,w,next;
};
bian c[200000];
bool check(int i,int k,int j)
{
  if(ying[j]<i)  return 1;
  if(xuan[j]>k)  return 1;
  return 0;
}
void add(int x,int y,int z)
{
  c[++jishu].qi=x;
  c[jishu].zhong=y;
  c[jishu].w=z;
  c[jishu].next=head[x];
  head[x]=jishu;
}
void spfa(int x)
{
  dis[x]=0;
  flag[x]=1;
  q.push(x);
  while(!q.empty())
  {
    int t=q.front();
    q.pop();
    for(int i=head[t];i;i=c[i].next)
      if(dis[c[i].zhong]>dis[t]+c[i].w)
      {
        dis[c[i].zhong]=dis[t]+c[i].w;
        if(!flag[c[i].zhong])
        {
          flag[c[i].zhong]=1;
          q.push(c[i].zhong);
        }
      }
    flag[t]=0;
  }
}     
int main()
{
  cin>>n>>m>>k>>e;
	for(int i=0;i<=n;++i)
	  f[i]=inf;
  for(int i=1;i<=e;++i)
    cin>>x[i]>>y[i]>>z[i];
  cin>>d;
  for(int j=1;j<=d;++j)
    cin>>num[j]>>xuan[j]>>ying[j];
  for(int i=1;i<=n;++i)
  {
  	for(int k=i;k<=n;++k)
  	{
	  	bool vis[50];
	  	for(int j=1;j<=m;++j)
	  	  vis[j]=1;
	  	vis[1]=vis[m]=1;
	  	jishu=0;
			for(int j=1;j<=m;++j)
			  dis[j]=inf;
	    memset(flag,0,sizeof(flag));
	    memset(head,0,sizeof(head));
	    for(int j=1;j<=d;++j)
	      if(!check(i,k,j))
	        vis[num[j]]=0;
	    for(int j=1;j<=e;++j)
	    {
	      if(vis[x[j]]&&vis[y[j]])
	      {
	        add(x[j],y[j],z[j]);
	        add(y[j],x[j],z[j]);
	      }
	    }
	    spfa(1);
	    cost[i][k]=dis[m];
    }
  }
  f[0]=0;
  for(int i=1;i<=n;++i)
    for(int j=0;j<=i-1;++j)
    {
    	if(cost[j+1][i]!=inf&&f[j]!=inf)
        f[i]=min(f[i],f[j]+cost[j+1][i]*(i-j)+k);
    }
  cout<<f[n]-k;
  return 0;
}