树和基环树上两点间最短距离问题

总时间限制: 

20000ms

单个测试点时间限制: 

1000ms

内存限制: 

524288kB

描述

给定 n (2≤n≤2*105) 个点 m (m=n 或 m=n-1) 条边的无向连通图，回答 q 组询问，每组询问给出点 u 和点 v，请计算 u 和 v 的最短距离。两点间的最短距离，是指连接它们的路径的最少边数。

输入

第一行包含两个整数 n, m (2≤n≤2*105, m=n 或 m=n-1）分别表示点数和边数。
接下去 m 行，每行两个整数 a, b (1≤a,b≤n, a≠b) 表示点 a 和点 b 之间连一条无向边。数据保证 n 个点构成连通图，且无自环、无重边。
下一行包含一个整数 q (1≤q≤2*105) 表示询问个数。
接下去 q 行，每行两个整数 u, v (1≤u,v≤n, u≠v) 表示一组询问。

输出

共 q 行，每行表示一个询问的答案，即点 u 到点 v 的最短距离。

样例输入

样例输入1
6 5 
1 3 
2 3 
2 5 
3 4 
4 6
5
1 3
3 5
2 6
5 6
4 2

样例输入2
6 6 
1 3 
2 3 
2 5 
3 4 
4 5 
4 6
5
1 6
3 4
2 4
1 5
6 5

样例输出

样例输出1
1
2
3
4
2

样例输出2
3
1
2
3
2

对于普通树的询问，我们只需要找到它们两点的LCA，可得最短距离为d[ x ] + d[ y ] - 2*d[ lca ]

对于基环树讲，我们可以将环上的每个点都当作根节点来对其子树进行一次dfs做预处理（同普通树）并记录所属根节点，在询问时，如果两点属于同一根节点的子树，那么按照普通树的询问来处理，否则我们将环上若干个相邻的点进行编号为1……tot，那么我们就可以得到两个根之间有两个走法，第一个是dis = abs( id[ x ] - id[ y ])，第二个则是tot-dis，最后取两者最小值加上d[ x ] + d[ y ]就是最短距离了

我的代码中找环使用的是tarjan，当然也可以拓扑，最后拓扑完degree>=2就在环上

代码：

#include<bits/stdc++.h>
#define ll long long
#define pii pair<int,int>
#define endl '\n'
using namespace std;
const int maxn=2e5+10,L=20;
int n,m,fa[maxn][L],d[maxn],q,k;
int deg[maxn],id[maxn],tot;
int bl[maxn],rt;
bool on[maxn],fst,vis[maxn];
vector<int> e[maxn];
void dfs(int x,int pa){
	fa[x][0]=pa,d[x]=d[pa]+1,bl[x]=rt;
	for(int i=1;i<=k;i++) fa[x][i]=fa[fa[x][i-1]][i-1];
	for(int y:e[x]) if(!on[y]&&y!=pa) dfs(y,x);
}
inline int lca(int x,int y){
	if(d[x]<d[y]) swap(x,y);
	for(int i=k;~i;i--) if(d[fa[x][i]]>=d[y]) x=fa[x][i];
	if(x==y) return x;
	for(int i=k;~i;i--) if(fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
	return fa[x][0];
}
void s1(){
	rt=1;
	dfs(1,0);
	cin>>q;
	while(q--){
		int x,y;cin>>x>>y;
		int l=lca(x,y);
		cout<<d[x]+d[y]-2*d[l]<<endl;
	}
}
int s,ed,pre[maxn],que[maxn];
void tarjan(int x,int pa){
	vis[x]=1,pre[x]=pa;
	for(int y:e[x]) if(y!=pa){
		if(vis[y]) s=x,ed=y;
		else tarjan(y,x);
		if(s) return;
	}
}
void s2(){
	d[0]=-1;
	tarjan(1,0);
	int cur=s;
	while(1){
		on[cur]=1,id[cur]=tot;
		que[tot++]=cur;
		if(cur==ed) break;
		cur=pre[cur];
	}
	for(int i=0;i<tot;i++){
		rt=que[i];
		dfs(rt,0);
	}
	cin>>q;
	while(q--){
		int x,y;cin>>x>>y;
		if(bl[x]==bl[y]){
			int l=lca(x,y);
			cout<<d[x]+d[y]-2*d[l]<<endl;
		}
		else{
			int dis=abs(id[bl[x]]-id[bl[y]]);
			cout<<min(dis,tot-dis)+d[x]+d[y]<<endl;
		}
	}
}
int main(){
	ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
	cin>>n>>m;k=log2(n);
	for(int i=1;i<=m;i++){
		int u,v;cin>>u>>v;
		e[u].push_back(v);
		e[v].push_back(u);
		++deg[u],++deg[v];
	}
	if(n==m) s2();
	else s1();
	return 0;
}