【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*2];
int first[Maxn],len=0;

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=1;
}

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]=1;
    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=0;
        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]=0;
        while(sta.top()!=x)
        {
            scc[sta.top()]=cl;
            inst[sta.top()]=0;
            sta.pop();
            siz[cl]++;
        }
        scc[x]=cl;siz[cl]++;
        inst[x]=0;
        sta.pop();
    }
}

int cd[Maxn];

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

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

2017-02-21 21:56:30