Pagini recente » Cod sursa (job #726615) | Cod sursa (job #1812090) | Cod sursa (job #1107388) | Cod sursa (job #1203461) | Cod sursa (job #1166453)
#include <fstream>
#include <vector>
#include <deque>
#define inf 200000000
using namespace std;
ifstream f("fmcm.in");
ofstream g("fmcm.out");
int N, M, x, y, z, c, cap[352][352], flux[352][352], cost[352][352], D[352], prec[352], sink, source, sol;
bool inQ[352];
vector < int > G[352];
deque < int > Q;
int bellman(int sursa)
{
for (int i=1; i<=N; ++i)
D[i]=prec[i]=inf;
Q.push_back(sursa); inQ[sursa]=1; D[sursa]=0; prec[sursa]=sursa;
while (Q.size())
{
int nod=Q.front(); Q.pop_front(); inQ[nod]=0;
vector <int>::iterator it=G[nod].begin();
for (; it!=G[nod].end(); ++it)
if (D[nod]+cost[nod][*it]<D[*it] && cap[nod][*it]>flux[nod][*it])
{
D[*it]=D[nod]+cost[nod][*it]; prec[*it]=nod;
if (!inQ[*it]) inQ[*it]=1, Q.push_back(*it);
}
}
if (D[sink]!=inf)
{
int minflow=inf;
for (int i=sink; i!=sursa; i=prec[i])
if (cap[prec[i]][i]-flux[prec[i]][i]<minflow)
minflow=cap[prec[i]][i]-flux[prec[i]][i];
for (int i=sink; i!=sursa; i=prec[i])
flux[prec[i]][i]+=minflow, flux[i][prec[i]]-=minflow;
return minflow;
}
return 0;
}
int main()
{
f>>N>>M>>source>>sink;
for (int i=1; i<=M; ++i)
{
f>>x>>y>>z>>c; cap[x][y]=z;
cost[x][y]=c; cost[y][x]=-c;
G[x].push_back(y); G[y].push_back(x);
}
for (int flow=bellman(source); flow; flow=bellman(source))
sol+=flow*D[sink];
g<<sol<<'\n';
return 0;
}