Pagini recente » Cod sursa (job #554018) | Cod sursa (job #2436768) | Cod sursa (job #3230743) | Cod sursa (job #2948892) | Cod sursa (job #2887130)
#include <bits/stdc++.h>
using namespace std;
ifstream fin("fmcm.in");
ofstream fout("fmcm.out");
const int NMAX = 355;
const int INF = (1<<29);
int capacitate[NMAX][NMAX],flux[NMAX][NMAX];
int cost[NMAX][NMAX],dist[NMAX],tata[NMAX];
int n,m,S,D,rasp,x,y,cap,z;
vector <int> v[NMAX];
bool ver[NMAX];
void citire(){
fin >> n >> m >> S >> D;
for(int i=1;i<=m;i++){
fin >> x >> y >> cap >> z;
v[x].push_back(y);
v[y].push_back(x);
capacitate[x][y]=cap;
cost[x][y]=z;
cost[y][x]=-z;
}
}
void init(){
for(int i=1;i<=n;i++) {
ver[i]=false;
dist[i]=INF;
}
}
bool muchie_disponibila(int node1,int node2){
if(capacitate[node1][node2]>flux[node1][node2])
return true;
return false;
}
bool Bellman(){
init();
queue <int> q;
q.push(S);
dist[S]=0;
while(!q.empty()){
int node = q.front();
q.pop();
ver[node]=false;
for(int i=0;i<v[node].size();i++){
int vecin=v[node][i];
if(muchie_disponibila(node,vecin)==true){
if(dist[vecin]>dist[node]+cost[node][vecin]){
dist[vecin]=dist[node]+cost[node][vecin];
if(ver[vecin]==false){
q.push(vecin);
ver[vecin]=true;
}
tata[vecin]=node;
}
}
}
}
if(dist[D]!=INF) return true;
return false;
}
void solve(){
while(Bellman()==true){
int aux = D;
int minim = INF;
while(aux!=S){
if(capacitate[tata[aux]][aux]-flux[tata[aux]][aux]<minim)
minim = capacitate[tata[aux]][aux]-flux[tata[aux]][aux];
aux=tata[aux];
}
rasp+=dist[D]*minim;
aux = D;
while(aux!=S){
flux[tata[aux]][aux]+=minim;
flux[aux][tata[aux]]-=minim;
aux=tata[aux];
}
}
}
void afis(){
fout << rasp;
}
int main()
{
citire();
solve();
afis();
return 0;
}