Cod sursa(job #2925411)

Utilizator DragosG12Ghinea Dragos DragosG12 Data 15 octombrie 2022 10:42:18
Problema Componente tare conexe Scor 10
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 2.18 kb
//complexitate timp O(n+m)
//complexitate memorie O(n+m)

#include <iostream>
#include<fstream>
#include<vector>
#include<stack>
using namespace std;

vector<int> muchii[100002];
stack<int> istoric;
int radacini[100002]={};
int pozitii[100002]={};
int vizitat[100002]={};
//index folosit pentru ordinea vizitarii nodurilor
int poz=1;
vector<vector<int>> componente;

//algoritmul lui Torjan
void cauta(int nod){
    //initial nodu are accesibil doar pe el insusi
    pozitii[nod]=poz;
    radacini[nod]=poz;
    poz++;
    istoric.push(nod);
    vizitat[nod]=1;

    //parcurgere in adancime, trecand prin vecini dar apeland recursiv functia asta
    for(int vecin : muchii[nod]){
        //nu are cea mai mica pozitie accesibila setata, deci putem intra pe el
        if(radacini[vecin]==0){
            cauta(vecin);
            //actualizam cea mai mica pozitie accesibila atunci cand revine din recursie
            radacini[nod] = min(radacini[nod], radacini[vecin]);
        }
        //daca este vizitat, inseamna ca este in istoric (resetam vizitarile dupa fiecare componenta conexa), incercam sa actualizam cea mai mica pozitie accesibila
        else if(vizitat[nod]==1){
            radacini[nod] = min(radacini[nod], pozitii[vecin]);
        }
    }

    //cand ajunge aici toti vecinii au fost vizitati, verificam daca din nodu asta incepe componenta conexa, si o construim cu istoricul
    if(pozitii[nod]==radacini[nod]){
        vector<int> component;
        int vecin;
        do{
            vecin=istoric.top(); istoric.pop();
            vizitat[vecin]=0;
            component.push_back(vecin);

        }while(vecin!=nod);
        componente.push_back(component);
    }
}

int main()
{
    int n,m;
    ifstream fin("ctc.in");
    fin>>n>>m;
    int x,y;
    for(int i=0;i<m;i++){
        fin>>x>>y;
        muchii[x].push_back(y);
    }

    for(int i=1;i<=n;i++)
        if(radacini[i]==0)
            cauta(i);

    ofstream fout("ctc.out");
    fout<<componente.size()<<endl;
    for(vector<int>& component : componente){
        for(int nod : component)
            fout<<nod<<" ";
        fout<<endl;
    }

    return 0;
}