hdu 3844 Mining Your Own Business

见刘汝佳训练指南P318
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <iostream>
#include <stack>
#include <queue>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <set>
#include <vector>
#include <cstring>
#include <algorithm>

#define INF 0x3fffffff
#define N 50010
#define M 1000010
#define LL long long
#define mod 95041567

using namespace std;

struct Node{
    int u, v;
    Node(int x = 0, int y = 0){ 
        u = x;
        v = y;
    }
};

Node G[N];
int dfs_clock, bcc_cnt, G_len;
int head[N], next[N * 2][2], pre[N], bccno[N];
bool iscut[N];
vector<int> bcc[N];

void init(int n){
    for(int i = 0; i <= n + 1; ++ i) 
        head[i] = -1, iscut[i] = 0, pre[i] = bccno[i] = 0;
    dfs_clock = bcc_cnt = G_len = 0;
}

void add(int u, int v){
    next[dfs_clock][1] = v;
    next[dfs_clock][0] = head[u];
    head[u] = dfs_clock ++;
}

int dfs(int u, int fa){
    int child = 0;
    pre[u] = ++ dfs_clock;
    int lowu = pre[u];
    for(int i = head[u]; i != -1; i = next[i][0]){
        int v = next[i][1];
        if(! pre[v]){
            G[G_len ++] = Node(u, v);
            ++ child;
            int lowv = dfs(v, u);
            lowu = min(lowu, lowv);
            if(lowv >= pre[u]){
                iscut[u] = true;
                bcc[++ bcc_cnt].clear();
                while(1){
                    Node k = G[-- G_len];
                    if(bccno[k.u] != bcc_cnt){
                        bcc[bcc_cnt].push_back(k.u);
                        bccno[k.u] = bcc_cnt;
                    }
                    if(bccno[k.v] != bcc_cnt){
                        bcc[bcc_cnt].push_back(k.v);
                        bccno[k.v] = bcc_cnt;
                    }
                    if(k.u == u && k.v == v) break;
                }
            }
        }
        else if(pre[u] > pre[v] && v != fa){
            G[G_len ++] = Node(u, v);
            lowu = min(lowu, pre[v]);
        }
    }
    if(fa < 0 && child == 1) iscut[u] = 0;
    return lowu;
}

int main()
{
  //  freopen("in.txt","r",stdin);
    int n, t = 0;
    while(scanf("%d", &n) != EOF){
        if(! n) break;
        init(n);
        int u, v;
        for(int i = 0; i < n; ++ i){
            scanf("%d %d", &u, &v);
            -- u, -- v;
            add(u, v);
            add(v, u);
        }
        dfs_clock = 0;
        dfs(0, -1);
        if(bcc_cnt == 1){
            LL c = bcc[1].size();
            printf("Case %d: 2 %I64d\n", ++t, c * (c - 1) / 2);
            continue;
        }
        v = 0;
        LL c = 1;
        for(int i = 1; i <= bcc_cnt; ++ i){
            int cnt = 0;
            for(int j = bcc[i].size() - 1; j >= 0 && cnt < 2; -- j)
                if(iscut[bcc[i][j]]) ++ cnt;
            if(cnt == 1) ++ v, c *= (LL)(bcc[i].size() - 1);
        }
        printf("Case %d: %d %I64d\n", ++t, v, c);
    }
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值