[BZOJ1051][HAOI2006]受欢迎的牛 做题笔记

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1051
可以想到一个朴素的暴力思路，那就是如果一头奶牛符合条件，那么它一定可以被所有其它奶牛结点访问到，所以可以用N^2的复杂度来枚举每一头牛所能到达的结点，同时记录每一个结点被访问了多少次（包括自己访问自己的一次）。最后枚举每一个结点，被访问的次数等于n的就符合条件。
可以用强联通缩点，新图上跑暴力。
注意图是有向的，邻接表不要写反向边，这是我被自己坑的第二次了。。。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <stack>
#include <vector>
using namespace std;
const int N=10009,M=60009,inf=0x3fffffff;
int head[N],e[M],ver[M],next[M];
int dfn[N],low[N],q[N<<1],hav[N],bel[N],gue[N];
bool inq[N],vis[N];
int tot=1,t=0,n,m,scc=0,K,ans=0,top=0;
vector<int> G[N];

void add (int u,int v) {
    ver[++tot]=v;next[tot]=head[u];head[u]=tot;
}

void Tarjan (int x) {
    low[x]=dfn[x]=++t;
    q[++top]=x;inq[x]=1;
    for (int i=head[x];i;i=next[i]) 
        if (!dfn[ver[i]]) {
            Tarjan(ver[i]);
            low[x]=min(low[x],low[ver[i]]);
        }
        else if (inq[ver[i]]) 
            low[x]=min(low[x],dfn[ver[i]]);
    if (dfn[x]==low[x]) {
        int now=0;++scc;
        while (x!=now&&top) {
            now=q[top--];
            bel[now]=scc;
            ++hav[bel[now]];
            inq[now]=0;
        }
    }
}

void dfs (int x) {
    vis[x]=1;gue[x]+=hav[K];
    for (int i=0;i<G[x].size();i++) 
        if (!vis[G[x][i]]) {
            dfs(G[x][i]);
        }
}

int main () {
    int u,v;
    memset(head,0,sizeof(head));
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(hav,0,sizeof(hav));
    memset(bel,0,sizeof(bel));
    memset(gue,0,sizeof(gue));
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++) {
        scanf("%d%d",&u,&v);
        add(u,v);
    }
    for (int i=1;i<=n;i++) 
        if (!dfn[i]) Tarjan(i);
    for (int x=1;x<=n;x++) 
        for (int i=head[x];i;i=next[i]) 
            if (bel[ver[i]]!=bel[x]) 
                G[bel[x]].push_back(bel[ver[i]]);
    for (K=1;K<=scc;K++) {
        memset(vis,0,sizeof(vis));
        dfs(K);
    }
    for (int i=1;i<=n;i++) 
        if (gue[bel[i]]==n) ans++;
    printf("%d",ans);
}