LA 3514

LA 3514

題意是給一個無向圖，若它是仙人掌，則求出它有多少生成子圖(包括自身)是仙人掌，否則輸出0。仙人掌被定義爲每條邊最多在一個簡單迴路上的連通無向圖。
關於Cactus問題，這裏有篇08年國家集訓隊周冬的論文，比較詳細地總結了所有的仙人掌問題。
以下是AC代碼:

#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define debug puts("Infinity is awesome!")
#define mm(a,b) memset(a,b,sizeof(a))
#define LS (root<<1)
#define RS (root<<1|1)
#define LSON LS,l,mid
#define RSON RS,mid+1,r
#define LL long long
const int Inf=1e9+7;
const int maxn=2e4+5;
vector<int> G[maxn];
struct solotion{
    int n, cylcnt;
    vector<int> G[maxn];
    int c[maxn], dfn[maxn];
    int ans[maxn], cycle[maxn];
    void init(int n){
        this->n=n;
        cylcnt=0;
        for(int i=0;i<=n;i++) G[i].clear();
        mm(c,0); mm(dfn,0);
        mm(cycle,0);
        mm(ans,0);
    }
    void add_edge(int u,int v){
        G[u].push_back(v);
        G[v].push_back(u);
    }
    void dfs(int u,int fa){
        for(int i=0;i<G[u].size();i++){
            int v=G[u][i];
            if(v==fa) continue;
            if(!dfn[v]){
                dfn[v]=dfn[u]+1;
                dfs(v, u);
                c[u]+=c[v];
            }else if(dfn[v]<dfn[u]){
                cycle[cylcnt++]=dfn[u]-dfn[v]+2;
                c[u]++, c[v]--;//此反向邊連向v,c[v]自減,防止重複加
            }
        }
    }
    void get_ans(){
        int len=0;
        ans[len]=1;
        for(int i=0;i<cylcnt;i++){
//            printf("cycle=%d\n",cycle[i]);
            for(int j=0;j<=len;j++)
                ans[j]*=cycle[i];
            for(int j=0;j<=len;j++){
                ans[j+1]+=ans[j]/10;
                ans[j]%=10;
            }
            while(ans[len+1]){
                ans[len+2]+=ans[len+1]/10;
                ans[++len]%=10;
            }
        }
        for(int i=len;i>=0;i--)
            printf("%d",ans[i]);
        puts("");
    }
    bool solve(){
        dfn[1]=1; dfs(1,-1);
        for(int i=1;i<=n;i++)
        if(dfn[i]==0||c[i]>1){
            puts("0");
            return false;
        }
        get_ans();
        return true;
    }
};
solotion sol;
int main(){
    int n, m, K;
    int u, v;
    int cas=0;

    while(~scanf("%d%d",&n, &m)){
        sol.init(n);
        if(cas++!=0) puts("");
        for(int i=0;i<m;i++){
            scanf("%d",&K);
            scanf("%d",&u);
            for(int i=1;i<K;i++){
                scanf("%d",&v);
                sol.add_edge(u,v);
                u=v;
            }
        }
        sol.solve();
    }
    return 0;
}