Cod sursa(job #2490170)

Utilizator maria15Maria Dinca maria15 Data 9 noiembrie 2019 20:43:44
Problema Componente biconexe Scor 0
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 2.95 kb
#include <fstream>
#include <vector>
#include <algorithm>
#include <deque>
#include <bitset>

using namespace std;


ifstream fin("biconex.in");
ofstream fout("biconex.out");

int n, i, t[100001], f[100001], low[100001], m, x, y, componente, j, p;
bitset<100002> pun;
vector<int> v[100001], punctedearticulatie;
vector<pair<int, int> > crit;
deque<pair<int, int> > stiva;
deque<int> sol[100001];

void dfs(int nod, int niv){
    f[nod] = niv;
    /// low[nod] = cel mai mic nivel la care pot ajunge dintr-un fiu al lui nod pe alta cale decat trecand iar prin tata
    /// cu siguranta, din nod pot ajunge in el insusi, deci low[nod] este macar egal cu niv
    low[nod] = niv;
    int fii = 0;    /// numarul de componente conexe care pornesc din nod
    for(int i = 0;i<v[nod].size();i++){
        int vecin = v[nod][i];
        if(f[vecin] == 0){
            t[vecin] = nod;
            fii++;
            stiva.push_back(make_pair(nod, vecin));
            dfs(vecin, niv+1);
            low[nod] = min(low[nod], low[vecin]);
            if(low[vecin] >= niv){   /// daca din unul dintre fii nu se poate ajunge mai sus de nod
                if(niv != 1 && pun[nod] == 0)
                    punctedearticulatie.push_back(nod), pun[nod] = 1;     /// atunci nod este punct de articulatie
                /// daca este punct de articulatie, inseamna ca delimiteaza o componenta conexa;
                /// asa ca sterg toate muchiile adaugate in stiva pana la aceasta muchie
                /// ele formand componenta biconexa
                componente++;
                y = stiva.back().second;
                x = stiva.back().first;
                if(stiva.back().first == nod){
                    sol[componente].push_back(stiva.back().second);
                    sol[componente].push_back(stiva.back().first);
                    stiva.pop_back();
                }
                else{
                    do{
                        sol[componente].push_back(stiva.back().second);
                        stiva.pop_back();
                    }
                    while(!stiva.empty() && stiva.back().first != nod);
                    stiva.pop_back();
                    sol[componente].push_back(nod);
                    sol[componente].push_back(vecin);
                }
            }
            if(low[vecin] > niv)
                crit.push_back(make_pair(nod, vecin));
        }
        else
            if(t[nod] != vecin)
                low[nod] = min(low[nod], f[vecin]);
    }
    if(fii >= 2 && nod == 1 && pun[nod] == 0)
        punctedearticulatie.push_back(nod), pun[nod] = 1;

}

int main(){
    fin>>n>>m;
    for(i=1;i<=m;i++){
        fin>>x>>y;
        v[x].push_back(y);
        v[y].push_back(x);
    }
    dfs(1, 1);
    fout<<componente<<"\n";
    for(i=1;i<=componente;i++){
        for(j=0;j<sol[i].size();j++)
            fout<<sol[i][j]<<" ";
        fout<<"\n";
    }

    return 0;
}