Cod sursa(job #265552)

Utilizator marcelcodreaCodrea Marcel marcelcodrea Data 24 februarie 2009 01:07:46
Problema Componente biconexe Scor 10
Compilator cpp Status done
Runda Arhiva educationala Marime 3.42 kb
#include<stdio.h>

struct Nod {
    int x;
    Nod *next;
};
Nod *a[100005];
int p[100005];
int viz[100005];
int stackx[100005];
int stacky[100005];
int nrBC;
int top;
int nrf;
int n,m,x,y;
int sn[100005];
int biconex[400001];
int art[100005];
int smin[100005];

int min(int a, int b)
{
    if (a > b) return b;
    return a;
}
int insert(Nod *&u, int val)
{
    Nod *q = new Nod;
    q -> x = val;
    q -> next = u;
    u = q;
    return 0;
}
void InitViz()
{
    for(int i = 1; i <= n; i++)
     viz[i] = 0;
}
int dfs(int s, int niv)
{
    sn[s] = niv;
    for(Nod *it = a[s]; it != NULL; it = it -> next)
     if (!sn[it -> x])
     {
      p[it -> x] = s;
      dfs(it -> x, niv + 1);
     }
    return 0;
}
int dfmin(int s)
{
    viz[s] = 1;
    int mini = smin[s];
    for(Nod *it = a[s]; it != NULL; it = it -> next)
     if ((it -> x) != p[s])
      if (mini > sn[it->x]) mini = sn[it->x];
    smin[s] = mini;
    for(Nod *it = a[s]; it != NULL; it = it -> next)
     if (!viz[it -> x])
      smin[s] = min(smin[s], dfmin(it -> x));

    return smin[s];
}

int df2conex(int s)
{
    viz[s] = 1;
    for(Nod *it = a[s]; it != NULL; it = it -> next)
     if (it -> x != p[s])
     {
        if (viz[it -> x] && sn[it -> x] <= sn[s])
         {
           stackx[++top] = s;
           stacky[top]= it -> x;
         }

        if (!viz[it->x])
         {
           stackx[++top] = s;
           stacky[top] = it -> x;
           df2conex(it -> x);
           if (smin[it -> x] >= sn[s])
            {

                nrBC++;
                InitViz();
                while (stackx[top] != s || stacky[top] != it->x)
                 {
                     if (!viz[stackx[top]])
                     {
                         viz[stackx[top]] = 1;
                         biconex[++biconex[0]] = stackx[top];
                     }
                     if (!viz[stacky[top]])
                     {
                         viz[stacky[top]] = 1;
                         biconex[++biconex[0]] = stacky[top];
                     }
                      top--;
                 }

                  if (!viz[stackx[top]])
                     {
                         viz[stackx[top]] = 1;
                         biconex[++biconex[0]] = stackx[top];
                     }
                     if (!viz[stacky[top]])
                     {
                         viz[stacky[top]] = 1;
                         biconex[++biconex[0]] = stacky[top];
                     }
                    top--;
                    biconex[0]++;
            }

            }
         }
        return 0;
     }

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

    scanf("%d %d", &n, &m);
    for(int i = 1; i <= m; i++)
    {
     scanf("%d %d", &x, &y);
     insert(a[x], y);
     insert(a[y], x);
    }
    // nivelul fiecarui nod
    dfs(1,1);
    // nivelul minim la care se poate ajunge prin intermediul subarborelui
    for(int i = 1; i <= n; i++)
     smin[i] = sn[i];
    InitViz();
    dfmin(1);



    // determina componente biconexe
    InitViz();
    df2conex(1);
    InitViz();
    printf("%d\n", nrBC);
    for(int i = 1; i <= biconex[0]; i++)
     {
         if (biconex[i])
           printf("%d ",biconex[i]);
          else
           printf("\n");
     }

    return 0;

}