传送门

费用流经典题。


按照题目要求建边。

为了方便我将所有格子拆点,三种情况下容量分别为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 题」数字梯形(费用流)的更多相关文章

  1. 2018.10.14 loj#6012. 「网络流 24 题」分配问题(费用流)

    传送门 费用流水题. 依然是照着题意模拟建边就行了. 为了练板子又重新写了一遍费用流. 代码: #include<bits/stdc++.h> #define N 305 #define ...

  2. 2018.10.14 loj#6011. 「网络流 24 题」运输问题(费用流)

    传送门 费用流入门题. 直接按照题意模拟. 把货物的数量当做容量建边. 然后跑一次最小费用流和最大费用流就行了. 代码: #include<bits/stdc++.h> #define N ...

  3. LOJ #6010. 「网络流 24 题」数字梯形

    #6010. 「网络流 24 题」数字梯形   题目描述 给定一个由 n nn 行数字组成的数字梯形如下图所示.梯形的第一行有 m mm 个数字.从梯形的顶部的 m mm 个数字开始,在每个数字处可以 ...

  4. 【刷题】LOJ 6010 「网络流 24 题」数字梯形

    题目描述 给定一个由 \(n\) 行数字组成的数字梯形如下图所示.梯形的第一行有 \(m\) 个数字.从梯形的顶部的 \(m\) 个数字开始,在每个数字处可以沿左下或右下方向移动,形成一条从梯形的顶至 ...

  5. 2018.10.15 loj#6013. 「网络流 24 题」负载平衡(费用流)

    传送门 费用流sb题. 直接从sss向每个点连边,容量为现有物品量. 然后从ttt向每个点连边,容量为最后库存量. 由于两个点之间可以互相任意运送物品,因此相邻的直接连infinfinf的边就行了. ...

  6. Libre 6010「网络流 24 题」数字梯形 (网络流,最大费用最大流)

    Libre 6010「网络流 24 题」数字梯形 (网络流,最大费用最大流) Description 给定一个由n 行数字组成的数字梯形如下图所示.梯形的第一行有m 个数字.从梯形的顶部的m 个数字开 ...

  7. 2018.10.14 loj#6003. 「网络流 24 题」魔术球(最大流)

    传送门 网络流好题. 这道题可以动态建图. 不难想到把每个球iii都拆点成i1i_1i1​和i2i_2i2​,每次连边(s,i1),(i2,t)(s,i_1),(i_2,t)(s,i1​),(i2​, ...

  8. 【刷题】LOJ 6227 「网络流 24 题」最长k可重线段集问题

    题目描述 给定平面 \(\text{xoy}\) 上 \(n\) 个开线段组成的集合 \(\text{I}\) ,和一个正整数 \(k\) ,试设计一个算法. 从开线段集合 \(\text{I}\) ...

  9. [luogu_P1251][LOJ#6008]「网络流 24 题」餐巾计划

    [luogu_P1251][LOJ#6008]「网络流 24 题」餐巾计划 试题描述 一个餐厅在相继的 \(N\) 天里,第 \(i\) 天需要 \(R_i\) 块餐巾 \((i=l,2,-,N)\) ...

随机推荐

  1. HTML5 借助http请求发送formdata对象,从而上传文件 XMLHttpRequest, FormData

    <!DOCTYPE html> <html xmlns="http://www.w3.org/1999/xhtml"> <head> <m ...

  2. python之model模块和包的介绍

    一,模块的概念:常见场景:一个模块就是一个包含了一组功能的Python文件,比如spam.py,模块名为spam,可以通过import spam使用 在Python中,模块的使用方式都是一样的,但其实 ...

  3. JSTL标签库学习记录2-fmt

    fmt的标签为辅助性功能标签 设置编码 <fmt:requestEncoding value=""> 国际化相关 <fmt:setLocale value=&qu ...

  4. mysql之MMM高可用方案简介

    MMM(Master-Master Replication Manager for MySQL) MySQL主主复制管理器,提供了MySQL主主复制配置的监控.故障转移和管理的一套可伸缩的脚本套件.在 ...

  5. mtr 命令详解

    一般在windows 来判断网络连通性用ping 和tracert,ping的话可以来判断丢包率,tracert可以用来跟踪路由,在Linux中有一个更好的网络连通性判断工具,它可以结合ping ns ...

  6. Linux下方便的块设备查看工具lsblk

    之前在Linux下看有什么块设备,通常都用fdisk什么的或者直接ls /dev/ 去看很不方便. 这个工具属于util-linux-ng包,在RHEL 6.1上是安装好的啦,直接用就好. ubunt ...

  7. Haskell语言学习笔记(66)Aeson

    Data.Aeson 安装 aeson $ cabal install aeson Installed aeson-1.2.3.0 Prelude> :m +Data.Aeson Prelude ...

  8. spring 整合 struts2 xml配置

    整合之前要搞清楚struts2是什么; struts2:表现层框架  增删改查  作用域  页面跳转   异常处理  ajax 上传下载  excel   调用service spring :IOC/ ...

  9. 疯狂JAVA——第五章 面向对象(上)

    5.1类和对象 构造器是一个类创建对象的根本途径,如果一个类没有构造器,这个类通常无法创建实例.通过new关键字来调用构造器,从而返回该类的实例. 类名:每个单词首字母大写,其他字母小写,单词之间不要 ...

  10. 迭代删除元素 并发bug 低级错误

    方法一:HashMap<String, Integer> myHashMap = new HashMap<>();myHashMap.put("1", 1) ...