先跑一下最短路图。
然后随便拎出来一条路径,把它叫做树边。
之后用不在这条路径的边去更新信息即可。
那么我们考虑这个路径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;
}