Pagini recente » Cod sursa (job #2069962) | Cod sursa (job #606197) | Cod sursa (job #1410336) | Cod sursa (job #692111) | Cod sursa (job #2959202)
#include <iostream>
#include <bits/stdc++.h>
#include <fstream>
using namespace std;
using namespace std::chrono;
//ifstream in("maxflow.in");
//ofstream out("maxflow.out");
ifstream in("harta.in");
ofstream out("harta.out");
vector<vector<int>> adjList;
vector<pair<int,int>> degrees;
unordered_map<int, int> cap[205];
unordered_map<int,int> sol;
vector<int> parents;
struct muchie{
};
int n,m,e;
void read(){
in >> n;
adjList.resize(2 * n + 2);
// cap.resize(n + m + 3, vector<int>(n + m + 3, 0));
int x, y;
degrees.emplace_back(make_pair(0,n));
while(in>>x>>y){
degrees.emplace_back(make_pair(x,y));
}
for(int i = 1; i <= n ;i++)
for(int j = n + 1; j<= n * 2; j++)
{
if(i + n != j){
adjList[i].emplace_back(j);
adjList[j].emplace_back(i);
cap[i][j]++;
}
}
for(int i = 1; i <= n; i++)
{
adjList[0].emplace_back(i);
// adjList[i].emplace_back(0);
cap[0][i] = degrees[i].first;
}
for(int i = 1; i + n <= n * 2; i++)
{
adjList[i + n].emplace_back(n * 2 + 1);
adjList[n * 2 + 1].emplace_back(i + n);
cap[i + n][n * 2 + 1] = degrees[i].second;
}
}
int BFS(int source, int dest, vector<int>& parent)
{
fill(parent.begin(), parent.end(), -1);
parent[source] = 0;
queue<int> q;
q.push(source);
while(!q.empty()) {
int pNode = q.front(); //nodul parinte
q.pop();
for (auto el: adjList[pNode])
{
int cNode = el; // copiii nodului parinte actual
if(cNode != pNode && parent[cNode] == -1 && cap[pNode][cNode] > 0)
{
parent[cNode] = pNode;
if(cNode == dest)
{
return 1;
}
//adaugam urmatorul nod ce urmeaza sa fie vizitat
q.push({cNode});
}
}
}
return 0;
}
int FordFulkerson(int source, int dest)
{
vector<int> parent(n * 2 + 2, -1);
parent[source] = -2;
int mxFlow = 0;
while(BFS(source,dest,parent)) {
for (auto el: adjList[dest]) {
if (parent[el] != -1) {
int mnCap = INT_MAX;
int cNode = dest;
// parent[dest] = el;
while (cNode != source) {
mnCap = min(mnCap, cap[parent[cNode]][cNode]);
cNode = parent[cNode];
}
cNode = dest;
while (cNode != source){
cap[cNode][parent[cNode]] += mnCap;
cap[parent[cNode]][cNode] -= mnCap;
cNode = parent[cNode];
}
mxFlow += mnCap;
}
}
}
return mxFlow;
}
int main()
{
auto start = high_resolution_clock::now();
read();
out<<FordFulkerson(0, 2*n+1)<<endl;
for(int i = 1; i<=n;i++)
for(int j = n + 1; j<= n * 2; j++){
if(i + n != j && cap[i][j] == 0)
out<<i<<" "<<j - n<<endl;
}
auto stop = high_resolution_clock::now();
auto duration = duration_cast<microseconds>(stop - start);
cout<<"Durata: "<< duration.count();
return 0;
}