题目
某一天,少年邂逅了同病相连的IA。见面后,IA一把牵起少年的手,决定和他一起逃离部落,离开这个无法容身的是非之地。
要逃离部落,少年和IA就需要先选择一条耗时最少的路线,从而避免被部落的大人们抓到。部落可以大致分为N个区域,少年和IA在区域1,部落的出口设在区域N。此外部落还有M条连接两个区域道路。道路是无向的,没有一条道路的两端连接相同的区域,也没有两条道路所连接的两个区域完全相同。对于其中前(M−1)条道路,其通过时间是确定的,但最后一条道路,由于地理因素,通过其的时间会不断变化。
现在,少年和IA得知了在K个不同的时段里,通过第M条道路的时间,请您分别计算出在这K 个时段中逃离部落的最少时间,以帮助他们确定行动的时刻。
分析
显然,答案有三种情况。设每次修改边(u,v)
①. 不经过(u,v)
②. S→u→v→T
③. S→v→u→T
如果当前要修改的边的权值为x,d(u,v)表示u,v在图中的最短路,那么答案
Ans=min{d(S,T), d(S,u)+x+d(v,T), d(S,v)+x+d(u,T)}
结合堆优化 Dijkstra,时间复杂度O(NlogN +M)。
我改了3个多小时,结果发现我以前的堆优化写错了。。。
学习了正确的姿势后马上就对了……
无力吐槽(ಥ﹏ಥ)
ps:这里是歌曲的地址
一首歌也能出题,%%%神牛
code
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<string>
#include<algorithm>
#include<stack>
#include<queue>
#include<vector>
using namespace std;
const long long maxx=1e16+10;
struct arr{
int x,y;
long long w;
int next;
}edge[2000000];
int ls[500000];
long long dis[500000];
long long dis1[500000];
long long dis2[500000];
int v[500000];
int n,m,k,s,t;
int edge_m=0;
struct cmp{
bool operator ()(int a,int b){
return dis[a]>=dis[b];
}
};
void add(int x,int y,long long w)
{
edge_m++;
edge[edge_m]=(arr){x,y,w,ls[x]};ls[x]=edge_m;
edge_m++;
edge[edge_m]=(arr){y,x,w,ls[y]};ls[y]=edge_m;
}
struct node
{
int num;
long long dis;
bool operator < (node b) const
{
return dis > b.dis;
}
};
priority_queue<node> Q;
void dij()
{
memset(v,0,sizeof(v));
dis[s]=0;
Q.push((node){s,0});
while (!Q.empty())
{
node x=Q.top();
Q.pop();
if (v[x.num]) continue;
v[x.num]=1;
//printf("%lld %d\n",dis2[x.num],x.num);
for (int j=ls[x.num];j;j=edge[j].next)
{
if (dis[edge[j].y]>dis[x.num]+edge[j].w)
{
dis[edge[j].y]=dis[x.num]+edge[j].w;
Q.push((node){edge[j].y,dis[edge[j].y]});
}
}
}
}
priority_queue<node> QQ;
void dij1()
{
memset(v,0,sizeof(v));
dis1[s]=0;
QQ.push((node){s,0});
while (!QQ.empty())
{
node x=QQ.top();
QQ.pop();
if (v[x.num]) continue;
v[x.num]=1;
//printf("%lld %d\n",dis2[x.num],x.num);
for (int j=ls[x.num];j;j=edge[j].next)
{
if (dis1[edge[j].y]>dis1[x.num]+edge[j].w)
{
dis1[edge[j].y]=dis1[x.num]+edge[j].w;
Q.push((node){edge[j].y,dis1[edge[j].y]});
}
}
}
}
priority_queue<node> QQQ;
void dij2()
{
memset(v,0,sizeof(v));
dis2[s]=0;
QQQ.push((node){s,0});
while (!QQQ.empty())
{
node x=QQQ.top();
QQQ.pop();
if (v[x.num]) continue;
v[x.num]=1;
//printf("%lld %d\n",dis2[x.num],x.num);
for (int j=ls[x.num];j;j=edge[j].next)
{
if (dis2[edge[j].y]>dis2[x.num]+edge[j].w)
{
dis2[edge[j].y]=dis2[x.num]+edge[j].w;
QQQ.push((node){edge[j].y,dis2[edge[j].y]});
}
}
}
}
long long min_(long long a,long long b)
{
if (a>=b) return b;
else return a;
}
int main()
{
freopen("monogatari.in","r",stdin);
freopen("monogatari.out","w",stdout);
scanf("%d%d%d",&n,&m,&k);
for (int i=1;i<=n;i++)
{
dis[i]=maxx;
dis1[i]=maxx;
dis2[i]=maxx;
}
for (int i=1;i<m;i++)
{
int x,y;
long long w;
scanf("%d%d%lld",&x,&y,&w);
add(x,y,w);
}
s=1; t=n;
int u,v;
/*for (int i=1;i<=m;i++)
{
if (edge[i].w==0) printf("%d %d %d",edge[i].w,edge[i].x,edge[i].y);
}*/
scanf("%d%d",&u,&v);
dij();
s=u;
dij1();
s=v;
dij2();
//int tt=t;
/*while (sum[tt])
{
printf("%d %lld\n",sum[tt],dis[tt]);
tt=sum[tt];
}
printf("%lld %lld %lld %lld",dis[13],dis[102],dis[260],dis[t]);*/
//printf("%lld %lld %lld %lld\n",dis[u],dis[v],dis1[t],dis2[t]);
for (int i=1;i<=k;i++)
{
long long w;
scanf("%lld",&w);
long long mi=min_(dis1[t]+dis[v]+w,dis2[t]+dis[u]+w);
//printf("%lld %lld %lld ",dis1[t]+dis[v]+w,dis2[t]+dis[u]+w,w);
if ((dis[t]==maxx)&&(mi>dis[t]))
{
printf("+Inf\n");
continue;
}
if (mi<dis[t])
printf("%lld\n",mi);
else
printf("%lld\n",dis[t]);
}
return 0;
}