Pagini recente » Cod sursa (job #2329594) | Cod sursa (job #2652957) | Cod sursa (job #2688389) | Cod sursa (job #2028903) | Cod sursa (job #2925411)
//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;
}