Pagini recente » Cod sursa (job #1418677) | Cod sursa (job #868959) | Cod sursa (job #2312881) | Cod sursa (job #1209206) | Cod sursa (job #906236)
Cod sursa(job #906236)
#include<cstdio>
#include<vector>
#include<queue>
#define vec G[nod][i]
#define NMAX 355
#define INF (1<<30)
using namespace std;
FILE *fin,*fout;
int CAPACITY[NMAX][NMAX],COST[NMAX][NMAX],n,source,destination;
int DIST[NMAX],k,FLUX[NMAX][NMAX],sw;
long long VOLMIN,sol;
short TT[NMAX];
vector<short> G[NMAX];
bool use[NMAX];
queue<short> Q;
void read()
{
fin=fopen("fmcm.in","r");
int m,z,c;
fscanf(fin,"%d %d %d %d",&n,&m,&source,&destination);
short x,y;
while(m--)
{
fscanf(fin,"%hi %hi %d %d",&x,&y,&z,&c);
G[x].push_back(y);
G[y].push_back(x);
CAPACITY[x][y]=z;
COST[x][y]=c;
COST[y][x]=-c;
}
fclose(fin);
}
void print()
{
fout=fopen("fmcm.out","w");
fprintf(fout,"%ld",sol);
fclose(fout);
}
void INIT()
{
for(int i=1;i<=n;i++)
{
DIST[i]=INF;
TT[i]=-1;
use[i]=0;
}
}
void bellman_ford()
{
int nod;
sw=0;
INIT();
DIST[source]=0;
Q.push(source);
use[source]=1;
while(!Q.empty())
{
nod=Q.front();
Q.pop();
for(size_t i=0;i<G[nod].size();i++)
if(CAPACITY[nod][vec]>FLUX[nod][vec] && DIST[nod]+COST[nod][vec]<DIST[vec])
{
DIST[vec]=DIST[nod]+COST[nod][vec];
TT[vec]=nod;
if(!use[vec])
{
use[vec]=1;
Q.push(vec);
}
}
use[nod]=0;
}
}
void flux_update()
{
int i;
if(DIST[destination]!=INF)
{
VOLMIN=INF;
sw=1;
for(i=destination;i!=source;i=TT[i])
VOLMIN=min(VOLMIN,1LL*(CAPACITY[TT[i]][i]-FLUX[TT[i]][i]));
for(i=destination;i!=source;i=TT[i])
{
FLUX[TT[i]][i]+=VOLMIN;
FLUX[i][TT[i]]-=VOLMIN;
}
sol+=VOLMIN*DIST[destination];
}
}
int main()
{
read();
sw=1;
while(sw)
{
bellman_ford();
flux_update();
}
print();
return 0;
}