POJ 3177 Redundant Paths

http://poj.org/problem?id=3177

先找双连通分量缩点，然后找桥重新构图成树，然后找树的叶子节点

答案是：(树的叶子节点+1)/2

顺便附上双连通分量的模板

//#pragma comment(linker, "/STACK:102400000,102400000")
#include<cstdio>
#include<iostream>
#include<cstring>
#include<string>
#include<cmath>
#include<set>
#include<list>
#include<map>
#include<iterator>
#include<cstdlib>
#include<vector>
#include<queue>
#include<stack>
#include<algorithm>
#include<functional>
using namespace std;
typedef long long LL;
#define ROUND(x) round(x)
#define FLOOR(x) floor(x)
#define CEIL(x) ceil(x)
const int maxn=5010;
const int maxm=20010;
const int inf=0x3f3f3f3f;
const LL inf64=0x3f3f3f3f3f3f3f3fLL;
const double INF=1e30;
const double eps=1e-6;

/**
*双连通分量：Tarjan算法($O(V+E)$)
*输入：图(链式前向星),n
*输出：bcc[](双连通分量)
*/
struct Edge
{
    int v;
    int next;
}edge[maxm];
int head[maxn],edgeNum;
void addSubEdge(int u,int v)
{
    edge[edgeNum].v=v;
    edge[edgeNum].next=head[u];
    head[u]=edgeNum++;
}
void addEdge(int u,int v)
{
    addSubEdge(u,v);
    addSubEdge(v,u);
}
int n,m;
int mark[maxn],low[maxn];    //mark初始化为0
int sstack[maxn],stop;        //top初始化为0
int dfn;                    //初始化为1
int bcc[maxn],bid;            //bid初始化为1
int dfs(int k,int p)
{
    int i,j;
    low[k]=mark[k]=dfn++;
    sstack[stop++]=k;
    for(i=head[k]; i!=-1; i=edge[i].next)
    {
        j=edge[i].v;
        if(mark[j]==0)
        {
            dfs(j,k);
            low[k]=min(low[k],low[j]);
        }
        else if(j!=p) low[k]=min(low[k],mark[j]);
    }
    if(mark[k]>low[k])return 0;
    while(sstack[--stop]!=k)// 导出一个双连通分量
    {
        bcc[sstack[stop]]=bid;
    }
    bcc[k]=bid;
    ++bid;
    return 0;
}
/*Dfs调用*/
void tarjan()
{
    memset(mark,0,sizeof(mark));
    dfn=bid=1;
    stop=0;
    for(int i=1; i<=n; ++i) if(mark[i]==0) dfs(i,-1);
}

Edge tredge[maxm];
int trhead[maxn],tredgeNum;
void addTrSubEdge(int u,int v)
{
    tredge[tredgeNum].v=v;
    tredge[tredgeNum].next=trhead[u];
    trhead[u]=tredgeNum++;
}
void addTrEdge(int u,int v)
{
    addTrSubEdge(u,v);
    addTrSubEdge(v,u);
}

int low2[maxn];
int mark2[maxn];
int dfn2;
void dfs2(int k,int p)
{
    mark2[k]=low2[k]=dfn2++;
    for(int i=head[k]; i!=-1; i=edge[i].next)
    {
        int j=edge[i].v;
        if(j==p) continue;
        if(mark2[j]==0)
        {
            dfs2(j,k);
            low2[k]=min(low2[k],low2[j]);
            if(low2[j]==mark2[j])
            {
                addTrEdge(bcc[k],bcc[j]);
                //发现(k,j)为桥
            }
        }
        else low2[k]=min(low2[k],mark2[j]);
    }
}
void bridge()
{
    memset(mark2,0,sizeof(mark2));
    dfn2=1;
    for(int i=1;i<=n;i++) if(mark2[i]==0) dfs2(i,-1);
}

void init()
{
    memset(head,-1,sizeof(head));
    edgeNum=0;
    memset(trhead,-1,sizeof(trhead));
    tredgeNum=0;
}
void input()
{
    scanf("%d%d",&n,&m);
    for(int i=0;i<m;i++)
    {
        int u,v;
        scanf("%d%d",&u,&v);
        addEdge(u,v);
    }
}
void solve()
{
    tarjan();
//    for(int i=1;i<=n;i++) cout<<bcc[i]<<" ";
//    cout<<endl;
//    cout<<bid<<endl;
    if(bid<=2)
    {
        puts("0");
        return;
    }
    bridge();
    int ans=0;
    for(int i=1;i<bid;i++)
    {
        int cnt=0;
        for(int j=trhead[i];j!=-1;j=tredge[j].next,cnt++);
        if(cnt==1) ans++;
    }
    printf("%d
",(ans+1)/2);
}
int main()
{
    std::ios_base::sync_with_stdio(false);
//    freopen("in.cpp","r",stdin);
    init();
    input();
    solve();
    return 0;
}