[CF231E]Cactus

Cactus

题解

很明显，每当经过一个边双连通分量时，如果这个分量里的点点数大于2，就会存在两种走法。

于是，我们很快就想到了缩点。将每个连通分量缩成一个点，size大于2的点权为1，否则为0。

通过lca求出两点间的路径上权值的和即可。

当然，我们也可以将所有询问离线下来，dfs时通过并查集维护lca，这样可以将时间复杂度降到，相比原复杂度减少了许多，当然，两种方法都不会T。

源码

#include<bits/stdc++.h>
using namespace std;
#define MAXN 200005
typedef long long LL;
typedef unsigned long long uLL;
const LL mo=1e9+7;
typedef pair<int,int> pii;
template<typename _T>
_T Fabs(_T x){return x<0?-x:x;}
template<typename _T>
void read(_T &x){
	_T f=1;x=0;char s=getchar();
	while(s>'9'||s<'0'){if(s=='-')f=-1;s=getchar();}
	while('0'<=s&&s<='9'){x=(x<<3)+(x<<1)+(s^48);s=getchar();}
	x*=f;
}
int n,m,q,head[MAXN],tot,idx,dep[MAXN],siz[MAXN],bin[MAXN];
int cnt,belong[MAXN],U[MAXN],V[MAXN],fa[MAXN],ans[MAXN];
int dfn[MAXN],low[MAXN],vis[MAXN],Head[MAXN],Tot;
int sta[MAXN],stak;bool insta[MAXN];
vector<pii>G[MAXN];
struct edge{int to,nxt,flag;}e[MAXN<<1],E[MAXN<<1];
void addEdge(int u,int v){e[++tot]=(edge){v,head[u]};head[u]=tot;}
void addedge(int u,int v){E[++Tot]=(edge){v,Head[u]};Head[u]=Tot;}
int add(int x,int y){return x+y>=mo?x+y-mo:x+y;}
void makeSet(int x){for(int i=1;i<=x;i++)fa[i]=i;}
int findSet(int x){return fa[x]==x?x:fa[x]=findSet(fa[x]);}
void unionSet(int a,int b){
	int u=findSet(a),v=findSet(b);
	if(u==v)return ;fa[u]=v;
} 
void tarjan(int u,int fa){
	dfn[u]=low[u]=++idx;sta[++stak]=u;insta[u]=1;
	for(int i=head[u];i;i=e[i].nxt){
		int v=e[i].to;if(v==fa)continue;
		if(!dfn[v])tarjan(v,u),low[u]=min(low[v],low[u]);
		else if(insta[u])low[u]=min(low[u],dfn[v]);
	}
	if(dfn[u]==low[u]){
		cnt++;int v;
		do{
			v=sta[stak--];siz[cnt]++;
			belong[v]=cnt;insta[v]=0;
		}while(u!=v);
	}
}
void getDis(int u,int dp){
	vis[u]=1;dep[u]=dp;
	for(int i=Head[u];i;i=E[i].nxt){
		int v=E[i].to;if(vis[v])continue;
		getDis(v,dp+(siz[v]>2));
	}
}
void solve(int u){
	vis[u]=1;
	for(int i=0;i<G[u].size();i++){
		int v=G[u][i].first,id=G[u][i].second;
		if(vis[v]){
			int com=findSet(v);
			ans[id]=dep[u]+dep[v]-2*dep[com];
			//printf("%d %d %d\n",dep[u],dep[v],dep[com]);
			if(siz[com]>2)ans[id]++;
		}
	}
	for(int i=Head[u];i;i=E[i].nxt){
		int v=E[i].to;if(vis[v])continue;
		solve(v);unionSet(v,u);
	}
}
signed main(){
	read(n);read(m);
	for(int i=1;i<=m;i++){
		int u,v;read(u);read(v);
		addEdge(u,v);addEdge(v,u);
		U[i]=u;V[i]=v;;
	}
	read(q);tarjan(1,1);memset(vis,0,sizeof(vis));
	for(int i=1;i<=m;i++)
		if(belong[U[i]]!=belong[V[i]])
			addedge(belong[U[i]],belong[V[i]]),
			addedge(belong[V[i]],belong[U[i]]);
	for(int i=0;i<q;i++){
		int u,v;read(u);read(v);
		G[belong[u]].push_back(make_pair(belong[v],i));
		G[belong[v]].push_back(make_pair(belong[u],i));
	}
	memset(vis,0,sizeof(vis));getDis(1,siz[1]>2);
	memset(vis,0,sizeof(vis));makeSet(n);solve(1);
	bin[0]=1;for(int i=1;i<=(int)1e5;i++)bin[i]=bin[i-1]*2%mo;
	for(int i=0;i<q;i++)printf("%d\n",bin[ans[i]]);
	return 0;
}

谢谢！！！