BZOJ 1834--网络扩容(最大流&费用流)
1834: [ZJOI2010]network 网络扩容
Time Limit: 3 Sec Memory Limit: 64 MB
Submit: 3351 Solved: 1750
[Submit][Status][Discuss]
Description
Input
Output
Sample Input
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
30%的数据中,N<=100
100%的数据中,N<=1000,M<=5000,K<=10
题目链接:
http://www.lydsy.com/JudgeOnline/problem.php?id=1834
Solution
值得一做的好题。。。
第一问是赤果果的最大流。。。上模板就好。。。
然后问扩流的最小费用。。。发现可以直接利用之前跑完最大流剩下的残余网络。
在每一条边上额外加一条容量无限大的边,但是这样不能控制流的大小,所以要额外添加一个汇点S。
从S向原起点连一条容量为k的边,再跑一遍费用流,就是第二问的答案。。。
代码
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<iostream>
#define N 1010
#define inf 2000000000
using namespace std;
inline int Read(){
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 n,m,k,ans=0,cnt=1,S,T;
int hed[N],dis[N],l[N],q[200000];
bool vis[N];
struct edge{
int l,r,nxt,v,c,t;
}e[50000];
void insert(int u,int v,int w,int c){
cnt++;e[cnt].l=u;e[cnt].r=v;e[cnt].nxt=hed[u];hed[u]=cnt;e[cnt].v=w;e[cnt].c=0;e[cnt].t=c;
cnt++;e[cnt].l=v;e[cnt].r=u;e[cnt].nxt=hed[v];hed[v]=cnt;e[cnt].v=0;e[cnt].c=0;e[cnt].t=-c;
}
bool BFS(){
int head=0,tail=1,now;
memset(dis,-1,sizeof(dis));
dis[S]=1;q[1]=S;
while(head!=tail){
head++;
now=q[head];
for(int i=hed[now];i;i=e[i].nxt)
if(e[i].v && dis[e[i].r]==-1){
dis[e[i].r]=dis[now]+1;
q[++tail]=e[i].r;
}
}
return dis[T]!=-1;
}
int DFS(int x,int F){
if(x==T)return F;
int w,used=0;
for(int i=hed[x];i;i=e[i].nxt)
if(dis[e[i].r]==dis[x]+1){
w=F-used;
if(w>e[i].v)w=e[i].v;
w=DFS(e[i].r,w);
e[i].v-=w;
e[i^1].v+=w;
used+=w;
if(used==F)return F;
}
if(!used)dis[x]=-1;
return used;
}
void dinic(){
while( BFS() )
ans+=DFS(S,inf);
}
void ins(int u,int v,int w,int c){
cnt++;e[cnt].l=u;e[cnt].r=v;e[cnt].v=w;e[cnt].nxt=hed[u];hed[u]=cnt;e[cnt].c=c;
cnt++;e[cnt].l=v;e[cnt].r=u;e[cnt].v=0;e[cnt].nxt=hed[v];hed[v]=cnt;e[cnt].c=-c;
}
bool spfa(){
int head=0,tail=1,now;
for(int i=0;i<=n;i++)dis[i]=inf;
q[1]=S;dis[S]=0;vis[S]=1;
while(head!=tail){
head++;
now=q[head];
for(int i=hed[now];i;i=e[i].nxt)
if(e[i].v>0 && dis[now]+e[i].c<dis[e[i].r]){
dis[e[i].r]=dis[now]+e[i].c;
l[e[i].r]=i;
if(!vis[e[i].r]){
tail++;
q[tail]=e[i].r;
vis[e[i].r]=1;
}
}
vis[now]=0;
}
return dis[T]!=inf;
}
int mincf(){
int now,used=inf;
now=l[T];
while(now){
if(used>e[now].v)used=e[now].v;
now=l[e[now].l];
}
now=l[T];
while(now){
e[now].v-=used;
e[now^1].v+=used;
ans+=used*e[now].c;
now=l[e[now].l];
}
}
void Cost_Flow(){
while( spfa() )
ans+=mincf();
}
int main(){
int u,v,w,c;
n=Read();m=Read();k=Read();
S=1;T=n;
for(int i=1;i<=m;i++){
u=Read();v=Read();w=Read();c=Read();
insert(u,v,w,c);
}
dinic();
printf("%d ",ans);
ans=0;
int tot=cnt;
for(int i=2;i<=tot;i+=2)
ins(e[i].l,e[i].r,inf,e[i].t);
S=0;
ins(0,1,k,0);
Cost_Flow();
printf("%d\n",ans);
return 0;
}
This passage is made by Iscream-2001.
BZOJ 1834--网络扩容(最大流&费用流)的更多相关文章
- BZOJ 1834 网络扩容(最大流+费用流)
对于第一问,直接求最大流. 对于第二问,建源点s和汇点t,s连1容量为INF,费用为0的边,n连t容量为最大流+k,费用为0的边.这样就把最大流限制为最多增加k了. 限制需要求扩充的最小费用,原图的边 ...
- BZOJ 1834 网络扩容 最大流+最小费用流
题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id=1834 题目大意: 给定一张有向图,每条边都有一个容量C和一个扩容费用W.这里扩容费用是 ...
- BZOJ 1834网络扩容题解
一道不算太难的题目 但是真的很恶心 显然,对于第一问,我们直接无脑打模板就好了 第二问也不是很难,我们将每条边再连一条容量为inf,费用为w的边 但是流量只要小于第一问的答案加k就行了 所以我们增加一 ...
- BZOJ 1834: [ZJOI2010]network 网络扩容(网络流+费用流)
一看就知道是模板题= = ,不说什么了= = PS:回去搞期末了,暑假再来刷题了 CODE: #include<cstdio> #include<iostream> #incl ...
- 洛谷P2604 网络扩容 拆点+费用流
原题链接 这题貌似比较水吧,最简单的拆点,直接上代码了. #include <bits/stdc++.h> using namespace std; #define N 1000 #def ...
- BZOJ 3876 支线剧情 | 有下界费用流
BZOJ 3876 支线剧情 | 有下界费用流 题意 这题题面搞得我看了半天没看懂--是这样的,原题中的"剧情"指的是边,"剧情点"指的才是点. 题面翻译过来大 ...
- 从多种角度看[BZOJ 1061] [NOI 2008]志愿者招募(费用流)
从多种角度看[BZOJ 1061] [NOI 2008]志愿者招募(费用流) 题面 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难题:为即将启动的奥运 ...
- 【BZOJ】1834: [ZJOI2010]network 网络扩容(最大流+费用流)
http://www.lydsy.com/JudgeOnline/problem.php?id=1834 我又思考人生了T_T,nd的数组开小了,一直wa,调了一个小时才发现啊!!!!!我一直以为我的 ...
- bzoj 1834: [ZJOI2010]network 网络扩容 -- 最大流+费用流
1834: [ZJOI2010]network 网络扩容 Time Limit: 3 Sec Memory Limit: 64 MB Description 给定一张有向图,每条边都有一个容量C和一 ...
随机推荐
- 迷你MVVM框架 avalonjs 0.84发布
本版本只要做了如下改进 重构ui绑定 fix scanTag bug ms-include 的值必须不为空值,否则不做任何操作.
- kafka启动报java.net.UnknownHostException
kafka启动报java.net.UnknownHostException 参考资料: 1.https://blog.csdn.net/zdxiq000/article/details/6258765 ...
- filebeat 笔记
认识Beats Beats是用于单用途数据托运人的平台.它们以轻量级代理的形式安装,并将来自成百上千台机器的数据发送到Logstash或Elasticsearch. (画外音:通俗地理解,就是采集数据 ...
- react-native 打包 出apk
先上步骤: 一. 生成签名文件(my-release-key.keystore文件) Android要求所有应用都有一个数字签名才会被允许安装在用户手机上 1. 在项目目录下运行如下命令: keyt ...
- Java项目的结构
-------siwuxie095 以 Hello World 为例 这个工程用一个文件夹表示,被放置在左侧的资源管理面板 Package Explorer 中 ...
- httpwebrequest详解
HttpWebRequest和HttpWebResponse类是用于发送和接收HTTP数据的最好选择.它们支持一系列有用的属性.这两个类位 于System.Net命名空间,默认情况下这个类对于控制台程 ...
- springMVC入门程序。使用springmvc实现商品列表的展示。
1.1 开发环境 本教程使用环境: Jdk:jdk1.7.0_72 Eclipse:mars Tomcat:apache-tomcat-7.0.53 Springmvc:4.1.3 1.2 需求 使用 ...
- DBA的工作职责和每日工作
DBA一般职责 1.安装和升级数据库服务器,以及应用程序工具构建和配置网络环境. 2.熟悉数据库系统的存储结构预测未来的存储需求,制订数据库的存储方案. 3.根据开发人员设计的应用系统需求创建数据库存 ...
- <<高级计算机网络>>(Advaned Computer Networks) 徐恪 徐明伟 陈文龙 马东超
目录 第1章 计算机网络与Internet1 1.1 引言1 1.2 Internet发展历史2 1.2.1 互联网发展的主要阶段4 1.2.2 互联网在中国的发展5 1.2.3 互联网主要创新5 1 ...
- 自动创建orcl表
using System;using System.Collections.Generic;using System.Data;using System.Linq;using System.Text; ...