Cod sursa(job #600529)

Utilizator nervousNervous nervous Data 2 iulie 2011 10:39:13
Problema Componente biconexe Scor 52
Compilator cpp Status done
Runda Arhiva educationala Marime 2.27 kb
#include <fstream>
#include <cstring>
#include <vector>
#include <stack>
#include <algorithm>

#define X1 100001

using namespace std;

ifstream in;
ofstream out;

vector <int> v[X1];
vector <int> biconex[X1];
stack < pair<int,int> > st;
int use[X1],T[X1],low[X1],level[X1],er[X1];
int cnt=0;

inline void df(int nod)
{
    int x,y;

    use[nod]=1;
    for(vector <int>::iterator it=v[nod].begin();it!=v[nod].end();++it)
    {
        if(T[nod]!=*it&&level[nod]>level[*it]) st.push(make_pair(nod,*it));
        if(!use[*it])
        {
            T[*it]=nod;
            level[*it]=level[nod]+1;
            low[*it]=level[*it];
            df(*it);
            if(low[nod]>low[*it]) low[nod]=low[*it];
            if(low[*it]>=level[nod])
            {
                ++cnt;
                memset(er,0,sizeof(er));

                if(st.size())
                    do
                    {
                        x=st.top().first;
                        y=st.top().second;
                        if(!er[x])
                        {
                            biconex[cnt].push_back(x);
                            er[x]=1;
                        }
                        if(!er[y])
                        {
                            biconex[cnt].push_back(y);
                            er[y]=1;
                        }
                        st.pop();
                    }
                    while((x!=*it||y!=nod)&&(x!=nod||y!=*it)&&st.size());
            }
        }
        else
        if(T[nod]!=*it)
            if(low[nod]>level[*it]) low[nod]=level[*it];
    }
}

int main()
{
    int M,N,x,y;

    in.open("biconex.in");
    in>>N>>M;
    for(;M;--M)
    {
        in>>x>>y;
        v[x].push_back(y);
        v[y].push_back(x);
    }
    in.close();

    memset(use,0,sizeof(use));
    memset(level,0,sizeof(level));
    memset(low,0,sizeof(low));
    memset(T,0,sizeof(T));

    for(int i=1;i<=N;++i)
        if(!use[i]) df(i);

    out.open("biconex.out");
    out<<cnt<<'\n';
    for(int nod=1;nod<=cnt;++nod)
    {
        for(vector <int>::iterator it=biconex[nod].begin();it!=biconex[nod].end();++it)
            out<<*it<<' ';
        out<<'\n';
    }
    out.close();

    return 0;
}