hdu 4582 (树上的贪心)

题意：一张无向图中告诉你一个dfs树，还有若干反向边。问你如何选取最小的边使得所有只包含一条反向边的环被覆盖。

思路：我们把有反向边的点与反向边连接的点在树上的路径称作一个区间，可以想到了从叶子结点开始向上到一个区间到不得不选边覆盖的时候在选边，同样能选的边越靠近根节点越好。关键是在程序实现上，自己实现的时候用的是vector<bool>来标记边的状态果断超时，网上搜到一种用bitset，以前从来没用过太神了，代码又短思路又清晰参考http://hi.baidu.com/zjut_dd/item/5d3543e55ba0ebebe0a5d4d7

主要就是从低至上每次遇到反向边则记录，若是当前节点到其父亲结点必须选了则选上（选上了则代表当前子树的所有结点连接的反向边全部搞定了），否则的话则向上一层传递还未选择的区间。

代码如下：

/**************************************************
 * Author     : xiaohao Z
 * Blog     : http://www.cnblogs.com/shu-xiaohao/
 * Last modified : 2014-05-13 22:22
 * Filename     : hdu_4582.cpp
 * Description     : 
 * ************************************************/

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#include <queue>
#include <stack>
#include <vector>
#include <set>
#include <map>
#define MP(a, b) make_pair(a, b)
#define PB(a) push_back(a)

using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
typedef pair<unsigned int,unsigned int> puu;
typedef pair<int, double> pid;
typedef pair<ll, int> pli;
typedef pair<int, ll> pil;

const int INF = 0x3f3f3f3f;
const double eps = 1E-6;
const int LEN = 2010;
vector<int> Map1[LEN], Map2[LEN];
int vex[LEN][1010];
int n, m, vis[LEN], ans;

void setvex(int pa, int pb){
    int tn = pb/30;
    vex[pa][tn] |= (1<<(pb%30));
}

bool getvex(int pa, int pb){
    int tn = pb/30;
    if(vex[pa][tn] & (1<<(pb%30))) return true;
    return false;
}

void dfs(int v, int fa){
    vis[v] = 1;
    memset(vex[v], 0, sizeof vex[v]);
    for(int i=0; i<Map2[v].size(); i++){
        int x = Map2[v][i];
        if(vis[x]) setvex(v, x);
    }
    for(int i=0; i<Map1[v].size(); i++){
        int x = Map1[v][i];
        if(x != fa) dfs(x, v);
    }
    if(fa != -1 && getvex(v, fa)) ans++;
    else if(fa != -1){
        for(int i=0; i<210; i++) vex[fa][i] |= vex[v][i];
    }
}

int main()
{
//    freopen("in.txt", "r", stdin);

    int a, b;
    while(scanf("%d%d", &n, &m)!=EOF){
        if(!n && !m) break;
        ans = 0;
        for(int i=0; i<LEN; i++) Map1[i].clear();
        for(int i=0; i<LEN; i++) Map2[i].clear();
        for(int i=0; i<m; i++){
            scanf("%d%d", &a, &b);
            if(i < n-1){
                Map1[a].PB(b);
                Map1[b].PB(a);
            }else{
                Map2[a].PB(b);
                Map2[b].PB(a);
            }
        }
        memset(vis, 0, sizeof vis);
        dfs(1, -1);
        printf("%d
", ans);
    }
    return 0;
}