Pagini recente » Cod sursa (job #965147) | Cod sursa (job #1103353) | Cod sursa (job #422468) | Cod sursa (job #2074519) | Cod sursa (job #1636993)
#include<fstream>
#include<vector>
#include<cstring>
#include<queue>
#define inf 0x7f7f7f7f
using namespace std;
vector<int> L[360];
queue<int> q;
bool u[360];
int n, m, s, d, i, x, y, z, c, D[360], C[360][360], F[360][360], Z[360][360], nod, vecin, cost, minim, T[360];
int bellmanford()
{
memset(D, 127, sizeof(D));
D[s]=0;
memset(u, 0, sizeof(u));
q.push(s);
u[s]=1;
while(!q.empty())
{
nod=q.front();
for(i=0; i<L[nod].size(); i++)
{
vecin=L[nod][i];
if(D[vecin]>D[nod]+Z[nod][vecin] && C[nod][vecin]-F[nod][vecin]>0)
{
D[vecin]=D[nod]+Z[nod][vecin];
T[vecin]=nod;
if(u[vecin]==0)
{
q.push(vecin);
u[vecin]=1;
}
}
}
u[nod]=0;
q.pop();
}
return (D[d]!=inf);
}
void fmcm()
{
while(bellmanford())
{
minim=C[T[d]][d]-F[T[d]][d];
for(nod=T[d]; T[nod]!=0 && minim!=0; nod=T[nod])
minim=min(minim, C[T[nod]][nod]-F[T[nod]][nod]);
if(minim!=0)
{
cost+=minim*Z[T[d]][d];
F[T[d]][d]+=minim;
F[d][T[d]]-=minim;
for(nod=T[d]; T[nod]!=0; nod=T[nod])
{
cost+=minim*Z[T[nod]][nod];
F[T[nod]][nod]+=minim;
F[nod][T[nod]]-=minim;
}
}
}
}
ifstream in("fmcm.in");
ofstream out("fmcm.out");
int main()
{
in>>n>>m>>s>>d;
for(i=1; i<=m; i++)
{
in>>x>>y>>c>>z;
L[x].push_back(y);
C[x][y]=c;
Z[x][y]=z;
L[y].push_back(x);
Z[y][x]=-z;
}
fmcm();
out<<cost;
return 0;
}