Tarjan缩环的好方法

我的观点

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