Cod sursa(job #1555563)

Utilizator Liviu98Dinca Liviu Liviu98 Data 23 decembrie 2015 09:19:00
Problema Componente tare conexe Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.49 kb
#include <iostream>
#include <stdio.h>
#include <stack>
#include <vector>
#define NMax 100005
using namespace std;
vector<int> Graf[NMax];
vector<int> GrafT[NMax],Solution[NMax];
int Comp[NMax];
int NumarComp,N,M,x,y;
stack<int> Stiva;
bool mark[NMax];

void Read()
{
    scanf("%d%d",&N,&M);
    for(int i=1;i<=M;i++)
    {
        scanf("%d%d",&x,&y);
        Graf[x].push_back(y);
        GrafT[y].push_back(x);
    }
}

void DFS1(int nod)
{
    mark[nod]=true;
    for(vector<int>::iterator it=Graf[nod].begin();it!=Graf[nod].end();it++)
        if(mark[*it]==false)
          DFS1(*it);
    Stiva.push(nod);
}

void DFS2(int nod)
{
    Comp[nod]=NumarComp;
    Solution[NumarComp].push_back(nod);
    for(vector<int>::iterator it=GrafT[nod].begin();it!=GrafT[nod].end();it++)
        if(Comp[*it]==0)
         DFS2(*it);
}

void Kosaraju()
{
    for(int i=1;i<=N;i++)
        if(mark[i]==false)
        DFS1(i);

    while(Stiva.empty()==false)
    {
        if(Comp[Stiva.top()]==0)
        {
            NumarComp++;
            DFS2(Stiva.top());
        }
        Stiva.pop();
    }
}

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

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