Cod sursa(job #1136156)

Utilizator gabrielinelusGabriel-Robert Inelus gabrielinelus Data 8 martie 2014 20:46:04
Problema Componente biconexe Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.68 kb
#include <cstdio>
#include <algorithm>
#include <vector>
#include <memory.h>
#include <bitset>
#include <stack>

#define Nmax 100005

using namespace std;
int N,M;
vector<int> G[Nmax],sol[Nmax];
bitset<Nmax> used;
int niv[Nmax],low[Nmax];
stack<int> S;

void read()
{
    scanf("%d%d",&N,&M);
    int a,b;
    for(int i = 1; i <= M; ++i)
    {
        scanf("%d%d",&a,&b);
        G[a].push_back(b);
        G[b].push_back(a);
    }
}
int nrbic;
void bix(int nodc,int nodn)
{
    ++nrbic;
    while(S.top()!= nodn)
    {
        sol[nrbic].push_back(S.top());
        S.pop();
    }
    sol[nrbic].push_back(nodn);
    sol[nrbic].push_back(nodc);
    S.pop();
}

void DFS(int k,int daddy)
{
    used[k] = 1;
    niv[k] = niv[daddy] + 1;
    low[k] = niv[k];
    for(vector<int>::iterator it = G[k].begin(); it != G[k].end(); ++it)
    {
        if(!used[*it])
        {
            S.push(*it);
            DFS(*it,k);
            low[k] = min(low[k],low[*it]);
            if(niv[k] <= low[*it])
                bix(k,*it);
            continue;
        }
        low[k] = min(low[k],niv[*it]);
    }
}

void solve()
{
    for(int i = 1; i <= N; ++i)
        if(!used[i])
        {
            S.push(i);
            DFS(1,0);
            S.pop();
        }
}
void print()
{
    printf("%d\n",nrbic);
    for(int i = 1; i <= nrbic; ++i)
    {
        for(vector<int>::iterator it = sol[i].begin(); it != sol[i].end(); ++it)
            printf("%d ",*it);
        printf("\n");
    }
}

int main()
{
    freopen("biconex.in","r",stdin);
    freopen("biconex.out","w",stdout);

    read();
    solve();
    print();

    return 0;
}