Cod sursa(job #606920)

Utilizator a_h1926Heidelbacher Andrei a_h1926 Data 10 august 2011 13:41:13
Problema Componente tare conexe Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.84 kb
#include <iostream>
#include <vector>
#include <stack>
#include <cstdio>

#define NMax 100005

using namespace std;

vector <int> G[NMax];
int N, Index[NMax], CIndex, MinIndex[NMax];
stack <int> Stack;
bool InStack[NMax];
vector < vector <int> > CTC;

void Read ()
{
    freopen ("ctc.in", "r", stdin);
    int M;
    scanf ("%d %d", &N, &M);
    for (; M>0; --M)
    {
        int X, Y;
        scanf ("%d %d", &X, &Y);
        G[X].push_back (Y);
    }
}

void Print ()
{
    freopen ("ctc.out", "w", stdout);
    printf ("%d\n", CTC.size ());
    for (int i=0; i<CTC.size (); ++i)
    {
        for (int j=0; j<CTC[i].size (); ++j)
        {
            printf ("%d ", CTC[i][j]);
        }
        printf ("\n");
    }
}

inline int Min (int a, int b)
{
    if (a<b)
    {
        return a;
    }
    return b;
}

void DetCTC (int X)
{
    vector <int> Component;
    int Y;
    do
    {
        Y=Stack.top ();
        Stack.pop ();
        InStack[Y]=false;
        Component.push_back (Y);
    }
    while (Y!=X);
    CTC.push_back (Component);
}

void DFS (int X)
{
    Index[X]=MinIndex[X]=CIndex;
    Stack.push (X);
    InStack[X]=true;
    for (vector <int> :: iterator V=G[X].begin (); V!=G[X].end (); ++V)
    {
        if (Index[*V]==0)
        {
            ++CIndex;
            DFS (*V);
            MinIndex[X]=Min (MinIndex[X], MinIndex[*V]);
        }
        else
        {
            if (InStack[*V])
            {
                MinIndex[X]=Min (MinIndex[X], MinIndex[*V]);
            }
        }
    }
    if (MinIndex[X]==Index[X])
    {
        DetCTC (X);
    }
}

int main()
{
    Read ();
    for (int i=1; i<=N; ++i)
    {
        if (Index[i]==0)
        {
            ++CIndex;
            DFS (i);
        }
    }
    Print ();
    return 0;
}