Poj 2186 强联通水题

题意：有向图中  有多少个点可以被其他所有点为起点到达

解法：求出强联通分量缩点以后 利用拓扑次序累加费用 最后检验每个点的费用是非是点数之和

如果是的话 说明这个点可以被所有点到达

#pragma comment(linker,"/STACK:102400000,102400000")
#include <iostream>
#include <string.h>
#include <stdio.h>
#include <algorithm>
#include <vector>
#include <string>
#include <math.h>
#include <queue>
#include <stack>
#include <map>
#include <set>
using namespace std;
typedef long long ll;   //记得必要的时候改成无符号
const int maxn=100005;   //点数
const int maxm=500005;   //边数
const int INF=1000000000;
struct EdgeNode
{
    int from;
    int to;
    int next;
}edge[maxm];
int head[maxn],cnt;
void add(int x,int y)
{
    edge[cnt].from=x;edge[cnt].to=y;edge[cnt].next=head[x];head[x]=cnt++;
}

void init()
{
    cnt=1;
    memset(head,-1,sizeof(head));
}

int dfn[maxn],low[maxn],ins[maxn];
int cixu,scc_count,n,ru[maxn],chu[maxn],sd[maxn];
stack<int>S;

void dfs(int u)
{
    int v,x;
    dfn[u]=low[u]=++cixu;
    S.push(u);
    ins[u]=1;
    for(int i=head[u];i!=-1;i=edge[i].next)
    {
        v=edge[i].to;
        if(!dfn[v]){
            dfs(v);
            low[u]=min(low[u],low[v]);
        }
        else if(ins[v])
            //如果已经搜过，并且在栈里面，如果不在栈里面，说明该点属于其他的强连通分量里面
            low[u]=min(low[u],dfn[v]);
             //此处low代替dfn也是可以的，但是在tarjan求无向图的时候不能
    }
    if(dfn[u]==low[u])
//表示u为此强连通分量第一个被发现的点，属于该强连通分量的所有的其他点都在栈中，他的上面
    {
        scc_count++;
        do
        {
            x=S.top();
            S.pop();
            sd[x]=scc_count;
            ins[x]=0;
        }while(x!=u);
    }
}

void tarjan()
{
    for(int i=1;i<=n;i++)sd[i]=i;   //缩点标记数组
    memset(ins,0,sizeof(ins));
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    cixu=scc_count=0;
    for(int i=1;i<=n;i++)        //可能所给的图为不连通的图，分块，因此需要遍历所有的点
    {
        if(!dfn[i])
            dfs(i);
    }
}

void SD()
{
    int x,y,len=cnt;
    init();
    memset(ru,0,sizeof(ru));
    memset(chu,0,sizeof(chu));
    for(int i=1;i<len;i++){
        x=sd[edge[i].from]; y=sd[edge[i].to];
        if(x==y)continue;
        add(x,y);
        chu[x]++; ru[y]++;
    }
}
int value[maxn],va[maxn];
// scc_count 1...
int main()
{
    int m,i,p,sum,x,y;
    while(~scanf("%d%d",&n,&m)){
        init(); sum=0; p=0;
        for(i=1;i<=m;i++){
            scanf("%d%d",&x,&y);
            add(x,y);
        }
        tarjan();

        
//        printf("%d\n",scc_count);
        for(int i=1;i<=scc_count;++i)value[i]=0;
        SD();
        for(int i=1;i<=n;++i)value[sd[i]]++;
        
        priority_queue<int>q;
        for(int i=1;i<=scc_count;++i)
            if(ru[i]==0)
                q.push(i);
        while(!q.empty()){
            int x=q.top();q.pop();
            for(int i=head[x];i!=-1;i=edge[i].next){
                int v=edge[i].to;
                value[v]+=value[x];
                ru[v]--;
                if(ru[v]==0)q.push(v);
            }
        }
        
        for(int i=1;i<=n;++i){
            if(value[sd[i]]==n){
                sum++;
            }
        }
        
        printf("%d\n",sum);
        
    }
    return 0;
}