Pagini recente » Cod sursa (job #1395845) | Cod sursa (job #2957712) | Cod sursa (job #2487138) | Cod sursa (job #1245336) | Cod sursa (job #1678566)
#include <iostream>
#include <fstream>
#include <vector>
#include <queue>
#define oo 2147483647
using namespace std;
int n, s, d;
vector<int> G[360];
int c[360][360];
int cost[360][306];
int ds[360];
int real[360];
int mn[360];
int tata[360];
bool vis[360];
priority_queue< pair< int, int> > pq;
int fluxCost;
int dijkstra() {
int aux;
for(int i = 1; i <= n; i ++) ds[i] = oo;
ds[s] = real[s] = 0;
pq.push({0, s});
while(!pq.empty()) {
int kost = -pq.top().first;
int act = pq.top().second;
pq.pop();
if(kost != ds[act]) continue;
for(auto it: G[act]) {
aux = ds[act] + cost[act][it] + mn[act] - mn[it];
if(c[act][it] && aux < ds[it]) {
ds[it] = aux;
real[it] = real[act] + cost[act][it];
tata[it] = act;
pq.push({-ds[it], it});
}
}
}
for(int i = 1; i <= n; i ++) mn[i] = real[i];
if(ds[d] == oo) return false;
int cmn = oo;
for(int i = d; tata[i]; i = tata[i])
cmn = min(cmn, c[tata[i]][i]);
for(int i = d; tata[i]; i = tata[i]) {
c[tata[i]][i] -= cmn;
c[i][tata[i]] += cmn;
}
fluxCost += cmn * real[d];
return true;
}
queue<int> q;
bool bellmanFord() {
int aux;
for(int i = 1; i <= n; i ++) mn[i] = oo;
mn[s] = 0;
vis[s] = true;
q.push(s);
while(!q.empty()) {
int act = q.front();
vis[act] = false;
q.pop();
for(auto it: G[act]) {
aux = mn[act] + cost[act][it];
if(c[act][it] && mn[it] > aux) {
mn[it] = aux;
if(!vis[it]) {
vis[it] = true;
q.push(it);
}
}
}
}
if(mn[d] == oo)
return false;
return true;
}
void flux() {
if(!bellmanFord()) return ;
while(dijkstra());
}
int main()
{
ifstream f("fmcm.in");
ofstream g("fmcm.out");
int m, x, y;
f >> n >> m >> s >> d;
for(; m; m --) {
f >> x >> y;
f >> c[x][y] >> cost[x][y];
cost[y][x] = -cost[x][y];
G[x].push_back(y);
G[y].push_back(x);
}
flux();
g << fluxCost << "\n";
return 0;
}