POJ 2553 强联通分量

点击打开链接

题意：给一个有向图，然后现在定义一个特殊点，定义为它能够到达的点也都能够到达它，问有哪些点是特殊点

思路：有向图问的是相互到达，那么强联通分量先缩点就好了，对于每一个强联通来说，这个分量内的点都是相互到达的，但是若这个联通分量里有一个点与分量外的点相连，那么这个联通分量内的所有点都不满足特殊点，所以我们需要找到那些分量没有出度的

#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3fll;
const int maxn=10010;
vector<int> G[maxn];
vector<int> rG[maxn];
vector<int> vs;
bool used[maxn];
int cmp[maxn],V;
void add_edge(int from,int to){
    G[from].push_back(to);
    rG[to].push_back(from);
}
void dfs(int v){
    used[v]=1;
    for(int i=0;i<G[v].size();i++){
        if(!used[G[v][i]]) dfs(G[v][i]);
    }
    vs.push_back(v);
}
void rdfs(int v,int k){
    used[v]=1;
    cmp[v]=k;
    for(int i=0;i<rG[v].size();i++){
        if(!used[rG[v][i]]) rdfs(rG[v][i],k);
    }
}
int scc(){
    memset(used,0,sizeof(used));
    vs.clear();
    for(int v=0;v<V;v++) if(!used[v]) dfs(v);
    memset(used,0,sizeof(used));
    int k=0;
    for(int i=vs.size()-1;i>=0;i--){
        if(!used[vs[i]]) rdfs(vs[i],k++);
    }
    return k;
}
int A[maxn*10],B[maxn*10],vis[maxn],tmp[maxn];
int main(){
    int m;
    while(scanf("%d",&V)!=-1){
        if(V==0) break;
        scanf("%d",&m);
        for(int i=0;i<maxn;i++){
            G[i].clear();
            rG[i].clear();
        }
        vs.clear();
        for(int i=0;i<m;i++){
            scanf("%d%d",&A[i],&B[i]);
            A[i]--;B[i]--;
            add_edge(A[i],B[i]);
        }
        memset(vis,0,sizeof(vis));
        int ans=scc();
        for(int i=0;i<m;i++){
            int a=cmp[A[i]],b=cmp[B[i]];
            if(a!=b) vis[a]=1;
        }
        int sum=0;
        for(int i=0;i<V;i++) if(vis[cmp[i]]==0) tmp[sum++]=i+1;
        sort(tmp,tmp+sum);
        for(int i=0;i<sum;i++){
            printf("%d",tmp[i]);
            if(i!=sum-1) printf(" ");
        }
        printf("\n");
    }
    return 0;
}