bzoj1003: [ZJOI2006]物流运输(DP+spfa)

1003: [ZJOI2006]物流运输

题目:传送门

题解:

   可以用spfa处理出第i天到第j都走这条路的花费,记录为cost

   f[i]表示前i天的最小花费:f[i]=min(f[i],f[j-1]+cost*(i-j+1)+k);

水一发代码:

 #include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define qread(x) x=read()
using namespace std;
inline int read()
{
int f=,x=;char ch;
while(ch<'' || ch>''){if(ch=='-')f=-;ch=getchar();}
while(ch>='' && ch<=''){x=x*+ch-'';ch=getchar();}
return f*x;
}
int n,m,k,e,D;
struct edge
{
int x,y,d,next;
}a[];int len,last[];
void ins(int x,int y,int d)
{
len++;a[len].x=x;a[len].y=y;a[len].d=d;
a[len].next=last[x];last[x]=len;
}
int flag[][],f[],d[],list[];
int spfa(int st,int ed)
{
for(int i=;i<=m;i++)d[i]=;
int head=,tail=;list[]=;d[]=;
while(head!=tail)
{
int x=list[head];
for(int k=last[x];k;k=a[k].next)
{
int y=a[k].y;
if((flag[y][ed]-flag[y][st-]==) && (d[y]>d[x]+a[k].d))
{
d[y]=d[x]+a[k].d;
list[tail++]=y;
if(tail==m+)tail=;
}
}
head++;if(head==m+)head=;
}
return d[m];
}
int main()
{
qread(n);qread(m);qread(k);qread(e);
len=;memset(last,,sizeof(last));
for(int i=;i<=e;i++)
{
int x,y,d;
qread(x);qread(y);qread(d);
ins(x,y,d);ins(y,x,d);
}
scanf("%d",&D);
memset(flag,,sizeof(flag));
for(int i=;i<=D;i++)
{
int p,x,y;
qread(p);qread(x);qread(y);
for(int j=x;j<=y;j++)flag[p][j]=;
}
for(int i=;i<=m;i++)
for(int j=;j<=n;j++)
flag[i][j]+=flag[i][j-];
f[]=;
for(int i=;i<=n;i++)
{
f[i]=;
for(int j=;j<=i;j++)
{
int cost=spfa(j,i);
f[i]=min(f[i],f[j-]+cost*(i-j+)+k);
}
}
printf("%d\n",f[n]-k);
return ;
}
上一篇:动易CMS - 设为首页代码和加入收藏代码(兼容各种浏览器)


下一篇:动易CMS之标签管理