Pagini recente » Cod sursa (job #200019) | Istoria paginii utilizator/lacatus | Cod sursa (job #1854787) | Cod sursa (job #2808090) | Cod sursa (job #1409664)
#include <cstdio>
#include <vector>
#include <cstring>
#include <queue>
#define NMAX 355
#define INF 0x3f3f3f3f3f
using namespace std;
bool ap[NMAX];
int boss[NMAX],best[NMAX];
int n,m,x,y,z,c,surs,dest,costmax;
int cap[NMAX][NMAX],cost[NMAX][NMAX],flow[NMAX][NMAX];
vector <int>v[NMAX];
priority_queue <pair <int ,int > >q;
inline bool bellman(int surs)
{
memset(ap,0,sizeof(ap));
for (int i=1;i<=n;i++)best[i]=INF,boss[i]=0;
ap[surs]=1;
q.push(make_pair(0,surs));
best[surs]=0;
for (;!q.empty();q.pop())
{
int nod=q.top().second;
ap[nod]=0;
for (auto it:v[nod])
if (best[it]>best[nod]+cost[nod][it] && cap[nod][it]>0)
{
boss[it]=nod;
best[it]=best[nod]+cost[nod][it];
if (!ap[it])q.push(make_pair(-best[it],it));
}
}
return boss[dest]!=0;
}
int main()
{
freopen("fmcm.in","r",stdin);
freopen("fmcm.out","w",stdout);
scanf("%d %d %d %d",&n,&m,&surs,&dest);
for (int i=1;i<=m;i++)
{
scanf("%d %d %d %d",&x,&y,&z,&c);
v[x].push_back(y);
v[y].push_back(x);
cost[x][y]=c;
cost[y][x]=-c;
cap[x][y]=z;
}
for (;bellman(surs);)
{
int maxflow=INF;
for (int nod=dest;boss[nod];nod=boss[nod])
maxflow=(maxflow>cap[boss[nod]][nod]?cap[boss[nod]][nod]:maxflow);
for (int nod=dest;boss[nod];nod=boss[nod])
{
cap[boss[nod]][nod]-=maxflow;
cap[nod][boss[nod]]+=maxflow;
flow[boss[nod]][nod]+=maxflow;
costmax+=cost[boss[nod]][nod]*maxflow;
}
//printf("%d\n",costmax);
}
printf("%d",costmax);
fclose(stdin);
fclose(stdout);
}