bzoj5108 数据_BZOJ-1834 网络扩容 最小费用最大流+最大流+乱搞

1834: [ZJOI2010]network 网络扩容

Time Limit: 3 Sec Memory Limit: 64 MB

Submit: 2269 Solved: 1136

[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

HINT

Source

Day1

这道题啊，似乎不是很复杂，起码省去了bt建图，充其量是个模板堆上，随便乱搞几行，直接正解。(ZJOI中最水的了吧？？)

先按照读入连边。(第一问最大流时不需要费用，可以先存下来，为第二问准备)

Dinic模板一套，第一问A

在第一问的参与网络上，建边。

每两个点相连，边权为inf，费用为之前存下来的。

最后再建一个源，连向1，容量为k，费用为0

zkw模板一套，第二问A

。。。。。

code:(写的冗余QwQ)

#include

#include

#include

#include

using namespace std;

#define inf 0x7fffffff

struct data{

int from,to,next,c,v,co;

}edge[100010];

int q[100010],h,t;

int dis[100010];

int head[100010]={0},cnt=1;

bool visit[100010],mark[100010];

int n,m,k;

int S,T;

int ans;

int read()

{

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

while(ch'9'){if(ch=='-')f=-1;ch=getchar();}

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

return x*f;

}

void add(int u,int v,int cap,int cost)

{

cnt++;edge[cnt].from=u;edge[cnt].to=v;

edge[cnt].v=cap;edge[cnt].co=cost;

edge[cnt].next=head[u];head[u]=cnt;

}

void insert(int u,int v,int cap,int cost)

{

add(u,v,cap,cost);add(v,u,0,-cost);

}

void add_edge(int u,int v,int cap,int cost)

{

cnt++;edge[cnt].from=u;edge[cnt].to=v;

edge[cnt].v=cap;edge[cnt].c=cost;

edge[cnt].next=head[u];head[u]=cnt;

}

void insert_edge(int u,int v,int cap,int cost)

{

add_edge(u,v,cap,cost);add_edge(v,u,0,-cost);

}

void init()

{

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

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

{

int u,v,c,w;

u=read();v=read();c=read();w=read();

insert(u,v,c,w);

}

S=1;T=n;

}

bool bfs()

{

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

q[1]=S;dis[S]=1;h=0;t=1;

while (h

{

int j=q[++h],i=head[j];

while (i)

{

if (dis[edge[i].to]<0 && edge[i].v>0)

{

dis[edge[i].to]=dis[j]+1;

q[++t]=edge[i].to;

}

i=edge[i].next;

}

}

if (dis[T]>0) return true;

else return false;

}

int dfs(int loc,int low)

{

if (loc==T) return low;

int w,used=0;

for (int i=head[loc]; i; i=edge[i].next)

{

if (edge[i].v>0 && dis[edge[i].to]==dis[loc]+1)

{

w=dfs(edge[i].to,min(low-used,edge[i].v));

edge[i].v-=w;edge[i^1].v+=w;

used+=w;if (used==low) return low;

}

}

if (!used) dis[loc]=-1;

return used;

}

int dinic()

{

int tmp=0;

while (bfs())

{

tmp+=dfs(S,inf);

}

return tmp;

}

void problem_1()

{

int tmp=dinic();

printf("%d ",tmp);

}

void make()

{

int num=cnt;

for (int i=2; i<=num; i+=2)

insert_edge(edge[i].from,edge[i].to,inf,edge[i].co);

insert(0,1,k,0);

S=0;

}

bool spfa()

{

memset(visit,0,sizeof(visit));

for (int i=S; i<=T; i++) dis[i]=inf;

h=0,t=1;

q[0]=T;visit[T]=1;dis[T]=0;

while (h

{

int now=q[h];h++;visit[now]=0;

for (int i=head[now]; i; i=edge[i].next)

if (edge[i^1].v && dis[now]-edge[i].c

{

dis[edge[i].to]=dis[now]-edge[i].c;

if (!visit[edge[i].to])

{

visit[edge[i].to]=1;

q[t++]=edge[i].to;

}

}

}

return dis[S]!=inf;

}

int dfs1(int loc,int low)

{

mark[loc]=1;

if (loc==T) return low;

int w,used=0;

for (int i=head[loc]; i; i=edge[i].next)

if (dis[edge[i].to]==dis[loc]-edge[i].c && edge[i].v && !mark[edge[i].to])

{

w=dfs1(edge[i].to,min(low-used,edge[i].v));

ans+=w*edge[i].c;

edge[i].v-=w;edge[i^1].v+=w;

used+=w;if (used==low) return low;

}

return used;

}

void zkw()

{

int tmp=0;

while (spfa())

{

mark[T]=1;

while (mark[T])

{

memset(mark,0,sizeof(mark));

tmp+=dfs1(S,inf);

}

}

}

void problem_2()

{

make();

zkw();

printf("%d",ans);

}

int main()

{

init();

problem_1();

problem_2();

return 0;

}