【BZOJ 1051】 1051: [HAOI2006]受欢迎的牛 （SCC）

1051: [HAOI2006]受欢迎的牛

Description

　　每一头牛的愿望就是变成一头最受欢迎的牛。现在有N头牛，给你M对整数(A,B)，表示牛A认为牛B受欢迎。 这

种关系是具有传递性的，如果A认为B受欢迎，B认为C受欢迎，那么牛A也认为牛C受欢迎。你的任务是求出有多少头

牛被所有的牛认为是受欢迎的。

Input

　　第一行两个数N,M。 接下来M行，每行两个数A,B，意思是A认为B是受欢迎的（给出的信息有可能重复，即有可

能出现多个A,B）

Output

　　一个数，即有多少头牛被所有的牛认为是受欢迎的。

Sample Input

3 3
1 2
2 1
2 3

Sample Output

1

HINT

100%的数据N<=10000,M<=50000

Source

【分析】

　　SCC缩点，判断出度为0的强联通分量是否<=1，如果是，那么就输出这个SCC的大小，否则ans为0。

 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #include<stack>
 using namespace std;
 #define Maxn 10010
 #define Maxm 50010

 int mymin(int x,int y) {return x<y?x:y;}

 struct node
 {
     int x,y,next,o;
     bool p;
 }t[Maxm*];
 int first[Maxn],len=;

 void ins(int x,int y)
 {
     t[++len].x=x;t[len].y=y;
     t[len].next=first[x];first[x]=len;
     // if(len^1) t[len].o=len+1;
     // else t[len].o=len-1;
     t[len].p=;
 }

 stack<int > sta;

 int dfn[Maxn],low[Maxn],scc[Maxn],cl;
 int siz[Maxn],tot;
 bool inst[Maxn];
 void dfs(int x)
 {
     dfn[x]=low[x]=++tot;
     sta.push(x);inst[x]=;
     for(int i=first[x];i;i=t[i].next) if(t[i].p)
     {
         // t[i].p=t[t[i].o].p=0;
         t[i].p=;
         int y=t[i].y;
         if(!dfn[y])
         {
             dfs(y);
             low[x]=mymin(low[x],low[y]);
         }
         else if(inst[y]) low[x]=mymin(low[x],dfn[y]);
     }
     if(low[x]==dfn[x])
     {
         cl++;siz[cl]=;
         while(sta.top()!=x)
         {
             scc[sta.top()]=cl;
             inst[sta.top()]=;
             sta.pop();
             siz[cl]++;
         }
         scc[x]=cl;siz[cl]++;
         inst[x]=;
         sta.pop();
     }
 }

 int cd[Maxn];

 int main()
 {
     int n,m;
     scanf("%d%d",&n,&m);
     for(int i=;i<=m;i++)
     {
         int x,y;
         scanf("%d%d",&x,&y);
         ins(x,y);
     }
     memset(dfn,,sizeof(dfn));
     memset(inst,,sizeof(inst));
     // sta.clear();
     cl=;tot=;
     for(int i=;i<=n;i++) if(!dfn[i])
     {
         dfs(i);
     }
     memset(cd,,sizeof(cd));
     for(int i=;i<=m;i++)
     {
         if(scc[t[i].x]==scc[t[i].y]) continue;
         cd[scc[t[i].x]]++;
     }
     int nw=,id;
     for(int i=;i<=cl;i++) if(cd[i]==) nw++,id=i;
     if(nw<=) printf("%d\n",siz[id]);
     else printf("0\n");
     return ;
 }

之前的SCC模板似乎有问题，不应该用fa判断的（有向图啊醉。。）

2017-02-21 21:56:30