Cod sursa(job #1121044)

Utilizator rares96cheseliRares Cheseli rares96cheseli Data 25 februarie 2014 11:17:27
Problema Componente biconexe Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.53 kb
#include <fstream>
#include <vector>
#include <algorithm>
#define pb push_back
using namespace std;
ifstream f("biconex.in");
ofstream g("biconex.out");

typedef pair< int, int > edge;
int N, M, x, y, dad[100002], low[100002], found[100002], order;
vector < int > G[100002];
vector < edge > st;
vector < vector < int > > BCC;

template < class T >
inline void EliminaDuplicate(vector < T > &V)
{
    sort(V.begin(), V.end());
    V.erase(unique(V.begin(), V.end()), V.end());
}

void getbcc(edge M)
{
    vector < int > newbcc; edge E;
    do
    {
        E=st.back(); st.pop_back();
        newbcc.pb(E.first); newbcc.pb(E.second);
    }while (E!=M);
    EliminaDuplicate(newbcc);
    BCC.pb(newbcc);
}

void dfs(int nod)
{
    low[nod]=found[nod]=++order;
    vector < int >:: iterator it=G[nod].begin();
    for (; it!=G[nod].end(); ++it)
        if (!dad[*it])
        {
            dad[*it]=nod; st.pb(edge(nod, *it));
            dfs(*it);
            low[nod]=min(low[nod], low[*it]);
            if (low[*it]>=found[nod])
                getbcc(edge(nod, *it));
        }
        else if (*it!=dad[nod])
                low[nod]=min(low[nod], found[*it]);
}

int main()
{
    f>>N>>M;
    for (int i=1; i<=M; ++i)
        f>>x>>y, G[x].pb(y), G[y].pb(x);
    for (int i=1; i<=N; ++i)
        if (!dad[i]) dad[i]=i, dfs(i);
    g<<BCC.size()<<'\n';
    for (int i=0; i<(int)BCC.size(); ++i && g<<'\n')
        for (int j=0; j<(int)BCC[i].size(); ++j)
            g<<BCC[i][j]<<' ';
    return 0;
}