Cod sursa(job #1356277)

Utilizator gabrielinelusGabriel-Robert Inelus gabrielinelus Data 23 februarie 2015 12:27:04
Problema Componente tare conexe Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.45 kb
#include <algorithm>
#include <cstdio>
#include <vector>
#include <bitset>
#include <stack>

#define Nmax 100005

using namespace std;
int N,M;
stack<int>S;
bitset<Nmax> used;
vector<int> G[Nmax],Gt[Nmax],scc[Nmax];
int nrscc;

void Read()
{
    int a,b;
    scanf("%d%d",&N,&M);
    for(int i = 1; i <= M; ++i){
        scanf("%d%d",&a,&b);
        G[a].push_back(b);
        Gt[b].push_back(a);
    }
}

void DFS(int k)
{
    used[k] = 1;
    for(vector<int>::iterator it = G[k].begin(); it != G[k].end(); ++it)
        if(!used[*it])
            DFS(*it);
    S.push(k);
}

void DFST(int k)
{
    scc[nrscc].push_back(k);
    used[k] = 1;
    for(vector<int>::iterator it = Gt[k].begin(); it != Gt[k].end(); ++it)
        if(!used[*it])
            DFST(*it);
}

void Solve()
{
    for(int i = 1; i <= N; ++i)
        if(!used[i])
            DFS(i);
    used = 0;
    int k;
    while(!S.empty())
    {
        k = S.top();
        S.pop();
        if(used[k])
            continue;
        ++nrscc;
        DFST(k);
    }

}

void Print()
{
    printf("%d\n",nrscc);
    for(int i = 1; i <= nrscc; ++i)
    {
        for(vector<int>::iterator it = scc[i].begin(); it != scc[i].end(); ++it)
            printf("%d ",*it);
        printf("\n");
    }
}

int main()
{
    freopen("ctc.in","r",stdin);
    freopen("ctc.out","w" ,stdout);

    Read();
    Solve();
    Print();

    return 0;
}