Forest Program HDU - 6736 dfs找无向图所有环(长度)

题意：

给定一个仙人掌，删去一些边可以让它变成一个森林（一棵树也是森林），求方案数。 n≤300000,m≤500000

思路：

对于一个有n条边的环，可以删去1、2…n条边，则共有C(1,N) + C(1,2) + … + C(1,N)=2^n-1 种方案；

对于不成环的边（若有p条），每一条可删可不删，即两种选择，则有2^p种方案

那么最终的答案就是每个环的方案数相乘，再乘上非成环的边的方案数

所以问题现在转变为找无向图中所有的环（长度），和非成环的边数；

代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
using namespace std;
typedef long long ll;
const ll mod=998244353;
vector<int>a[300010];
int v[300010],fa[300010],d[300010],num=0;
ll cnt[300010];
void dfs(int u)
{
    v[u]=1;
    d[u]=d[fa[u]]+1;
    for(int i=0; i<a[u].size(); i++)
    {
        int vv=a[u][i];
        if(fa[u]==vv)
            continue;
        if(!v[vv])
        {
            fa[vv]=u;
            dfs(vv);

        }
        else if(d[vv]<d[u])///找到起点
        {
            cnt[++num]=d[u]-d[vv]+1;
            cnt[0]+=cnt[num];

        }
    }
}
ll mul(ll a, ll b)
//求a*b%m
{
    ll ans = 0;
    a %= mod;
    while(b)
    {
        if(b & 1)
            ans = (ans + a) % mod;
        b /= 2;
        a = (a + a) % mod;
    }
    return ans;
}
ll pp(ll a, ll b)
//a^b % m
{
    ll ans = 1;
    a %= mod;
    while(b)
    {
        if(b & 1)
            ans = mul(a, ans);
        b /= 2;
        a = mul(a, a);
    }
    ans %= mod;
    return ans;
}
int main()
{
    int n,m,u,vv;
    while(~scanf("%d%d",&n,&m))
    {

        for(int i=0; i<=n; i++)
            a[i].clear();
        num=0;
        memset(v,0,sizeof(v));
        memset(fa,0,sizeof(fa));
        memset(d,0,sizeof(d));
        memset(cnt,0,sizeof(cnt));
        ll ans=1;
        for(int i=1; i<=m; i++)
        {
            scanf("%d%d",&u,&vv);
            a[u].push_back(vv);
            a[vv].push_back(u);
        }
        for(int i=1; i<=n; i++)
        {
            if(!v[i])
                dfs(i);
        }
        ll z;

        for(int i=1; i<=num; i++)
        {
            z=pp(2,cnt[i])-1;
            ans=mul(ans,z);
        }
        ll s=m-cnt[0];
        if(s>0)
        {
            z=pp(2,s);
            ans=mul(ans,z);
        }
        printf("%lld\n",ans);
    }
    return 0;
}