Pagini recente » Cod sursa (job #2499962) | Cod sursa (job #1219242) | Cod sursa (job #2737975) | Cod sursa (job #1823138) | Cod sursa (job #2047318)
#include <iostream>
#include <fstream>
#include <vector>
#include <queue>
using namespace std;
const int MAX_N = 350;
int cost[1 + MAX_N][1 + MAX_N];
int flow[1 + MAX_N][1 + MAX_N];
vector<int> edges[1 + MAX_N];
int dist[1 + MAX_N];
int from[1 + MAX_N];
const int INF = 1e9;
int main() {
ifstream fin( "fmcm.in" );
ofstream fout( "fmcm.out" );
int n, m, S, D;
fin >> n >> m >> S >> D;
for ( int i = 0; i < m; i ++ ) {
int u, v, c, z;
fin >> u >> v >> c >> z;
edges[u].push_back( v );
edges[v].push_back( u );
cost[u][v] = z;
cost[v][u] = -z;
flow[u][v] = c;
}
for ( int i = 1; i <= n; i ++ )
if ( i != S )
dist[i] = INF;
queue<int> bellman;
bellman.push( S );
while ( bellman.size() ) {
int u = bellman.front();
bellman.pop();
for ( int v : edges[u] )
if ( flow[u][v] && dist[u] + cost[u][v] < dist[v] ) {
dist[v] = dist[u] + cost[u][v];
bellman.push( v );
}
}
int min_cost = 0, max_flow = 0;
bool improve = true;
do {
for ( int i = 1; i <= n; i ++ )
from[i] = 0;
from[S] = -1;
priority_queue<pair<int, int>> dijkstra;
dijkstra.push( make_pair( 0, S ) );
while ( dijkstra.size() && dijkstra.top().second != D ) {
int u = dijkstra.top().second, z = dijkstra.top().first;
dijkstra.pop();
for ( int v : edges[u] )
if ( flow[u][v] > 0 && !from[v] ) {
from[v] = u;
dijkstra.push( make_pair( z + dist[u] - dist[v] + cost[u][v], v ) );
}
}
if ( dijkstra.size() ) {
int f = 1e9;
for ( int i = D; i != S; i = from[i] )
f = min( f, flow[from[i]][i] );
for ( int i = D; i != S; i = from[i] ) {
flow[from[i]][i] -= f;
flow[i][from[i]] += f;
min_cost += cost[from[i]][i] * f;
}
max_flow += f;
} else
improve = false;
} while ( improve );
fout << min_cost << endl;
return 0;
}