POJ 3177 Redundant Paths(Tarjan_BCC+缩点)

题意:

给你一个联通无向图,让你添加最少的边,使得这个图双联通。题意很清
晰,只要用tarjan_bcc算法缩点就能形成一棵树。
然后怎么计算最少的边呢?
我们用一种贪心的方法,我们想尽量让一条边能够让最多的点形成一个环
,那么我们就找LCA最远的一对叶节点,连起来,然后再找一对LCA最远
的,这样一直下去,就行了。如果叶节点为偶数,那么答案就是leaf/2,
如果为奇数,会多一个,那个只要和根节点连就好了,答案就是
(leaf+1)/2。统一起来就是(leaf+1)/2。

代码:

//
//  Created by  CQU_CST_WuErli
//  Copyright (c) 2015 CQU_CST_WuErli. All rights reserved.
//
// #include
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#define CLR(x) memset(x,0,sizeof(x))
#define OFF(x) memset(x,-1,sizeof(x))
#define MEM(x,a) memset((x),(a),sizeof(x))
#define ALL(x) x.begin(),x.end()
#define AT(i,v) for (auto &i:v)
#define BUG cout << "I am here" << endl
#define lookln(x) cout << #x << "=" << x << endl
#define SI(a) scanf("%d",&a)
#define SII(a,b) scanf("%d%d",&a,&b)
#define SIII(a,b,c) scanf("%d%d%d",&a,&b,&c)
#define Lson l,mid,rt<<1
#define Rson mid+1,r,rt<<1|1
#define Root 1,n,1
#define BigInteger bign
const int MAX_L=2005;// For BigInteger
const int INF_INT=0x3f3f3f3f;
const long long INF_LL=0x7fffffff;
const int MOD=1e9+7;
const double eps=1e-9;
const double pi=acos(-1);
typedef long long  ll;
using namespace std;

const int N=3e4+200;

int n,m;
int pnt[N],nxt[N],head[N];
int cnt;

void add_edge(int u,int v) {
    pnt[cnt]=v;nxt[cnt]=head[u];head[u]=cnt++;
}

int low[N],pre[N],bccno[N],in[N];
int dfs_clock,bcc_cnt;
stack<int> S;

void dfs(int fa,int u) {
    low[u]=pre[u]=++dfs_clock;
    S.push(u);in[u]=1;
    for (int i=head[u];~i;i=nxt[i]) {
        int v=pnt[i];
        if (i==(fa^1)) continue;
        if (!pre[v]) {
            dfs(i,v);
            low[u]=min(low[u],low[v]);
        }
        else if (in[v]) low[u]=min(low[u],pre[v]);
    }
    if (low[u]==pre[u]) {
        int v=-1;bcc_cnt++;
        do {
            v=S.top();S.pop();
            in[v]=0;bccno[v]=bcc_cnt;
        }while(u!=v);
    }
}

int deg[N];


int main(){
#ifdef LOCAL
    freopen("C:\\Users\\john\\Desktop\\in.txt","r",stdin);
//  freopen("C:\\Users\\john\\Desktop\\out.txt","w",stdout);
#endif
    while (SII(n,m)==2) {
        dfs_clock=bcc_cnt=0;while(S.size()) S.pop();
        CLR(low);CLR(pre);CLR(bccno);CLR(in);CLR(deg);
        OFF(head);cnt=0;
        for (int i=1;i<=m;i++) {
            int u,v;
            scanf("%d%d",&u,&v);
            add_edge(u,v);
            add_edge(v,u);
        }
        for (int i=1;i<=n;i++) if (!pre[i]) {
            dfs(-1,i);
        }
//      lookln(bcc_cnt);
        for (int u=1;u<=n;u++) {
            for (int i=head[u];~i;i=nxt[i]) {
                int v=pnt[i];
                if (bccno[u]!=bccno[v]) {
                    deg[bccno[u]]++;
                }
            }
        }
        int ans=0;
        for (int i=1;i<=n;i++) if (deg[i]==1) ans++;
        ans=(ans+1)/2 ;
        cout << ans << endl;
    }
    return 0;
}

你可能感兴趣的:(图论)