Cod sursa(job #1145092)

Utilizator uagamagaMatei Rogoz uagamaga Data 17 martie 2014 21:14:21
Problema Componente biconexe Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.8 kb
#include <fstream>
#include <vector>
#include <stack>
#include <set>
#define min(a,b) ( ((a)<(b))?(a):(b) )
#define DIM 100005
using namespace std;

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

int m,n;
vector <int> G[DIM];
set <int> compBiconexe[DIM];
stack < pair<int,int> > mStack;
int h[DIM],returnNode[DIM],timpDescoperire;
int cbx;


void dfsPct(int u,int parinte)
{
    h[u] = timpDescoperire;
    timpDescoperire++;
    returnNode[u] = h[u];

    for(int i=0;i<G[u].size();++i)
    {
        int v = G[u][i];

        if(!h[v])
        {
            mStack.push(make_pair(u,v));
            dfsPct(v,u);
            returnNode[u] = min(returnNode[u],returnNode[v]);

            if(returnNode[v]>=h[u])
            {
                cbx++;
                pair <int,int> node;
                do
                {
                    node = mStack.top();
                    mStack.pop();
                    compBiconexe[cbx].insert(node.first);
                    compBiconexe[cbx].insert(node.second);

                }while((node.first!=u || node.second!=v) && (node.second!=u || node.first!=v));
            }
        }
        else if(v!=parinte) returnNode[u] = min(returnNode[u],h[v]);
    }

}
void comp_biconexe()
{
    for(int u=1;u<=n;++u)
        if(!h[u])
        {
            timpDescoperire = 1;
            dfsPct(u,0);
        }
}

int main()
{
    int a,b;
    fin>>n>>m;

    while(m--)
    {
        fin>>a>>b;
        G[a].push_back(b);
        G[b].push_back(a);
    }

    comp_biconexe();

    fout<<cbx<<"\n";

    for(int i=1;i<=cbx;++i){
        for(set<int>::iterator j=compBiconexe[i].begin();j!=compBiconexe[i].end();++j){
            fout<<*j<<" ";
        }
        fout<<"\n";
    }

    return 0;
}