当前位置: 代码迷 >> 综合 >> bzoj-1003 物流运输
  详细解决方案

bzoj-1003 物流运输

热度:19   发布时间:2023-12-01 22:08:07.0

我用的是Dijkstra+DP(当然是参考大佬们的啦),Dijkstra用来进行预处理,处理出dis[start_time][end_time] --[start_time,end_time]时间段的1->m的最短路,其实很简单,我们首先要把在某个时间段不能装卸货物的码头给预处理出来,可以用前缀和来处理,具体看代码,所以当到达这样点时,我们直接continue。DP阶段,我们初始dp[0]=-k,为什么?因为 0 时间没有通过任何码头,并且我们到第一次运行的航线还要交换一次航线,所以是-k,因为要抵消交换一次航线的代价。转移方程 是:dp[i]=min(dp[i],dp[j]+dis[j+1][i]*(i-j)+k)。
还有就是 当dis[j+1][i]是INF时,不要进行转移(INF很大的情况,我的INF=1e18,习惯问题),因为爆long long了。

#include<cstdio>
#include<vector>
#include<algorithm>
#include<queue>
using namespace std;typedef long long ll;
typedef pair<ll,int>P;
const int maxn=100;
const int day=200;
const ll INF=1e18;struct Edge{int to;ll len;Edge(int to_=0,ll len_=0){to=to_;len=len_;}
};
vector<Edge>G[maxn];
ll distan[maxn];
ll dis[day][day];
ll dp[day];
int notok[maxn][day];
bool vis[maxn];
int n,m,k,e,d;void Dijkstra(int stime,int etime){for(int i=0;i<maxn;i++) distan[i]=INF,vis[i]=false;distan[1]=0;priority_queue<P,vector<P>,greater<P> >que;que.push(P(0,1));int u,v,len;while(!que.empty()){P p=que.top();que.pop();u=p.second;if(vis[u]) continue;vis[u]=true;len=G[u].size();for(int i=0;i<len;i++){Edge e=G[u][i];v=e.to;if(notok[v][etime]-notok[v][stime-1]>0) continue;if(!vis[v]&&distan[v]>distan[u]+e.len){distan[v]=distan[u]+e.len;que.push(P(distan[v],v));}}}dis[stime][etime]=distan[m];
}int main(){scanf("%d%d%d%d",&n,&m,&k,&e);for(int i=0;i<e;i++){int u,v;ll len;scanf("%d%d%lld",&u,&v,&len);G[u].push_back(Edge(v,len));G[v].push_back(Edge(u,len));} scanf("%d",&d);for(int i=0;i<d;i++){int t,stime,etime;scanf("%d%d%d",&t,&stime,&etime);for(int j=stime;j<=etime;j++) notok[t][j]=1;}for(int i=1;i<=m;i++){for(int j=1;j<=n;j++){notok[i][j]+=notok[i][j-1];}}for(int i=1;i<=n;i++){for(int j=1;j<=i;j++){Dijkstra(j,i);}}dp[0]=(ll)-k;for(int i=1;i<=n;i++){dp[i]=INF;for(int j=0;j<i;j++){if(dis[j+1][i]!=INF) dp[i]=min(dp[i],dp[j]+dis[j+1][i]*(ll)(i-j)+(ll)k);}}printf("%lld\n",dp[n]);
}