BZOJ_1834_[ZJOI2010]network 网络扩容_费用流
BZOJ_1834_[ZJOI2010]network 网络扩容_费用流
题意:
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <queue>
using namespace std;
#define N 2050
#define M 60050
#define inf 100000000
#define S (n+1)
#define T (n+2)
int head[N],to[M],nxt[M],flow[M],val[M],path[M],cnt=1;
int n,m,k,Q[N],l,r,dep[N],inq[N];
int xx[M],yy[M],cc[M],ww[M],dis[N];
inline void add1(int u,int v,int f)
{
to[++cnt]=v;nxt[cnt]=head[u];head[u]=cnt;flow[cnt]=f;
}
inline void add2(int u,int v,int f,int c)
{
to[++cnt]=v;nxt[cnt]=head[u];head[u]=cnt;flow[cnt]=f;val[cnt]=c;
}
bool bfs()
{
memset(dep,0,sizeof(dep));
l=r=0;
dep[1]=1;
Q[r++]=1;
int i;
while(l^r)
{
int x=Q[l++];if(l==n+1)l=0;
for(i=head[x];i;i=nxt[i])if(!dep[to[i]]&&flow[i])
{
dep[to[i]]=dep[x]+1;
if(to[i]==n)return 1;
Q[r++]=to[i];if(r==n+1)r=0;
}
}
return 0;
}
int dfs(int x,int mf)
{
if(x==n)return mf;
int nf=0;
int i;
for(i=head[x];i;i=nxt[i])if(dep[to[i]]==dep[x]+1&&flow[i])
{
int tmp=dfs(to[i],min(flow[i],mf-nf));
nf+=tmp;
flow[i]-=tmp;
flow[i^1]+=tmp;
if(nf==mf)break;
}
dep[x]=0;
return nf;
}
bool spfa()
{
memset(dis,0x3f,sizeof(dis));
memset(path,0,sizeof(path));
l=r=0;
int i;
Q[r++]=S;dis[S]=0;inq[S]=1;
while(l^r)
{
int x=Q[l++];if(l==n+4)l=0;inq[x]=0;
for(i=head[x];i;i=nxt[i])
{
if(flow[i]>0&&dis[to[i]]>dis[x]+val[i])
{
dis[to[i]]=dis[x]+val[i];
path[to[i]]=i^1;
if(!inq[to[i]])
{
inq[to[i]]=1;
Q[r++]=to[i];
if(r==n+4)r=0;
}
}
}
}
return dis[T]<inf;
}
void mcmf()
{
int mnc=0,mxf=0,i;
while(spfa())
{
int nf=1<<30;
for(i=T;i!=S;i=to[path[i]])
{
nf=min(nf,flow[path[i]^1]);
}
for(i=T;i!=S;i=to[path[i]])
{
flow[path[i]]+=nf;
flow[path[i]^1]-=nf;
mnc+=val[path[i]^1]*nf;
}
mxf+=nf;
}
printf("%d\n",mnc);
}
void dinic()
{
int ans=0,f;
while(bfs())
{
while(f=dfs(1,inf))
ans+=f;
}
k+=ans;
printf("%d ",ans);
}
int main()
{
scanf("%d%d%d",&n,&m,&k);
int i;
for(i=1;i<=m;i++)
{
scanf("%d%d%d%d",&xx[i],&yy[i],&cc[i],&ww[i]);
add1(xx[i],yy[i],cc[i]);
add1(yy[i],xx[i],0);
}
dinic();
cnt=1;
memset(head,0,sizeof(head));
add2(S,1,k,0);
add2(1,S,0,0);
add2(n,T,k,0);
add2(T,n,0,0);
for(i=1;i<=m;i++)
{
add2(xx[i],yy[i],cc[i],0);
add2(yy[i],xx[i],0,0);
add2(xx[i],yy[i],inf,ww[i]);
add2(yy[i],xx[i],0,-ww[i]);
}
mcmf();
}
BZOJ_1834_[ZJOI2010]network 网络扩容_费用流的更多相关文章
- BZOJ 1834 [ZJOI2010]network 网络扩容(费用流)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=1834 [题目大意] 给定一张有向图,每条边都有一个容量C和一个扩容费用W. 这里扩容费 ...
- BZOJ1834 [ZJOI2010]network 网络扩容 【最大流,费用流】
1834: [ZJOI2010]network 网络扩容 Time Limit: 3 Sec Memory Limit: 64 MB Submit: 3394 Solved: 1774 [Subm ...
- 【BZOJ】1834: [ZJOI2010]network 网络扩容(最大流+费用流)
http://www.lydsy.com/JudgeOnline/problem.php?id=1834 我又思考人生了T_T,nd的数组开小了,一直wa,调了一个小时才发现啊!!!!!我一直以为我的 ...
- 2018.10.13 bzoj1834: [ZJOI2010]network 网络扩容(最大流+费用流)
传送门 网络流水题啊. 第一问直接放心跑最大流(本来还以为有什么tricktricktrick). 第二问就直接把原来的边(u,v,c,w)(u,v,c,w)(u,v,c,w)变成(u,v,c,0)( ...
- bzoj 1834: [ZJOI2010]network 网络扩容【最大流+最小费用最大流】
第一问直接跑最大流即可.建图的时候按照费用流建,费用为0. 对于第二问,在第一问dinic剩下的残量网络上建图,对原图的每条边(i,j),建(i,j,inf,cij),表示可以用c的花费增广这条路.然 ...
- BZOJ 1834: [ZJOI2010]network 网络扩容(最大流+最小费用最大流)
第一问直接跑最大流.然后将所有边再加一次,费用为扩容费用,容量为k,再从一个超级源点连一条容量为k,费用为0的边到原源点,从原汇点连一条同样的边到超级汇点,然 后跑最小费用最大流就OK了. ---- ...
- bzoj 1834: [ZJOI2010]network 网络扩容 -- 最大流+费用流
1834: [ZJOI2010]network 网络扩容 Time Limit: 3 Sec Memory Limit: 64 MB Description 给定一张有向图,每条边都有一个容量C和一 ...
- [BZOJ1834][ZJOI2010]network 网络扩容 最大流+费用流
1834: [ZJOI2010]network 网络扩容 Time Limit: 3 Sec Memory Limit: 64 MB Submit: 3330 Solved: 1739 [Subm ...
- 【BZOJ1834】[ZJOI2010]network 网络扩容 最大流+最小费用流
[BZOJ1834][ZJOI2010]network 网络扩容 Description 给定一张有向图,每条边都有一个容量C和一个扩容费用W.这里扩容费用是指将容量扩大1所需的费用.求: 1. 在不 ...
随机推荐
- 如何通过jQuery获取一个没有定高度的元素---------的自适应高度(offsetHeight的正确使用方法)
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/ ...
- 洛谷P1919 【模板】A*B Problem升级版 题解(FFT的第一次实战)
洛谷P1919 [模板]A*B Problem升级版(FFT快速傅里叶) 刚学了FFT,我们来刷一道模板题. 题目描述 给定两个长度为 n 的两个十进制数,求它们的乘积. n<=100000 如 ...
- 窗口函数解决数据岛问题(mysql暂无窗口函数,可以通过用户变量解决窗口函数问题)
数据岛问题: 有表: create table dataisland (id int) insert into dataisland values(1),(2),(3),(7),(11),(12) ...
- JavaScript头像上传器的实现
最近做这方面的东西,刚开始准备用一个开源项目:https://github.com/yueyoum/django-upload-avatar 后来发现这个开源组件的原设计者的定制化选项设计略显复杂,发 ...
- linux 系统centOS 7 怎么在线安装mysql
以下操作,注意空格,为避免出错,可以选择直接copy 第一步: 输入命令行 yum list installed | grep mysql 检查是否已经安装mysql 已安装输入 ...
- MFC学习问题总结
1.学习MFC添加位图,无法获取其ID 1).点击视图->其他窗口->资源视图,你会发现“无法在此窗口显示”,找到resource.h文件,关闭即可重新走一遍上面的过程就会发现可以打开了. ...
- JS实现鼠标放在文字上面显示全部内容
web中当我们把text等的宽固定后如果文本框中内容过多就只能看到前面部分的内容,这时我们可以用样式控制当鼠标移到文本框时显示全部内容. var pointX; var pointY; $(funct ...
- java面试总结
一.java的集合框架 HashMap.HashTable.CurrentHashMap的底层数据结构与区别? CurrentHashMap与HashTable是如何保证线程安全的? ArrayLis ...
- XShell上传文件到Linux服务器上
在学习Linux过程中,我们常常需要将本地文件上传到Linux主机上,这里简单记录下使用Xsheel工具进行文件传输 1:首先连接上一台Linux主机 2:输入rz命令,看是否已经安装了lrzsz,如 ...
- Python_ jiba、snownlp中文分词、pypinyin中文转拼音
import jieba #导入jieba模块 x = '分词的准确度直接影响了后续文本处理和挖掘算法的最终效果.' jieba.cut(x) #使用默认词库进行分词 print(list(jieba ...