Pagini recente » Cod sursa (job #441803) | Cod sursa (job #2547807) | Cod sursa (job #2842026) | Cod sursa (job #2120557) | Cod sursa (job #1865567)
#include <cstdio>
#include <queue>
#define ll long long
#define maxn 505
#define maxx 1000010
using namespace std;
vector <int> ls[maxn];
ll cap[maxn][maxn], f[maxn][maxn], cs[maxn][maxn], dist[maxn], fmin, sol;
const ll f_mare = 2e11;
int n, m, inc, sf, x, y, z, t, tt[maxn], g[maxn], i;
int coada[maxx], st, dr;
bool viz[maxn], ok;
ll bf() {
for (i = 1; i <= n; i++) {
dist[i] = f_mare;
tt[i] = -1;
viz[i] = 0;
}
dist[inc] = 0;
st=dr=0, coada[dr] = inc;
viz[inc] = 1;
while (st <= dr) {
x = coada[st++];
for (i = 0; i < g[x]; i++) {
y = ls[x][i];
if (dist[y] > dist[x] + cs[x][y] && cap[x][y]-f[x][y] > 0) {
if (!viz[y]) {coada[++dr] = y; viz[y] = 1;}
dist[y] = dist[x] + cs[x][y];
tt[y] = x;
}
}
viz[x]=0;
}
if (dist[y] < f_mare) {
ok = 1;
fmin = f_mare;
y = sf;
while (y != inc) {
fmin = min(fmin, cap[tt[y]][y] - f[tt[y]][y]);
y = tt[y];
}
y = sf;
while (y != inc) {
f[tt[y]][y] += fmin;
f[y][tt[y]] -= fmin;
y = tt[y];
}
return dist[sf]*fmin;
}
return 0;
}
void calc() {
ok=1;
while (ok) {
ok = 0;
sol += bf();
}
}
int main() {
FILE *f = fopen("fmcm.in", "r");
FILE *G = fopen("fmcm.out", "w");
fscanf(f, "%d %d %d %d\n", &n, &m, &inc, &sf);
while (m--) {
fscanf(f, "%d %d %d %d\n", &x, &y, &z, &t);
ls[x].push_back(y);
ls[y].push_back(x);
cap[x][y] = z;
cs[x][y] = t;
cs[y][x] = -t;
}
for (i = 1; i <= n; i++)
g[i] = ls[i].size();
calc();
fprintf(G,"%d",sol);
return 0;
}