题解
先从全都是文明城市情况考虑:显然就是 x x x 到 y y y 的路径上最大的边的最小值,直接重构树即可。
考虑非文明城市的加入:
因为从 x x x 到 y y y 的路径一定是经过若干个文明城市,所以这时路径的权值就定义为:两个文明城市直接的最短路。
然后同上,问题是这样的边是 n 2 n^2 n2 的,考虑从性质入手优化:
对于 M S T MST MST 上的一条边,对应原图上两个文明城市 i , j i,j i,j 的一条最短路的权值和,因为边 i → j i \to j i→j 在 M S T MST MST 上,所以有更优的选择使得不需要连 i → j i \to j i→j 而让 i , j i,j i,j 连通,所以记录每个点到所有文明城市的最短距离以及是哪个城市, i , j i,j i,j 就一定会被路径上的一条边枚举到。
注意,该题卡spfa
题解
#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int N=2e5+10;
char c[N];
LL d[N],val[N<<1];
int n,v[N],m,be[N],cnt,fa[N<<1],f[N<<1][22],dep[N<<1],now,Q;
int tot,ver[N<<1],fst[N],nxt[N<<1],w[N<<1];
inline void add(int x,int y,int z){ver[++tot]=y;nxt[tot]=fst[x];fst[x]=tot;w[tot]=z;}
vector<int> to[N<<1];
priority_queue<pair<int,int> >q;
struct Node{
int x,y;
LL w;
}e[N<<2];
bool cmp(const Node x,const Node y){return x.w<y.w;}
inline int get(int x){return (fa[x]==x)?x:fa[x]=get(fa[x]);}
void dfs1(int x){
for(auto y:to[x]){
dep[y]=dep[x]+1;
f[y][0]=x;
dfs1(y);
}
}
inline int LCA(int x,int y){
if(dep[x]<dep[y])swap(x,y);
int cha=dep[x]-dep[y];
for(int i=20;i>=0;--i)if(cha&(1<<i))x=f[x][i];
if(x==y)return x;
for(int i=20;i>=0;--i)if(f[x][i]!=f[y][i])x=f[x][i],y=f[y][i];
return f[x][0];
}
int main(){
scanf("%d%d",&n,&m);
scanf("%s",c+1);
for(int i=1;i<=m;++i){
int x,y,z;
scanf("%d%d%d",&x,&y,&z);
add(x,y,z);add(y,x,z);
}
for(int i=1;i<=n;++i)d[i]=1e17;
for(int i=1;i<=n;++i)if(c[i]=='1')q.push(make_pair(0,i)),v[i]=1,d[i]=0,be[i]=i;
while(q.size()){
int x=q.top().second,y=-q.top().first;q.pop();
if(d[x]!=y)continue;
// cout<<"FAQ "<<x<<" "<<d[x]<<" "<<be[x]<<endl;
v[x]=1;
for(int i=fst[x];i;i=nxt[i]){
int y=ver[i];
if(v[y])continue;
if(d[y]>d[x]+w[i]){
d[y]=d[x]+w[i];
be[y]=be[x];
q.push(make_pair(-d[y],y));
}
}
}
// for(int i=1;i<=n;++i)cout<<"FAQ "<<i<<" "<<d[i]<<" "<<be[i]<<endl;
for(int i=1;i<=n;++i){
for(int j=fst[i];j;j=nxt[j]){
int y=ver[j];
if(be[i]!=be[y])e[++cnt]=(Node){be[i],be[y],d[i]+d[y]+w[j]};
}
}
sort(e+1,e+cnt+1,cmp);
now=n;
for(int i=1;i<=n;++i)fa[i]=i;
for(int i=1;i<=cnt;++i){
int x=get(e[i].x),y=get(e[i].y);
if(x==y)continue;
now++;
fa[now]=now;
fa[x]=fa[y]=now;
to[now].push_back(x);
to[now].push_back(y);
val[now]=e[i].w;
}
dfs1(now);
for(int j=1;j<=20;++j)for(int i=1;i<=now;++i)f[i][j]=f[f[i][j-1]][j-1];
scanf("%d",&Q);
while(Q--){
int x,y;
scanf("%d%d",&x,&y);
printf("%lld\n",val[LCA(x,y)]);
}
return 0;
}