1003: [ZJOI2006]物流运输trans
思路:最短路+DP
每一天都走最短路显然是不现实的
不过我们发现天数很小
设cost[i][j]——i到j天内都可以走的最短路距离(若不存在则为inf)
然后我们在进行一个DP
设f[i]——前i天运输的最小费用
显然,状态转移方程为:
f[i]=min(f[j]+cost[j+1][i]*(i-j)+K) (1<=j< i)
为了方便,初始状态为f[0]=-K
感觉做了好多事啊!
代码:
这应该是我改代码风格前的作品
难怪dijkstra都打错了……
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <iostream>
#define MIN(a,b) ((a)<(b)?(a):(b))
using namespace std;
const int MAXDAY=110,MAXN=25,INF=1000000;
struct Front
{
int x,y;
int num;
}front[2*MAXN*MAXN];
int orderx[2*MAXN*MAXN],stdx[MAXN],num[MAXN];
bool use[MAXDAY][MAXN];
bool check[MAXN];
int cost[MAXDAY][MAXDAY];
int f[MAXDAY];
int day,n,K,e,m,d,maxroad;
inline void add(int x,int y,int num)
{
front[++m].x=x;
front[m].y=y;
front[m].num=num;
return ;
}
inline void sortx()
{
for (int i=1;i<=m;i++)
num[front[i].x]++;
for (int i=1;i<=n;i++)
{
num[i]+=num[i-1];
stdx[i]=num[i];
}
for (int i=1;i<=m;i++)
orderx[++num[front[i].x-1]]=i;
return ;
}
inline int dijstra(int l,int r)
{
int ans[MAXN];
memset(check,false,sizeof(check));
for (int i=0;i<=n;i++)
ans[i]=INF;
ans[1]=0;
for (int i=1;i<n;i++)
{
int expend=0;
for (int j=1;j<=n;j++)
if (!check[j] && ans[j]<ans[expend])
expend=j;
check[expend]=true;
for (int j=stdx[expend-1]+1;j<=stdx[expend];j++)
{
int next=front[orderx[j]].y,COST=front[orderx[j]].num;
bool pass=true;
for (int k=l;k<=r;k++)
if (!use[k][next])
{
pass=false;
continue;
}
if (!pass)
continue;
ans[next]=MIN(ans[expend]+COST,ans[next]);
}
}
return ans[n];
}
inline void init()
{
cin>>day>>n>>K>>e;
for (int i=1;i<=e;i++)
{
int x,y,num;
scanf("%d%d%d",&x,&y,&num);
add(x,y,num);
add(y,x,num);
}
memset(use,true,sizeof(use));
sortx();
int d;
cin>>d;
while (d--)
{
int x,l,r;
scanf("%d%d%d",&x,&l,&r);
for (int i=l;i<=r;i++)
use[i][x]=false;
}
for (int i=1;i<=day;i++)
for (int j=i;j<=day;j++)
cost[i][j]=dijstra(i,j);
for (int i=1;i<=day;i++)
f[i]=INF;
f[0]=-K;
for (int i=1;i<=day;i++)
for (int j=0;j<i;j++)
f[i]=MIN(f[i],f[j]+cost[j+1][i]*(i-j)+K);
cout<<f[day];
return ;
}
int main()
{
init();
return 0;
}