2018.10.15 loj#6010. 「网络流 24 题」数字梯形(费用流)
传送门
费用流经典题。
按照题目要求建边。
为了方便我将所有格子拆点,三种情况下容量分别为111,infinfinf,infinfinf,费用都为validi,jval_{id_{i,j}}validi,j。
然后从源点向第一排的mmm个点连边,三种情况下容量都为111,费用都为0。
然后从最后一排的m+n−1m+n-1m+n−1个点向汇点连边,三种情况下容量为111,infinfinf,infinfinf,费用都为0。
至于格子之间的路径,三种情况下容量为111,111,infinfinf,费用都为0。
跑三次最大费用流就行了。
代码:
#include<bits/stdc++.h>
#define N 1005
#define M 20005
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 n,m,id[25][25],a[25][25],tot=0,all=0;
struct edge{int v,next,c,w;};
struct MCMF{
edge e[M<<1];
bool in[N];
int first[N],d[N],flow[N],pos[N],pred[N],cnt,s,t;
inline void init(){memset(first,-1,sizeof(first)),cnt=-1,s=0,t=all*2+1;}
inline void addedge(int u,int v,int c,int w){e[++cnt].v=v,e[cnt].w=w,e[cnt].c=c,e[cnt].next=first[u],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 bool spfa(){
queue<int>q;
for(int i=0;i<=t;++i)d[i]=-0x3f3f3f3f;
in[s]=1,d[s]=0,pred[t]=-1,flow[s]=0x3f3f3f3f,q.push(s);
while(!q.empty()){
int x=q.front();
q.pop(),in[x]=0;
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,flow[v]=min(flow[x],e[i].c),pred[v]=x,pos[v]=i;
if(!in[v])in[v]=1,q.push(v);
}
}
}
return d[t]!=-0x3f3f3f3f;
}
inline int solve(){
int ret=0;
for(int w=t;spfa();w=t){
ret+=d[t]*flow[t];
while(w!=s)e[pos[w]].c-=flow[t],e[pos[w]^1].c+=flow[t],w=pred[w];
}
return ret;
}
}mcmf;
int main(){
m=read(),n=read(),all=(m*2+n-1)*n/2;
for(int i=1;i<=n;++i)for(int j=1;j<=m+i-1;++j)id[i][j]=++tot,a[i][j]=read();
mcmf.init();
for(int i=1;i<=m;++i)mcmf.add(mcmf.s,id[1][i],1,0);
for(int i=1;i<n;++i)for(int j=1;j<=m+i-1;++j){
mcmf.add(id[i][j],id[i][j]+all,1,a[i][j]);
mcmf.add(id[i][j]+all,id[i+1][j],1,0);
mcmf.add(id[i][j]+all,id[i+1][j+1],1,0);
}
for(int i=1;i<=m+n-1;++i)mcmf.add(id[n][i],id[n][i]+all,1,a[n][i]),mcmf.add(id[n][i]+all,mcmf.t,1,0);
printf("%d\n",mcmf.solve());
mcmf.init();
for(int i=1;i<=m;++i)mcmf.add(mcmf.s,id[1][i],1,0);
for(int i=1;i<n;++i)for(int j=1;j<=m+i-1;++j){
mcmf.add(id[i][j],id[i][j]+all,0x3f3f3f3f,a[i][j]);
mcmf.add(id[i][j]+all,id[i+1][j],1,0);
mcmf.add(id[i][j]+all,id[i+1][j+1],1,0);
}
for(int i=1;i<=m+n-1;++i)mcmf.add(id[n][i],id[n][i]+all,0x3f3f3f3f,a[n][i]),mcmf.add(id[n][i]+all,mcmf.t,0x3f3f3f3f,0);
printf("%d\n",mcmf.solve());
mcmf.init();
for(int i=1;i<=m;++i)mcmf.add(mcmf.s,id[1][i],1,0);
for(int i=1;i<n;++i)for(int j=1;j<=m+i-1;++j){
mcmf.add(id[i][j],id[i][j]+all,0x3f3f3f3f,a[i][j]);
mcmf.add(id[i][j]+all,id[i+1][j],0x3f3f3f3f,0);
mcmf.add(id[i][j]+all,id[i+1][j+1],0x3f3f3f3f,0);
}
for(int i=1;i<=m+n-1;++i)mcmf.add(id[n][i],id[n][i]+all,0x3f3f3f3f,a[n][i]),mcmf.add(id[n][i]+all,mcmf.t,0x3f3f3f3f,0);
printf("%d",mcmf.solve());
return 0;
}
2018.10.15 loj#6010. 「网络流 24 题」数字梯形(费用流)的更多相关文章
- 2018.10.14 loj#6012. 「网络流 24 题」分配问题(费用流)
传送门 费用流水题. 依然是照着题意模拟建边就行了. 为了练板子又重新写了一遍费用流. 代码: #include<bits/stdc++.h> #define N 305 #define ...
- 2018.10.14 loj#6011. 「网络流 24 题」运输问题(费用流)
传送门 费用流入门题. 直接按照题意模拟. 把货物的数量当做容量建边. 然后跑一次最小费用流和最大费用流就行了. 代码: #include<bits/stdc++.h> #define N ...
- LOJ #6010. 「网络流 24 题」数字梯形
#6010. 「网络流 24 题」数字梯形 题目描述 给定一个由 n nn 行数字组成的数字梯形如下图所示.梯形的第一行有 m mm 个数字.从梯形的顶部的 m mm 个数字开始,在每个数字处可以 ...
- 【刷题】LOJ 6010 「网络流 24 题」数字梯形
题目描述 给定一个由 \(n\) 行数字组成的数字梯形如下图所示.梯形的第一行有 \(m\) 个数字.从梯形的顶部的 \(m\) 个数字开始,在每个数字处可以沿左下或右下方向移动,形成一条从梯形的顶至 ...
- 2018.10.15 loj#6013. 「网络流 24 题」负载平衡(费用流)
传送门 费用流sb题. 直接从sss向每个点连边,容量为现有物品量. 然后从ttt向每个点连边,容量为最后库存量. 由于两个点之间可以互相任意运送物品,因此相邻的直接连infinfinf的边就行了. ...
- Libre 6010「网络流 24 题」数字梯形 (网络流,最大费用最大流)
Libre 6010「网络流 24 题」数字梯形 (网络流,最大费用最大流) Description 给定一个由n 行数字组成的数字梯形如下图所示.梯形的第一行有m 个数字.从梯形的顶部的m 个数字开 ...
- 2018.10.14 loj#6003. 「网络流 24 题」魔术球(最大流)
传送门 网络流好题. 这道题可以动态建图. 不难想到把每个球iii都拆点成i1i_1i1和i2i_2i2,每次连边(s,i1),(i2,t)(s,i_1),(i_2,t)(s,i1),(i2, ...
- 【刷题】LOJ 6227 「网络流 24 题」最长k可重线段集问题
题目描述 给定平面 \(\text{xoy}\) 上 \(n\) 个开线段组成的集合 \(\text{I}\) ,和一个正整数 \(k\) ,试设计一个算法. 从开线段集合 \(\text{I}\) ...
- [luogu_P1251][LOJ#6008]「网络流 24 题」餐巾计划
[luogu_P1251][LOJ#6008]「网络流 24 题」餐巾计划 试题描述 一个餐厅在相继的 \(N\) 天里,第 \(i\) 天需要 \(R_i\) 块餐巾 \((i=l,2,-,N)\) ...
随机推荐
- HTML5 借助http请求发送formdata对象,从而上传文件 XMLHttpRequest, FormData
<!DOCTYPE html> <html xmlns="http://www.w3.org/1999/xhtml"> <head> <m ...
- python之model模块和包的介绍
一,模块的概念:常见场景:一个模块就是一个包含了一组功能的Python文件,比如spam.py,模块名为spam,可以通过import spam使用 在Python中,模块的使用方式都是一样的,但其实 ...
- JSTL标签库学习记录2-fmt
fmt的标签为辅助性功能标签 设置编码 <fmt:requestEncoding value=""> 国际化相关 <fmt:setLocale value=&qu ...
- mysql之MMM高可用方案简介
MMM(Master-Master Replication Manager for MySQL) MySQL主主复制管理器,提供了MySQL主主复制配置的监控.故障转移和管理的一套可伸缩的脚本套件.在 ...
- mtr 命令详解
一般在windows 来判断网络连通性用ping 和tracert,ping的话可以来判断丢包率,tracert可以用来跟踪路由,在Linux中有一个更好的网络连通性判断工具,它可以结合ping ns ...
- Linux下方便的块设备查看工具lsblk
之前在Linux下看有什么块设备,通常都用fdisk什么的或者直接ls /dev/ 去看很不方便. 这个工具属于util-linux-ng包,在RHEL 6.1上是安装好的啦,直接用就好. ubunt ...
- Haskell语言学习笔记(66)Aeson
Data.Aeson 安装 aeson $ cabal install aeson Installed aeson-1.2.3.0 Prelude> :m +Data.Aeson Prelude ...
- spring 整合 struts2 xml配置
整合之前要搞清楚struts2是什么; struts2:表现层框架 增删改查 作用域 页面跳转 异常处理 ajax 上传下载 excel 调用service spring :IOC/ ...
- 疯狂JAVA——第五章 面向对象(上)
5.1类和对象 构造器是一个类创建对象的根本途径,如果一个类没有构造器,这个类通常无法创建实例.通过new关键字来调用构造器,从而返回该类的实例. 类名:每个单词首字母大写,其他字母小写,单词之间不要 ...
- 迭代删除元素 并发bug 低级错误
方法一:HashMap<String, Integer> myHashMap = new HashMap<>();myHashMap.put("1", 1) ...