图论 无向图的双连通分量（模板） LA5135

#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <stack>
using namespace std;
typedef long long ll;
const int maxn=50005;
int pre[maxn],iscut[maxn],bccno[maxn],dfs_clock,bcc_cnt;
vector<int> G[maxn],bcc[maxn];

struct Edge
{
   int u,v;
};
stack<Edge> S;

int dfs(int u,int fa)//求无向图的双连通分量
{
   int lowu=pre[u]=++dfs_clock;
   int child=0;
   for(int i=0;i<G[u].size();i++)
   {
      int v=G[u][i];
      Edge e=(Edge){u,v};
      if(!pre[v])
      {
         S.push(e);
         child++;
         int lowv=dfs(v,u);
         lowu=min(lowu,lowv);
         if(lowv>=pre[u])
         {
            iscut[u]=true;
            bcc_cnt++;
            bcc[bcc_cnt].clear();
            for(;;)
            {
               Edge x=S.top();
               S.pop();
               if(bccno[x.u]!=bcc_cnt)
               {
                  bcc[bcc_cnt].push_back(x.u);
                  bccno[x.u]=bcc_cnt;
               }
               if(bccno[x.v]!=bcc_cnt)
               {
                  bcc[bcc_cnt].push_back(x.v);
                  bccno[x.v]=bcc_cnt;
               }
               if(x.u==u && x.v==v) break;
            }
         }
      }
      else if(pre[v]<pre[u] && v!=fa)
      {
         S.push(e);
         lowu=min(lowu,pre[v]);
      }
   }
   if(fa<0 && child==1) iscut[u]=0;
   return lowu;
}

void find_bcc(int n)
{
   memset(pre,0,sizeof(pre));
   memset(iscut,0,sizeof(iscut));
   memset(bccno,0,sizeof(bccno));
   dfs_clock=bcc_cnt=0;
   for(int i=0;i<n;i++)
      if(!pre[i]) dfs(i,-1);
}

int main()
{
   long long ans1=0,ans2=1;
   int kase=0,n,i,x,y;
   while(scanf("%d",&n)!=EOF)
   {
      ans1=0;ans2=1;
      if(n==0) break;
      while(!S.empty()) S.pop();
      for(i=0;i<50005;i++) bcc[i].clear();
      for(i=0;i<50005;i++) G[i].clear();
      for(i=0;i<n;i++)
      {
         scanf("%d%d",&x,&y);
         G[x].push_back(y);
         G[y].push_back(x);
      }
      find_bcc(n);
      for(i=1;i<=bcc_cnt;i++)
      {
         int cnt_cnt=0;
         for(int j=0;j<bcc[i].size();j++)
         {
            if(iscut[bcc[i][j]]) cnt_cnt++;
         }
         if(cnt_cnt==1)
         {
            ans1++;
            ans2*=(ll) (bcc[i].size()-cnt_cnt);
         }
      }
      if(bcc_cnt==1)
      {
         ans1=2;
         ans2=bcc[1].size()*(bcc[1].size()-1)/2;
      }
      printf("Case %d: %lld %lld\n",++kase,ans1,ans2);
   }
   return 0;
}