Cod sursa(job #851107)

Utilizator Catah15Catalin Haidau Catah15 Data 9 ianuarie 2013 15:54:56
Problema Componente tare conexe Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.48 kb
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>

using namespace std;

#define maxN 100010
#define PB push_back

vector <int> list[maxN], listT[maxN], ans[maxN], Q;
bool Viz[maxN];
int N, sol;


void DFP (int nod)
{
    Viz[nod] = true;

    for (int i = 0; i < list[nod].size(); ++ i)
    {
        int nod2 = list[nod][i];
        if (Viz[nod2]) continue;
        DFP (nod2);
    }

    Q.PB (nod);
}


void DFM (int nod)
{
    Viz[nod] = true;
    ans[sol].PB (nod);

    for (int i = 0; i < listT[nod].size(); ++ i)
    {
        int nod2 = listT[nod][i];
        if (Viz[nod2]) continue;

        DFM (nod2);
    }
}


void Kosaraju ()
{
    for (int i = 1; i <= N; ++ i)
    {
        if (Viz[i]) continue;
        DFP (i);
    }

    memset (Viz, false, sizeof (Viz));

    for (int i = Q.size() - 1; i >= 0; -- i)
    {
        int nod = Q[i];
        if (Viz[nod]) continue;

        ++ sol;
        DFM (nod);
    }
}


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

    int M;
    scanf ("%d %d", &N, &M);

    while (M --)
    {
        int x, y;
        scanf ("%d %d", &x, &y);
        list[x].PB (y);
        listT[y].PB (x);
    }

    Kosaraju ();

    printf ("%d\n", sol);
    for (int i = 1; i <= sol; ++ i)
    {
        for (int j = 0; j < ans[i].size(); ++ j) printf ("%d ", ans[i][j]);
        printf ("\n");
    }

    return 0;
}