bzoj2725

先跑一下最短路图。
然后随便拎出来一条路径,把它叫做树边。
之后用不在这条路径的边去更新信息即可。
那么我们考虑这个路径P。
其中(u,v)在路径P上。
如果我们不删除(u,v)而是删除非路径上的边,那么我们就会直接输出。
如果(u,v)在路径上,那么我们肯定要选这种东西去规避(u,v):
原路径是:S->x’->..->(u->v)->..->y’->T
现在的路径:S->x’->x->y->y’->T
其中,x,y在最短路图上。
这样的话我们就绕过了(u,v)对吧。
定义front[x]为从S开始走最短路走到x准备在哪个点离开原本的最短路。
同理定义next[x]是回到最短路。
那么我们处理每一条非树边,显然它对答案的贡献:
把是树边且开头在[pos[front[x]],pos[next[y]] - 1]这段区间里所有答案都取min{ds[x] + dt[y] + val(x -> y)]}。
问题在于怎么去算front[]数组什么的。
1.最短路(树)的建立:
肯定是在最短路取等号的时候。
我们随便用dfs拎出一条路径之后break掉。
2.front/next数组的建立:
把在路径p上的点u塞到队列里,那么u能遍历到的在最短路图上的点v的front都更新成x的front。注意一件事情,就是如果front[v]已经有了,那么肯定你就不要乱更新了,因为前面的那个比x优秀。
剩下就没啥好说的了hh

#include <bits/stdc++.h>
#define Rep(i,n) for(int i = 1;i <= n;++ i)
#define Dwn(i,n) for(int i = n;i;-- i)
#define v edge[i].to
#define RepG(i,x) for(int i = head[x];~ i;i = edge[i].next)
using namespace std;
const int N = 131072 * 2;
const int Mod = 131072 * 2 - 1;
typedef long long LL;
const LL inf = 1ll << 62;
int n,m,head[N],cnt,S,T,Q,pos[N],ft[N],fs[N],q[N],tim,Rank[N];
LL ds[N],dt[N],ans[N];
bool vis[N];
struct Edge{int next,to,w;}edge[N << 1];
void save(int a,int b,int c){edge[cnt] = (Edge){head[a],b,c},head[a] = cnt ++;}
priority_queue<pair<LL,int> > QAQ;
int read()
{
    char ch = getchar();
    while(ch > '9' || ch < '0')ch = getchar();
    int x = 0;
    while(ch >= '0' && ch <= '9')x = 10 * x + ch - '0',ch = getchar();
    return x;
}
int Power(int x,int p,int Mod)
{
    int cur = 1;
    while(p)
    {
        if(p & 1)cur = 1ll * cur * x % Mod;
        p >>= 1;
        x = 1ll * x * x % Mod;
    }
    return cur;
} 
void spfa(int x,LL *dis)
{
    Rep(i,n)vis[i] = 0;
    QAQ.push(make_pair(0ll,x));
    dis[x] = 0;
    while(!QAQ.empty())
    {
        x = QAQ.top().second;QAQ.pop();
        if(vis[x])continue;
        vis[x] = 1;
        RepG(i,x)
        {
            if(dis[v] > dis[x] + edge[i].w)
            {
                dis[v] = dis[x] + edge[i].w;
                QAQ.push(make_pair(-dis[v],v));
            }
        }
    }
}
bool in(int x,int y){return pos[x] && pos[y] && (abs(pos[x] - pos[y]) == 1);}
void dfs(int x,int fa = -1)//建最短路树
{
    pos[x] = ++ tim;
    Rank[tim] = x;
    if(x == T)return;
    RepG(i,x)if(v != fa && (ds[x] + edge[i].w + dt[v] == ds[T])){dfs(v,x);break;}
}
void Rebfs(int x,LL *dis,int *fr)
{
    int h = 0,t = 0;
    q[0] = x;
    fr[x] = x;
    while(h <= t)
    {
        x = q[h ++];
        RepG(i,x)
        {
            if(!pos[v] && !fr[v] && dis[v] == dis[x] + edge[i].w)//如果这个玩意在最短路图&&不在最短路树上
                fr[q[++ t] = v] = fr[x];
        }
    }
}
//~~~~~~~~~~~~~~~~~~~~~~~~~~~Segment Beats~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
#define u t[x]
#define mid (l + r >> 1)
#define lson x << 1,l,mid
#define rson x << 1 | 1,mid + 1,r
struct Seg{LL cov,val;Seg(){cov = val = inf;}}t[N << 2];
LL ccm;int ql,qr;
LL min(const LL &x,const LL &y){return x > y ? y : x;}
void Rlx(LL &x,LL y){x = min(x,y);}
void Upd(int x){Rlx(u.val,t[x << 1].val);Rlx(u.val,t[x << 1 | 1].val);}
void Dw(int x)
{
    LL &st = u.cov;
    Rlx(t[x << 1].val,st);Rlx(t[x << 1 | 1].val,st);
    Rlx(t[x << 1].cov,st);Rlx(t[x << 1 | 1].cov,st);
    st = inf;
}
void Modify(int x,int l,int r)
{
    if(l >= ql && r <= qr)
    {
        u.cov = min(u.cov,ccm);
        u.val = min(u.val,ccm);
        return;
    }
    if(u.cov != inf)Dw(x);
    if(ql <= mid)Modify(lson);
    if(mid < qr)Modify(rson);
    Upd(x);
}
void dfs(int x,int l,int r)
{
    if(l == r){ans[l] = u.val;return;}
    Dw(x);
    dfs(lson);dfs(rson);
}
//~~~~~~~~~~~~~~~~~~~~Segment Beats Ended~~~~~~~~~~~~~~~~~~~~~~~~~~~~
int e[N][3];
int main()
{
    scanf("%d%d",&n,&m);
    memset(head,-1,sizeof(head));
    Rep(i,m)
    {
        int a,b,c;
        //scanf("%d%d%d",&a,&b,&c);
        a = read(),b = read(),c = read();
        save(a,b,c),save(b,a,c);
        e[i][0] = a,e[i][1] = b,e[i][2] = c;
    }
    //fs[x]表示从S到x最早在什么地方离开的
    scanf("%d%d",&S,&T);
    Rep(i,n)ds[i] = dt[i] = 1ll << 62;
    spfa(S,ds);spfa(T,dt);
    dfs(S,-1);
    Rep(i,tim)Rebfs(Rank[i],ds,fs);
    Dwn(i,tim)Rebfs(Rank[i],dt,ft);
    Rep(i,m)
    {
        int a = e[i][0],b = e[i][1],c = e[i][2];
        if(!fs[a] || !ft[b] || in(a,b))continue;
        if(pos[fs[a]] < pos[ft[b]]){ql = pos[fs[a]],qr = pos[ft[b]] - 1;ccm = ds[a] + dt[b] + c;Modify(1,1,tim - 1);}
        swap(a,b);
        if(pos[fs[a]] < pos[ft[b]]){ql = pos[fs[a]],qr = pos[ft[b]] - 1;ccm = ds[a] + dt[b] + c;Modify(1,1,tim - 1);}
    }
    dfs(1,1,tim - 1);
    scanf("%d",&Q);
    while(Q --)
    {
        int a,b;
        //scanf("%d%d",&a,&b);
        a = read(),b = read();
        if(!in(a,b))printf("%lld\n",ds[T]);
        else 
        {
            if(pos[a] > pos[b])swap(a,b);
            if(ans[pos[a]] >= inf)puts("Infinity");
            else printf("%lld\n",ans[pos[a]]);
        }
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值