Cod sursa(job #2793556)

Utilizator PopescuMihneaPopescu Mihnea-Valentin PopescuMihnea Data 3 noiembrie 2021 18:52:57
Problema Componente biconexe Scor 90
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 6.68 kb
#include <iostream>
#include <fstream>
#include <vector>
#include <stack>
#include <algorithm>
#define PROBLEMA 3
using namespace std;
ifstream f;
ofstream g;
class graf
{
    int n,m,s;
    vector< vector<unsigned int> > lista_vecini;
public:
    graf(unsigned int);
    int setn(const int x);
    int setm(const int x);
    int getn();
    int getm();
    int gets();
    void getvecini();
    void sortvecini();
    void BFS();
    void DFS(unsigned int);
    void biconex(unsigned int);
};


graf::graf(unsigned int opt)
{

    switch (opt)
    {
    case 1:
    {
        int x,y;
        f>>this->n>>this->m>>this->s;
        this->lista_vecini.resize(n+1);
        for (int i=0; i<m; ++i)
        {
            f>>x>>y;
            this->lista_vecini[x].push_back(y);
        }
        break;
    }
    case 2:
    {
        int x,y;
        f>>this->n>>this->m;
        this->lista_vecini.resize(n+1);
        for (int i=0; i<m; ++i)
        {
            f>>x>>y;
            this->lista_vecini[x].push_back(y);
            this->lista_vecini[y].push_back(x);
        }
        break;
    }
    default:
        break;
    }
}

int graf::getm()
{
    return this->m;
}

int graf::getn()
{
    return this->n;
}

int graf::gets()
{
    return this->s;
}

void graf::getvecini()
{
    for (unsigned int i=1; i<this->lista_vecini.size(); ++i)
    {
        cout<<i<<": ";
        for (unsigned int j=0; j<this->lista_vecini[i].size(); ++j)
            cout<<this->lista_vecini[i][j]<<" ";
        cout<<"\n";
    }
}

void graf::sortvecini()
{
    int j,k;
    unsigned int ap[this->n+1];
    for (j=0; j<=this->n; ++j)
        ap[j]=0;
    for (int i=1; i<this->lista_vecini.size(); ++i)
    {
        k=0;
        for (j=0; j<lista_vecini[i].size(); ++j)
            ++ap[lista_vecini[i][j]];
        for (j=1; j<=this->n; ++j)
        {
            if (ap[j])
            {
                lista_vecini[i][k++]=j;
                ap[j]=0;
            }
        }
    }
}

void graf::BFS()
{
    int d[this->n+1];
    for (int i=1; i<=n; ++i)
        d[i]=-1;
    bool viz[this->n+1];
    for (int i=1; i<=n; ++i)
        viz[i]=false;
    unsigned int index=0;
    vector <int> coada;
    coada.push_back(this->s);
    d[this->s]=0;
    viz[this->s]=true;
    while (index<coada.size())
    {
        int nod_curent=coada[index++];
        for (unsigned int i=0; i<this->lista_vecini[nod_curent].size(); ++i)
            if (!viz[this->lista_vecini[nod_curent][i]])
            {
                coada.push_back(this->lista_vecini[nod_curent][i]);
                viz[this->lista_vecini[nod_curent][i]]=true;
                d[this->lista_vecini[nod_curent][i]]=d[nod_curent]+1;
            }
    }
    for (int i=1; i<=n; ++i)
    {
        g<<d[i]<<" ";
    }
}

void graf::DFS(unsigned int nod_curent)
{
    static bool *viz=new bool[this->n+1];
    static unsigned int nrconex;
    if (nod_curent==0)
    {
        for (int i=1; i<=n; ++i)
            viz[i]=false;
        nrconex=0;
        for (int i=1; i<=n; ++i)
            if (!viz[i])
            {
                nrconex++;
                DFS(i);
            }
        g<<nrconex;
        delete [] viz;
    }
    else
    {
        viz[nod_curent]=true;
        for (unsigned int i=0; i<this->lista_vecini[nod_curent].size(); ++i)
            if (!viz[lista_vecini[nod_curent][i]])
                DFS(lista_vecini[nod_curent][i]);
    }
}

void graf::biconex(unsigned int nod_curent)
{
    static int *niv=new int[this->n+1];
    static int *niv_int=new int[this->n+1];
    static unsigned int nrbiconex=0;
    static stack <pair<unsigned int,unsigned int>> stiva;
    static vector <vector<unsigned int>> afisare;
    unsigned int i;
    if (nod_curent==0)
    {
        for (i=1; i<=n; ++i)
        {
            niv[i]=-1;
            niv_int[i]=0;
        }
        niv[0]=0;
        biconex(1);
        g<<nrbiconex<<"\n";
        for (i=0;i<afisare.size();++i)
        {
            sort(afisare[i].begin(), afisare[i].end());
        afisare[i].erase(unique(afisare[i].begin(), afisare[i].end()), afisare[i].end());
        for (unsigned int j = 0; j < afisare[i].size(); ++ j)
             g<<afisare[i][j]<<" ";
             g<<"\n";
        }
        delete [] niv;
        delete [] niv_int;
    }
    else
    {
        for (unsigned int i=0; i<this->lista_vecini[nod_curent].size(); ++i)
            if (niv[lista_vecini[nod_curent][i]]==-1)
            {
                niv[lista_vecini[nod_curent][i]]=niv[nod_curent]+1;
                niv_int[lista_vecini[nod_curent][i]]=niv[nod_curent]+1;
                stiva.push(make_pair(nod_curent,lista_vecini[nod_curent][i]));
                biconex(lista_vecini[nod_curent][i]);
                niv_int[nod_curent]=min(niv_int[nod_curent],niv_int[lista_vecini[nod_curent][i]]);
                if (niv_int[lista_vecini[nod_curent][i]]>=niv[nod_curent])
                {
                    nrbiconex++;
                    vector<unsigned int> componente;
                    bool first=true;
                    unsigned int aux1;
                    unsigned int aux2;
                    do
                    {

                        aux1=stiva.top().first;
                        aux2=stiva.top().second;
                        componente.push_back(aux1);
                        componente.push_back(aux2);
                        stiva.pop();
                    }
                    while(aux1!=nod_curent || aux2!=lista_vecini[nod_curent][i]);
                    afisare.push_back(componente);
                }
            }
            else if (niv[nod_curent]-1!=niv[lista_vecini[nod_curent][i]])
                niv_int[nod_curent]=min(niv_int[nod_curent],niv[lista_vecini[nod_curent][i]]);
        //scadem 1 ca sa nu luam in considerare si tatal nodului
    }
}

int main()
{
    switch (PROBLEMA)
    {
    case 1:
    {
        f.open ("bfs.in", std::ifstream::in);
        g.open ("bfs.out", std::ifstream::out);
        graf a(1);
        a.sortvecini();
        a.BFS();
        f.close();
        g.close();
        break;
    }
    case 2:
    {
        f.open ("dfs.in", std::ifstream::in);
        g.open ("dfs.out", std::ifstream::out);
        graf a(2);
        a.sortvecini();
        a.DFS(0);
        f.close();
        g.close();
        break;
    }
    case 3:
    {
        f.open("biconex.in",std::fstream::in);
        g.open("biconex.out",std::fstream::out);
        graf a(2);
        a.biconex(0);
        f.close();
        g.close();
    }
    default:
        break;
    }
}