六兆年と一夜物語 jzoj3501 最短路

题目

某一天，少年邂逅了同病相连的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;}