bzoj1834(网络流加费用流)

1834: [ZJOI2010]network 网络扩容

Time Limit: 3 Sec Memory Limit: 64 MB

Submit: 2390 Solved: 1206

[Submit][Status][Discuss]

Description

给定一张有向图，每条边都有一个容量C和一个扩容费用W。这里扩容费用是指将容量扩大1所需的费用。求： 1、 在不扩容的情况下，1到N的最大流； 2、 将1到N的最大流增加K所需的最小扩容费用。

Input

输入文件的第一行包含三个整数N,M,K，表示有向图的点数、边数以及所需要增加的流量。 接下来的M行每行包含四个整数u,v,C,W，表示一条从u到v，容量为C，扩容费用为W的边。

Output

输出文件一行包含两个整数，分别表示问题1和问题2的答案。

Sample Input

5 8 2

1 2 5 8

2 5 9 9

5 1 6 2

5 1 1 8

1 2 8 7

2 5 4 9

1 2 1 1

1 4 2 1

Sample Output

13 19

30%的数据中，N<=100

100%的数据中，N<=1000，M<=5000，K<=10

解题思路：第一问就是最大流，然后第二问，可以对在建m条边，只是产生费用，再来一个源点，到一的容量为K。跑一边费用流。

#include<cstdio>

#include<cstring>

#include<algorithm>

#include<iostream>

using namespace std;

int n,m,k,len,ans;

int from[21000],to[21000],cc[21000],ww[21000],next[21000];

int h[1100];

int dis[1100];

int q[1000000];

bool og[1100];

int pre[1100];

int x[5100],y[5100],c[5100],w[5100];

inline int read()

{

char y; int x=0,f=1; y=getchar();

while (y<'0' || y>'9') {if (y=='-') f=-1; y=getchar();}

while (y>='0' && y<='9') {x=x*10+int(y)-48; y=getchar();}

return x*f;

}

void insert(int x,int y,int c,int w)

{

++len;

from[len]=x; to[len]=y; cc[len]=c; ww[len]=w; next[len]=h[x]; h[x]=len;

}

int dicnic(int sum,int now)

{

if (now==n) return sum;

int u=h[now]; int sug=0;

while (u!=0)

{

if (cc[u]>0 && dis[now]+1==dis[to[u]])

{

int s=dicnic(min(cc[u],sum-sug),to[u]);

if (s>0)

{

sug+=s;

cc[u]-=s; cc[u^1]+=s;

if (sug==sum) return sum;

}

}

u=next[u];

}

if (sug==0) dis[now]=-1;

return sug;

}

bool bfs()

{

memset(dis,-1,sizeof(dis));

dis[1]=1; int tail=0,head=0; ++tail; q[tail]=1;

while (head<tail)

{

++head;

int u=h[q[head]];

while (u!=0)

{

if (cc[u]>0 && dis[to[u]]==-1)

{

dis[to[u]]=dis[q[head]]+1;

++tail; q[tail]=to[u];

}

u=next[u];

}

}

if (dis
!=-1) return true; else return false;

}

bool spfa()

{

memset(dis,0x7e,sizeof(dis));

memset(og,true,sizeof(og));

og[0]=false;

dis[0]=0; int tail=1,head=0; q[tail]=0;

while (head<tail)

{

++head;

int u=h[q[head]];

while (u!=0)

{

if (cc[u]>0 && dis[to[u]]>dis[q[head]]+ww[u])

{

pre[to[u]]=u;

dis[to[u]]=dis[q[head]]+ww[u];

if (og[to[u]])

{

og[to[u]]=false;

++tail; q[tail]=to[u];

}

}

u=next[u];

}

og[q[head]]=true;

}

if (dis
<1000000) return true;else return false;

}

void mcf()

{

int now=n; int zan=0x7fffffff;

while (now!=0)

{

zan=min(zan,cc[pre[now]]);

now=from[pre[now]];

}

now=n;

while (now!=0)

{

cc[pre[now]]-=zan; cc[pre[now]^1]+=zan;

ans+=zan*ww[pre[now]];

now=from[pre[now]];

}

}

int main()

{

n=read(); m=read(); k=read();

len=1;

for (int i=1;i<=m;++i)

{

x[i]=read(); y[i]=read(); c[i]=read(); w[i]=read();

insert(x[i],y[i],c[i],0); insert(y[i],x[i],0,0);

}

int sum=0;

while (bfs())

{

sum+=dicnic(0x7fffffff,1);

}

printf("%d ",sum);

insert(0,1,k,0); insert(1,0,-k,0);

for (int i=1;i<=m;++i)

{

insert(x[i],y[i],k,w[i]); insert(y[i],x[i],0,-w[i]);

}

ans=0;

while (spfa())

{

mcf();

}

printf("%d",ans);

}