Pagini recente » Cod sursa (job #1171559) | Cod sursa (job #959805) | Cod sursa (job #234170) | Cod sursa (job #1413063) | Cod sursa (job #403010)
Cod sursa(job #403010)
#include <algorithm>
#include <vector>
#include <bitset>
#include <queue>
using namespace std;
#define pb push_back
#define mp make_pair
#define INF 0x3f3f
#define sc second
#define fs first
#define DIM 355
vector <pair <short,short> > g[DIM];
short c[DIM][DIM],f[DIM][DIM];
short dst[DIM],t[DIM];
bitset <DIM> viz;
short n,m,s,d;
int cost;
struct cmp
{
bool operator () (const short &a,const short &b) const
{
return dst[a]>dst[b];
}
}; priority_queue <short,vector <short>,cmp> q;
void read ()
{
short i,x,y,z,w;
scanf ("%hd%hd%hd%hd",&n,&m,&s,&d);
for (i=1; i<=m; ++i)
{
scanf ("%hd%hd%hd%hd",&x,&y,&z,&w);
c[x][y]=z;
g[x].pb (mp (y,w));
g[y].pb (mp (x,-w));
}
}
bool dijkstra ()
{
vector <pair <short,short> > :: iterator it;
short nod;
memset (dst,INF,sizeof (dst));
dst[s]=0;
for (q.push (s); !q.empty (); )
{
viz[nod=q.top ()]=0; q.pop ();
for (it=g[nod].begin (); it!=g[nod].end (); ++it)
if (dst[it->fs]>dst[nod]+it->sc && c[nod][it->fs]>f[nod][it->fs])
{
dst[it->fs]=dst[nod]+it->sc;
t[it->fs]=nod;
if (!viz[it->fs])
{
q.push (it->fs);
viz[it->fs]=1;
}
}
}
if (dst[d]==INF)
return 0;
return 1;
}
void solve ()
{
int nrmin;
short i;
for (nrmin=INF; dijkstra (); nrmin=INF)
{
for (i=d; i!=s; i=t[i])
nrmin=min (nrmin,c[t[i]][i]-f[t[i]][i]);
for (i=d; i!=s; i=t[i])
{
f[t[i]][i]+=nrmin;
f[i][t[i]]-=nrmin;
}
cost+=nrmin*dst[d];
}
printf ("%d",cost);
}
int main ()
{
freopen ("fmcm.in","r",stdin);
freopen ("fmcm.out","w",stdout);
read ();
solve ();
return 0;
}