vijos 1769 割边 ***

题意：考虑一个连通的无向图，可以知道，任意两个节点都可以通过一条路径连接起来。在所有节点中，某些节点向所有与它连通的节点提供A服务（包括向它自己），同时某些节点向所有与它连通的节点提供B服务（也包括向它自己）。注意一个节点也可能同时提供A、B两种服务。

当图中的某条边E被去掉的时候，如果图中有任何一个点无法接受A服务或者接受B服务，我们称E边为关键边。

那么，你需要做的事情就是：
1、输出图中存在多少关键边；
2、从小到大输出所有这样的关键边的编号。

链接：点我

首先缩点，缩点的过中找出割边。然后重新建一棵树，在树上跑一边dfs每棵子树的A的个数和B的个数。最后再跑一边dfs找那些子树中A[v] == 0 || B[v] == 0｜｜ A[v] == k || B[v] ==l

 1 #include <cstdio>
 2 #include <cstring>
 3 #include <cmath>
 4 #include <string>
 5 #include <iostream>
 6 #include <algorithm>
 7 #include <queue>
 8 using namespace std;
 9 #define rep(i, n) for(int i=0; i<(n); ++i)
10 #define for1(i,a,n) for(int i=(a);i<=(n);++i)
11 #define for2(i,a,n) for(int i=(a);i<(n);++i)
12 #define for3(i,a,n) for(int i=(a);i>=(n);--i)
13 #define for4(i,a,n) for(int i=(a);i>(n);--i)
14 #define CC(i,a) memset(i,a,sizeof(i))
15 #define read(a) a=getint()
16 #define print(a) printf("%d", a)
17 #define dbg(x) cout << (#x) << " = " << (x) << endl
18 #define printarr2(a, b, c) for1(_, 1, b) { for1(__, 1, c) cout << a[_][__]; cout << endl; }
19 #define printarr1(a, b) for1(_, 1, b) cout << a[_] << '\t'; cout << endl
20 inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<'0'||c>'9'; c=getchar()) if(c=='-') k=-1; for(; c>='0'&&c<='9'; c=getchar()) r=r*10+c-'0'; return k*r; }
21 inline const int max(const int &a, const int &b) { return a>b?a:b; }
22 inline const int min(const int &a, const int &b) { return a<b?a:b; }
23 
24 const int N=100005;
25 int ihead[N], cnt, n, m, l, k, FF[N], LL[N], ed[N], dep[N], tot, TM, X[N], Y[N], vis[N], a[N], b[N];
26 struct ED { int to, next, id; }e[N<<1];
27 void add(int u, int v, int id) {
28     e[++cnt].next=ihead[u]; ihead[u]=cnt; e[cnt].to=v; e[cnt].id=id;
29     e[++cnt].next=ihead[v]; ihead[v]=cnt; e[cnt].to=u; e[cnt].id=id;
30 }
31 void tarjan(int x, int fa) {
32     FF[x]=LL[x]=++TM;
33     int y;
34     for(int i=ihead[x]; i; i=e[i].next) {
35         y=e[i].to;
36         if(!FF[y]) {
37             tarjan(y, x);
38             if(LL[y]>FF[x]) ed[++tot]=e[i].id;
39             LL[x]=min(LL[x], LL[y]);
40         }
41         else if(FF[y]<FF[x] && y!=fa) LL[x]=min(LL[x], FF[y]);
42     }
43 }
44 void dfs(int x, int fa) {
45     dep[x]=dep[fa]+1;
46     vis[x]=1;
47     for(int i=ihead[x]; i; i=e[i].next) if(!vis[e[i].to]) {
48         dfs(e[i].to, x);
49         a[x]+=a[e[i].to];
50         b[x]+=b[e[i].to];
51     }
52 }
53 int main() {
54     #ifndef ONLINE_JUDGE
55     freopen("1.in","r",stdin);
56     #endif
57     read(n); read(m); read(k); read(l);
58     for1(i, 1, k) a[getint()]=1;
59     for1(i, 1, l) b[getint()]=1;
60     for1(i, 1, m) {
61         read(X[i]); read(Y[i]);
62         add(X[i], Y[i], i);
63     }
64     for1(i, 1, n) if(!FF[i]) tarjan(i, 0);
65     dfs((n+1)>>1, 0); CC(vis, 0);
66     sort(ed+1, ed+1+tot); ed[tot+1]=-1; int num=0;
67     /*for(int i=1;i<=tot;i++)
68     {
69         printf("%d ",ed[i]);
70     }
71     printf("\n");*/
72     for1(i, 1, tot) if(ed[i]!=ed[i+1]) ed[++num]=ed[i]; tot=num;
73     int ans=0;
74     for1(i, 1, tot) {
75         int x=X[ed[i]], y=Y[ed[i]];
76         //printf("x=%d  y=%d\n",X[ed[i]],Y[ed[i]]);
77         if(dep[y]<dep[x]) swap(x, y);
78         if(a[y]==0 || a[y]==k || b[y]==0 || b[y]==l) vis[i]=1, ++ans;   //子树中A和B的个数为0，说明到不了子树，即为所求割边
79     }
80     printf("%d\n", ans);
81     for1(i, 1, tot) if(vis[i]) printf("%d\n", ed[i]);
82     return 0;
83 }
2015/5/29

 

posted @ 2015-06-01 00:03  miao_a_miao  阅读(112)  评论(0)    收藏  举报