Pagini recente » Cod sursa (job #3267273) | Cod sursa (job #1505586) | Cod sursa (job #1577962) | Cod sursa (job #782478) | Cod sursa (job #3125013)
#include <bits/stdc++.h>
using namespace std;
ifstream fin("fmcm.in");
ofstream fout("fmcm.out");
typedef pair<int,int> pii;
int n,m,S,D,Fans,Fcost;
int cap[355][355],cost[355][355];
vector<int> muchii[355];
int dist[355],from[355],real_d[355],old_d[355],use[355];
bool inq[355];
priority_queue<pii,vector<pii>,greater<pii>> pq;
void bellman()
{
for(int i=1;i<=n;i++)
old_d[i]=1e9;
old_d[S]=0;
queue<int> coada;
coada.push(S);
while(!coada.empty())
{
int nod=coada.front();
coada.pop();
for(int i:muchii[nod])
if(cap[nod][i]>0&&old_d[i]>old_d[nod]+cost[nod][i])
{
old_d[i]=old_d[nod]+cost[nod][i];
coada.push(i);
}
}
}
bool flux()
{
for(int i=1;i<=n;i++)
{
dist[i]=1e9;
real_d[i]=1e9;
from[i]=0;
inq[i]=0;
use[i]=0;
}
dist[S]=0;
real_d[S]=0;
pq.push({dist[S],S});
while(!pq.empty())
{
int nod=pq.top().second;
pq.pop();
if(use[nod])
continue;
use[nod]=1;
for(int i:muchii[nod])
if(!use[i]&&cap[nod][i]>0&&dist[i]>dist[nod]+cost[nod][i]+old_d[nod]-old_d[i])
{
dist[i]=dist[nod]+cost[nod][i]+old_d[nod]-old_d[i];
//assert(dist[i]>=0);
real_d[i]=real_d[nod]+cost[nod][i];
from[i]=nod;
pq.push({dist[i],i});
}
}
if(real_d[D]==1e9)
return 0;
for(int i=1;i<=n;i++)
old_d[i]=real_d[i];
int minim=1e9;
for(int x=D;x!=S;x=from[x])
minim=min(minim,cap[from[x]][x]);
Fans+=minim;
Fcost+=real_d[D]*minim;
for(int x=D;x!=S;x=from[x])
{
cap[from[x]][x]-=minim;
cap[x][from[x]]+=minim;
}
return 1;
}
int main()
{
ios_base::sync_with_stdio(false);
fin.tie(0);
fin>>n>>m>>S>>D;
for(int i=1;i<=m;i++)
{
int a,b,c,z;
fin>>a>>b>>c>>z;
cap[a][b]=c;
cost[a][b]=z;
cost[b][a]=-z;
muchii[a].push_back(b);
muchii[b].push_back(a);
}
bellman();
while(flux());
fout<<Fcost;
return 0;
}