您的位置:首页 > 其它

UVA 1108 - Mining Your Own Business

2016-04-01 18:54 881 查看
刘汝佳书上都给出了完整的代码

在这里理一下思路:

由题意知肯定存在一个或者多个双连通分量;

假设某一个双连通分量有割顶。那太平井一定不能打在割顶上。

而是选择割顶之外的随意一个点;

假设没有割顶,则要在该双连通分量上打两个井

至于打井方案。见代码

#include <cstdio>
#include <cstring>
#include <vector>
#include <stack>
#include <map>
using namespace std;

const int N = 50005;

struct Edge {
int u, v;
Edge() {}
Edge(int u, int v) {
this->u = u;
this->v = v;
}
};

int pre
, bccno
, dfs_clock, bcc_cnt;
bool iscut
;

vector<int> g
, bcc
;
stack<Edge> S;

int dfs_bcc(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_bcc(v, u);
lowu = min(lowu, lowv);
if (lowv >= pre[u]) {
iscut[u] = true;
bcc_cnt++; bcc[bcc_cnt].clear(); //start from 1
while(1) {
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] = false;
return lowu;
}

int st;

void find_bcc() {
memset(pre, 0, sizeof(pre));
memset(iscut, 0, sizeof(iscut));
memset(bccno, 0, sizeof(bccno));
dfs_clock = bcc_cnt = 0;
dfs_bcc(0, -1);
}

int n, m;

typedef long long ll;

void solve() {
ll ans1 = 0, ans2 = 1;
for (int i = 1; i <= bcc_cnt; i++) {
int cut_cnt = 0;
for (int j = 0; j < bcc[i].size(); j++)
if (iscut[bcc[i][j]]) cut_cnt++;
if (cut_cnt == 1) {
ans1++;
ans2 *= (ll)(bcc[i].size() - cut_cnt);
}
}
if (bcc_cnt == 1) {
ans1 = 2;
ans2 = (ll)bcc[1].size() * (bcc[1].size() - 1) / 2;
}
printf(" %lld %lld\n", ans1, ans2);
}

int main() {
int cas = 0;
while (~scanf("%d", &m) && m) {
int u, v, Max = 0;
while (m--) {
scanf("%d%d", &u, &v);
u--; v--;
g[u].push_back(v);
g[v].push_back(u);
Max = max(Max, u);
Max = max(Max, v);
}
find_bcc();
printf("Case %d:", ++cas);
solve();
for (int i = 0; i <= Max; i++)
g[i].clear();
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: