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;
}