传送门

网络流水题啊。


第一问直接放心跑最大流(本来还以为有什么tricktricktrick)。

第二问就直接把原来的边(u,v,c,w)(u,v,c,w)(u,v,c,w)变成(u,v,c,0)(u,v,c,0)(u,v,c,0)和(u,v,inf,w)(u,v,inf,w)(u,v,inf,w),然后把ttt拆点限制流量跑费用流就行了。

代码:

#include<bits/stdc++.h>
#define N 1005
#define M 10005
using namespace std;
inline int read(){
	int ans=0;
	char ch=getchar();
	while(!isdigit(ch))ch=getchar();
	while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
	return ans;
}
int s,t,n,m,k,first[N],cur[N],d[N],flow[N],pos[N],pred[N],cnt;
bool in[N];
struct Edge{int u,v,w,c;}tt[M];
struct edge{int v,next,w,c;}e[M<<1];
inline void addedge(int u,int v,int c,int w){e[++cnt].v=v,e[cnt].next=first[u],e[cnt].c=c,e[cnt].w=w,e[cnt].v=v,first[u]=cnt;}
inline void add(int u,int v,int c,int w){addedge(u,v,c,w),addedge(v,u,0,-w);}
inline void init(){memset(first,-1,sizeof(first)),cnt=-1,s=1,t=n;}
inline bool bfs(){
	queue<int>q;
	memset(d,-1,sizeof(d)),d[s]=0,q.push(s);
	while(!q.empty()){
		int x=q.front();
		q.pop();
		for(int i=first[x];~i;i=e[i].next){
			int v=e[i].v;
			if(!e[i].c||~d[v])continue;
			d[v]=d[x]+1,q.push(v);
		}
	}
	return ~d[t];
}
inline int dfs(int x,int f){
	if(!f||x==t)return f;
	int flow=f;
	for(int&i=cur[x];~i;i=e[i].next){
		int v=e[i].v;
		if(!flow)return f;
		if(e[i].c&&d[v]==d[x]+1){
			int tmp=dfs(v,min(flow,e[i].c));
			if(!tmp)d[v]=-1;
			e[i].c-=tmp,e[i^1].c+=tmp,flow-=tmp;
		}
	}
	return f-flow;
}
inline bool spfa(){
	queue<int>q;
	memset(d,0x3f3f3f3f,sizeof(d)),memset(in,false,sizeof(in)),q.push(s),d[s]=0,in[s]=true,pred[t]=-1,flow[s]=0x3f3f3f3f;
	while(!q.empty()){
		int x=q.front();
		q.pop(),in[x]=false;
		for(int i=first[x];~i;i=e[i].next){
			int v=e[i].v;
			if(e[i].c&&d[v]>d[x]+e[i].w){
				d[v]=d[x]+e[i].w,pred[v]=x,pos[v]=i,flow[v]=min(flow[x],e[i].c);
				if(!in[v])in[v]=true,q.push(v);
			}
		}
	}
	return d[t]!=0x3f3f3f3f;
}
inline int dinic(){
	int ret=0;
	while(bfs()){
		for(int i=1;i<=n;++i)cur[i]=first[i];
		ret+=dfs(s,0x3f3f3f3f);
	}
	return ret;
}
inline int mcmf(){
	int ret=0;
	for(int w=t;spfa();w=t){
		ret+=flow[t]*d[t];
		while(w!=s)e[pos[w]].c-=flow[t],e[pos[w]^1].c+=flow[t],w=pred[w];
	}
	return ret;
}
int main(){
	n=read(),m=read(),k=read(),init();
	for(int i=1;i<=m;++i){
		tt[i].u=read(),tt[i].v=read(),tt[i].c=read(),tt[i].w=read();
		add(tt[i].u,tt[i].v,tt[i].c,0);
	}
	int tmp=dinic();
	printf("%d ",tmp);
	memset(first,-1,sizeof(first)),cnt=-1;
	for(int i=1;i<=m;++i)add(tt[i].u,tt[i].v,tt[i].c,0),add(tt[i].u,tt[i].v,0x3f3f3f3f,tt[i].w);
	add(n,n+1,k+tmp,0),t=n+1,printf("%d",mcmf());
	return 0;
}

2018.10.13 bzoj1834: [ZJOI2010]network 网络扩容(最大流+费用流)的更多相关文章

  1. BZOJ1834 [ZJOI2010]network 网络扩容(最小费用最大流)

    挺直白的构图..最小费用最大流的定义. #include<cstdio> #include<cstring> #include<queue> #include< ...

  2. BZOJ 1834: [ZJOI2010]network 网络扩容(网络流+费用流)

    一看就知道是模板题= = ,不说什么了= = PS:回去搞期末了,暑假再来刷题了 CODE: #include<cstdio> #include<iostream> #incl ...

  3. BZOJ1834 [ZJOI2010]network 网络扩容 【最大流,费用流】

    1834: [ZJOI2010]network 网络扩容 Time Limit: 3 Sec  Memory Limit: 64 MB Submit: 3394  Solved: 1774 [Subm ...

  4. [BZOJ1834][ZJOI2010]network 网络扩容 最大流+费用流

    1834: [ZJOI2010]network 网络扩容 Time Limit: 3 Sec  Memory Limit: 64 MB Submit: 3330  Solved: 1739 [Subm ...

  5. 【最大流】【费用流】bzoj1834 [ZJOI2010]network 网络扩容

    引用题解: 最大流+费用流. 第一问最大流即可. 第二问为“最小费用最大流”. 由题意,这一问的可转化为在上一问的“残量网络”上,扩大一些边的容量,使能从新的图中的最大流为k. 那么易得:对于还有剩余 ...

  6. bzoj1834: [ZJOI2010]network 网络扩容 费用流

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

  7. 【费用流】bzoj1834: [ZJOI2010]network 网络扩容

    还是稍微记一下这个拆点模型吧 Description 给定一张有向图,每条边都有一个容量C和一个扩容费用W.这里扩容费用是指将容量扩大1所需的费用. 求:  1.在不扩容的情况下,1到N的最大流:  ...

  8. bzoj1834 [ZJOI2010]network 网络扩容

    第一问跑最大流,第二问新建一条边连接0和1,流量为上第一问的答案+k,费用为0,接下来图中每条边拆成两条边,第一条容量为C费用为0,第二条容量无穷费用为W,再跑一遍费用流即可. 代码 #include ...

  9. BZOJ 1834: [ZJOI2010]network 网络扩容 最小费用流_最大流_残量网络

    对于第一问,跑一遍最大流即可. 对于第二问,在残量网络上的两点间建立边 <u,v>,容量为无限大,费用为扩充费用. 跑一遍最小费用流即可. Code: #include <vecto ...

随机推荐

  1. 下载bilibili视频

    http://www.urlgot.com/zh_CN/

  2. RAID0+1 RAID5 性能比较

    我想大家都很清楚,存储设备性能的好坏除了与处理器(CPU).缓存等有关之外,还与RAID组中的磁盘数量有很大的关系.按RAID技术的特点,相同磁盘数量下的RAID0性能高于RAID5,RAID1的性能 ...

  3. Spring 集成Hibernate的三种方式

    首先把hibernate的配置文件hibernate.cfg.xml放入spring的src目录下,并且为了便于测试导入了一个实体类Student.java以及它的Student.hbm.xml文件 ...

  4. 2基本类型数组和枚举类型——重拾Java

    2.1 标识符和关键字 2.1.1标识符 标识符:用来标志类名.变量名.方法名.类型名.数组名.文件名的有效字符序列称为标识符.简单地说,标识符就是一个名字. Java关于标识符的语法规则 标识符由字 ...

  5. with as (转)

    sql with as 用法(适用sqlserver,好像oracle也适用) Server 2005中提供了公用表表达式(CTE),使用CTE,可以使SQL语句的可维护性,同时,CTE要比表变量的效 ...

  6. 好用的模板引擎NVelocity

    CastleNVelocity-1.1.1,使用方法: 把dll放到项目中,添加引用,修改配置的文件夹以及数据模型,最后在逻辑代码中调用即可. 封装到CommonHelper.cs using Sys ...

  7. 项目中调试SQLServer 方便的查看SQL语句的执行时间的方法

    第一种方法,先记录执行前的时间,然后在记录执行Sql后的时间,然后做减法 1 第一种方法: 2 declare @begin_date datetime 3 declare @end_date dat ...

  8. 第八章 高级搜索树 (a1)伸展树:逐层伸展

  9. js字符串解析成数字

    parseInt() 先把参数转换成字符串:左边有连续的数字则返回数值,若没有则返回NaN. console.log('parseInt(null)',parseInt(null)); // NaN ...

  10. 41-json.decoder.JSONDecodeError: Invalid control character at: line 6894 column 12 (char 186418)

    在使用python中将单词本的单词用正则匹配成字典后,以json存储,仪json读入,但是一直报错: 原因是: 正则处理后的数据有的出了点问题,导致一个字典的 有多个相同的键!!!,则肯定会报错啊!! ...