我的观点
将Tarjan中的每一个最大的边双(即和其他非子边双有交集的边双)内的所有点连上一个新点,然后新点连上旧店连着的点。然后重新建树。
其他更好的方法
大家有更好的方法吗?可以通过评论或者找我讨论,谢了各位!
代码
题目求2的x到y的路径中的简单环的个数次方。(每个点最多属于一个简单环)
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define N 100005
#define M 150005
#define LL long long
#define mo 1000000007
#define fo(i,a,b) for(i=a;i<=b;i++)
#define fd(i,a,b) for(i=a;i>=b;i--)
#define eg(i,x) for(i=head[x];i;i=edge[i].next)
using namespace std;
struct note{
int to,next,id;
};note edge[M*2];
int i,j,k,l,m,n,q,en,cnt;
int head[N];
int u,v,w,tot,tr_tot,temp;
int fa[N][18],sum[N][18];
int dfn[N],low[N],stack[N],dep[N];
int bel[N];
LL _2[N];
bool vis[M],bz[N],bz1[N];
void lb(int x,int y,int z){
edge[++tot].to=y;
edge[tot].next=head[x];
edge[tot].id=z;
head[x]=tot;
}
int getlca(int x,int y){
if (dep[x]<dep[y]) swap(x,y);
int i;
fd(i,17,0)if(dep[fa[x][i]]>=dep[y])temp+=sum[x][i],x=fa[x][i];
if (x==y) return x;
fd(i,17,0) if (fa[x][i]!=fa[y][i]){
temp+=sum[x][i]+sum[y][i];
x=fa[x][i],y=fa[y][i];
}
if (x!=y) temp+=sum[x][0]+sum[y][0];
return fa[x][0];
}
void tarjan(int x){
dfn[x]=low[x]=++tr_tot;
bz[x]=1;
stack[++stack[0]]=x;
int i;
eg(i,x)
if (!vis[edge[i].id] && edge[i].id<=m){
vis[edge[i].id]=1;
if (dfn[edge[i].to]==-1){
tarjan(edge[i].to);
low[x]=min(low[x],low[edge[i].to]);
} else
if (bz[edge[i].to]){
low[x]=min(low[x],dfn[edge[i].to]);
}
}
if (dfn[x]==low[x]){
if (stack[stack[0]]!=x)sum[x][0]=1;else sum[x][0]=0;
en=stack[0];
while (stack[stack[0]+1]!=x){
bel[stack[stack[0]]]=x;
if (stack[stack[0]]!=x) bz1[stack[stack[0]]]=1;
stack[0]--;
}
fo(i,stack[0]+2,en){
eg(j,stack[i]){
cnt++;
lb(x,edge[j].to,cnt);
lb(edge[j].to,x,cnt);
}
}
}
}
void dfs(int x){
int i;
eg(i,x)
if (!bz1[edge[i].to] && fa[x][0]!=edge[i].to
&& x!=edge[i].to){
fa[edge[i].to][0]=x;
dep[edge[i].to]=dep[x]+1;
dfs(edge[i].to);
}
}
int main(){
_2[0]=1;fo(i,1,100000) _2[i]=(_2[i-1]*2)%mo;
scanf("%d%d",&n,&m);
fo(i,1,m){
scanf("%d%d",&u,&v);
lb(u,v,i);lb(v,u,i);
}
memset(dfn,255,sizeof(dfn));
memset(low,255,sizeof(low));
tarjan(1);
dep[1]=1;
dfs(1);
fo(j,1,17)fo(i,1,n){
fa[i][j]=fa[fa[i][j-1]][j-1];
sum[i][j]=sum[fa[i][j-1]][j-1]+sum[i][j-1];
}
scanf("%d",&q);
fo(i,1,q){
scanf("%d%d",&u,&v);
u=bel[u],v=bel[v];
temp=0;
w=getlca(u,v);
temp+=sum[w][0];
printf("%lld\n",_2[temp]);
}
return 0;
}