bzoj 1834: [ZJOI2010]network 網路擴容 -- 最大流+費用流

標籤：span   最大流   擴大   網路   network   mat   stream   ++   efi   

1834: [ZJOI2010]network 網路擴容Time Limit: 3 Sec  Memory Limit: 64 MBDescription給定一張有向圖，每條邊都有一個容量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 Input5 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 Output13 19
30%的資料中，N<=100
100%的資料中，N<=1000，M<=5000，K<=10
HINT 先跑一遍最大流然後在殘留網路建邊對於每條邊重建立一條容量inf的邊（注意原邊不變且費用為0）然後就直接跑費用流 

#include<map>#include<cmath>#include<queue>#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>using namespace std;#define inf 100000007#define ll long long#define N 50010inline int rd(){    int x=0,f=1;char ch=getchar();    while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();}    while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();}    return x*f;}int T;int fro[N],to[N],lj[N],v[N],w[N],fa[N],cnt=1,tw[N];void add(int a,int b,int c,int d){fro[++cnt]=lj[a];to[cnt]=b;fa[cnt]=a;v[cnt]=c;tw[cnt]=d;lj[a]=cnt;}void ins(int a,int b,int c,int d){add(a,b,c,d);add(b,a,0,-d);}int dis[N],q[N],from[N],ans,h,t;bool vs[N];bool bfs(){    memset(dis,0,sizeof(dis));    dis[1]=q[0]=t=1;h=0;    int tp;    while(h!=t)    {        tp=q[h++]; if(h==N) h=0;        for(int i=lj[tp];i;i=fro[i])        {            if(!dis[to[i]]&&v[i])            {                dis[to[i]]=dis[tp]+1;                q[t++]=to[i]; if(t==N) t=0;            }        }    }    return dis[T]?1:0;}int dfs(int x,int p){    if(x==T) return p;    int tp,res=0;    for(int i=lj[x];i;i=fro[i])    {        if(v[i]&&dis[to[i]]==dis[x]+1)        {            tp=dfs(to[i],min(p-res,v[i]));            v[i]-=tp;             v[i^1]+=tp;            res+=tp;            if(res==p) return p;        }    }    if(res==0) dis[x]=0;    return res;}void dinic(){while(bfs()) ans+=dfs(1,inf);}bool spfa(){    memset(dis,0x3f,sizeof(dis));    int x;    h=q[0]=dis[0]=0;t=vs[0]=1;    while(h!=t)    {        x=q[h++]; if(h==N) h=0;        for(int i=lj[x];i;i=fro[i])        {            if(v[i]&&dis[to[i]]>dis[x]+w[i])            {                dis[to[i]]=dis[x]+w[i];                from[to[i]]=i;                if(!vs[to[i]])                {                    vs[to[i]]=1;                    q[t++]=to[i];if(t==N) t=0;                }            }        }        vs[x]=0;    }    return dis[T]<inf;}void qaz(){    int tmp=inf;    for(int i=from[T];i;i=from[fa[i]]) tmp=min(tmp,v[i]);    for(int i=from[T];i;i=from[fa[i]])    {        v[i]-=tmp;v[i^1]+=tmp;        ans+=w[i]*tmp;    }}void add2(int a,int b,int c,int d){fro[++cnt]=lj[a];to[cnt]=b;fa[cnt]=a;v[cnt]=c;w[cnt]=d;lj[a]=cnt;}void ins2(int a,int b,int c,int d){add2(a,b,c,d);add2(b,a,0,-d);}void build(){    int tp=cnt;    for(int i=2;i<=tp;i+=2) ins2(fa[i],to[i],inf,tw[i]);}int n,m,k,a,b,c,d;int main(){    T=n=rd();m=rd();k=rd();    for(int i=1;i<=m;i++)    {        a=rd();b=rd();c=rd();d=rd();        ins(a,b,c,d);    }    dinic();    printf("%d ",ans);    ans=0;build();    ins(0,1,k,0);    while(spfa()) qaz();    printf("%d\n",ans);    return 0;}

 

 

bzoj 1834: [ZJOI2010]network 網路擴容 -- 最大流+費用流