Cod sursa(job #2333178)

Utilizator Valentin0709Datcu George Valentin Valentin0709 Data 31 ianuarie 2019 19:09:44
Problema Obiective Scor 15
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 2.82 kb
#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <set>
#include <iostream>
using namespace std;

#define NMAX 32005
#define foreach BOOST_FOREACH

FILE * f = fopen("obiective.in", "r");
FILE * g = fopen("obiective.out", "w");

int n, m, t, x, y, indexx, ssize, nrctc, ctc[NMAX], viz[NMAX], cost[NMAX], s[NMAX], idx[NMAX], lowlink[NMAX], onstack[NMAX];

set<int> l[NMAX], ls[NMAX], sol[NMAX];

int exista_muchie(int x, int y, set<int> a[]) {
    if(a[x].count(y)) return 1;
    else return 0;
}

void scrie(set<int> a) {
    for(set<int>::iterator it=a.begin(); it!=a.end(); ++it)
        cout << *it << '\n';
}

void citire_graf() {
    int x, y;

    fscanf(f, "%d%d", &n, &m);
    for(int i = 1; i <= m; i++) {
        fscanf(f, "%d%d", &x, &y);
        l[x].insert(y);
    }
}

void tarjan(int v) {
    int x;

    s[++ssize] = v; onstack[v] = 1;
    idx[v] = lowlink[v] = ++indexx;

    for(set<int>::iterator it = l[v].begin(); it != l[v].end(); ++it) {
        int u = *it;
        if(!idx[u]) {
            tarjan(u);
            lowlink[v] = min(lowlink[v], lowlink[u]);
        }
        else if(onstack[u]) lowlink[v] = min(lowlink[v], lowlink[u]);
    }

    if(idx[v] == lowlink[v]) {
        nrctc++;
        do {
            x = s[ssize];
            sol[nrctc].insert(x);
            ctc[x] = nrctc;
            onstack[x] = 0; ssize--;
        } while(x != v);
    }
}

void drum_ctc(int i, int j) {
    for(set<int>::iterator it = sol[i].begin(); it != sol[i].end(); ++it) {
        int x = *it;
        for(set<int>::iterator it2 = sol[j].begin(); it2 != sol[j].end(); ++it2) {
            int y = *it2;
            if(exista_muchie(x, y, l)) {
                ls[i].insert(j);
                return;
            }
            if(exista_muchie(y, x, l)) {
                ls[j].insert(i);
                return;
            }
        }
    }
}

void simplifica_graf() {
    for(int i = 1; i < nrctc; i++)
        for(int j = i + 1; j <= nrctc; j++) drum_ctc(i, j);
}

void df(int x) {
    viz[x] = 1;

    for(int i = 1; i <= nrctc; i++) {
        if(exista_muchie(x, i, ls) && (!viz[i] || cost[i] > cost[x])) {
            cost[i] = cost[x];
            df(i);
        }
        if(exista_muchie(i, x, ls) && (!viz[i] || cost[i] > cost[x] + 1)) {
            cost[i] = cost[x] + 1;
            df(i);
        }
    }
}

int main() {

    citire_graf();

    for(int i = 1; i <= n; i++)
        if(!idx[i]) tarjan(i);

    simplifica_graf();

    fscanf(f, "%d", &t);
    for(int i = 1; i <= t; i++) {
        fscanf(f, "%d%d", &x, &y);

        memset(cost, 0, sizeof(cost));
        memset(viz, 0, sizeof(viz));
        df(ctc[x]);
        fprintf(g, "%d\n", cost[ctc[y]]);
    }

    fclose(f); fclose(g);

    return 0;
}