Cod sursa(job #1508999)

Utilizator TibixbAndrei Tiberiu Tibixb Data 23 octombrie 2015 13:28:55
Problema Taramul Nicaieri Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.26 kb
#include<fstream>
#include<vector>
#include<bitset>
#include<deque>
using namespace std;
vector <int> L[203];
bitset<203> v;
deque <int> q;
int T[203], n, i, dout, din, j, C[203][203], nod, F[203][203], minim, ok, vecin, nrsol;
int bfs()
{
    ok=0;
    q.clear();
    v.reset();
    q.push_back(201);
    v[201]=1;
    while(!q.empty())
    {
        nod=*q.begin();
        if(nod==202)
        {
            ok=1;
            break;
        }
        for(int i=0; i<L[nod].size(); i++)
        {
            vecin=L[nod][i];
            if(v[vecin]==0 && C[nod][vecin]-F[nod][vecin]>0)
            {
                q.push_back(vecin);
                T[vecin]=nod;
                v[vecin]=1;
            }
        }
        q.pop_front();
    }
    return ok;
}
ifstream in("harta.in");
ofstream out("harta.out");
int main()
{
    in>>n;
    for(i=1; i<=n; i++)
    {
        in>>dout>>din;
        nrsol+=dout;
        L[201].push_back(i);
        L[i].push_back(201);
        C[201][i]=dout;
        L[i+n].push_back(202);
        L[202].push_back(i+n);
        C[i+n][202]=din;
    }
    for(i=1; i<=n; i++)
        for(j=1; j<=n; j++)
            if(i!=j){
                L[i].push_back(j+n);
                L[j+n].push_back(i);
                C[i][j+n]=1;
            }
    while(bfs())
    {
        for(i=0; i<L[202].size(); i++)
        {
            nod=L[202][i];
            if(v[nod]==1 && C[nod][202]-F[nod][202]>0)
            {
                minim=C[nod][202]-F[nod][202];
                for(; T[nod]!=0 && minim!=0; nod=T[nod])
                {
                    minim=min(minim, C[T[nod]][nod]-F[T[nod]][nod]);
                }
                if(minim!=0)
                {
                    nod=L[202][i];
                    F[nod][202]+=minim;
                    F[202][nod]-=minim;
                    for(; T[nod]!=0; nod=T[nod])
                    {
                        F[T[nod]][nod]+=minim;
                        F[nod][T[nod]]-=minim;
                    }
                }
            }
        }
    }
    out<<nrsol<<"\n";
    for(i=1; i<=n; i++)
        for(j=n+1; j<=2*n; j++)
            if(F[i][j]==1)
                out<<i<<" "<<j-n<<"\n";
    return 0;
}