[题解] Atcoder ARC 142 E Pairing Wizards 最小割
题目
建图很妙,不会。
考虑每一对要求合法的巫师(x,y),他们两个的\(a\)必须都大于\(min(b_x,b_y)\)。所以在输入的时候,如果\(a_x\)或者\(a_y\)小于\(min(b_x,b_y)\),可以先把\(a_x\)和\(a_y\)提升到\(min(b_x,b_y)\)(以后的a数组都指做过这步操作的)。接下来如果\(max(a_x,a_y)\geq max(b_x,b_y)\),那么这一对已经符合要求,可以直接不管。接下来只考虑需要管的。
发现每一对需要管的(x,y),两个巫师中有且仅有一个满足\(a<b\)。因此可以把所有需要管的(x,y)涉及到的巫师分成左右两组,左边的\(a<b\),右边的\(a \geq b\)。把每一对(x,y)都当成连接x,y的无向边,发现连成了一张二分图。把左边的点称为X类点,右边的称为Y类点。现在我们想让每一个X类点i满足:\(b_i \leq max(a_i,max_{j和i相邻}a_j)\),容易发现这与原题的要求等价。
提升X类点和Y类点的strength分别有什么用呢?发现提升一个X类点i,只对i自己有用,也就是使得\(a_i \geq b_i\),因为它的所有相邻的点j都是Y类点,都已经满足了\(a_j\geq b_j\)。提升一个Y类点则可能会帮助满足相邻的X类点的要求
考虑最小割。对每一个X类点i建立一个点\(c_i\),对于每一个Y类点i和一个[1,100]的整数j建立一个点\(d_{i,j}\)。从源点向每一个\(c_i\)连\(b_i-a_i\)的边,割掉这条边意味着点i的要求靠提升\(a_i\)满足;从每一个\(d_{i,j}向d_{i,j-1}连inf的边\);每一个\(d_{i,j}\)向汇点连边,\(j\leq a_i\)则流量为0,否则为1。最后,对于每一对(x,y),从\(c_x\)向\(d_{y,b_x}\)连inf的边。这个图的最小割加上输入的时候提升的值就是答案。
正确性证明:发现对于每条最后一种类型的边,这条边是不会割的,为了保证没有剩余流量,要么\(s \to c_x\)的边被割掉,代表提升\(a_x\)来满足要求;要么割掉\(d_{y,p}(p \leq b_x)\)到t的所有边(根据建图,这些边的代价和为\(b_x-a_y\)),代表\(a_y\)至少为\(b_x\)。证毕。
时间复杂度:O(能过)。
点击查看代码
#include <bits/stdc++.h>
#define rep(i,n) for(int i=0;i<n;++i)
#define repn(i,n) for(int i=1;i<=n;++i)
#define LL long long
#define pii pair <int,int>
#define pb push_back
#define fi first
#define se second
#define mpr make_pair
using namespace std;
struct node
{
LL to,flow,rev;
node(LL a,LL b,LL c):to(a),flow(b),rev(c){}
};
namespace maxFlow
{
LL n,dist[100010],cur[100010];
vector <node> g[100010];
queue <LL> q;
void add_edge(LL x,LL y,LL z)
{
g[x].pb(node(y,z,g[y].size()));
g[y].pb(node(x,0,g[x].size()-1));
}
void bfs(LL s)
{
rep(i,n+5) dist[i]=-1;
dist[s]=0;q.push(s);
while(!q.empty())
{
int f=q.front();q.pop();
rep(i,g[f].size())
{
if(dist[g[f][i].to]!=-1||g[f][i].flow==0) continue;
dist[g[f][i].to]=dist[f]+1;
q.push(g[f][i].to);
}
}
}
LL dfs(LL pos,LL t,LL flow)
{
if(pos==t) return flow;
for(int i=cur[pos];i<g[pos].size();++i,++cur[pos])
{
if(g[pos][i].flow==0||dist[g[pos][i].to]!=dist[pos]+1) continue;
LL flow2=dfs(g[pos][i].to,t,min(flow,g[pos][i].flow));
if(flow2>0)
{
g[pos][i].flow-=flow2;
g[g[pos][i].to][g[pos][i].rev].flow+=flow2;
return flow2;
}
}
return 0;
}
LL max_flow(LL s,LL t)
{
LL ret=0;
while(true)
{
bfs(s);
if(dist[t]==-1) return ret;
rep(i,n+5) cur[i]=0;
while(true)
{
LL add=dfs(s,t,1e12);
if(add==0) break;
ret+=add;
}
}
}
}
int n,m,a[110],b[110],x[10010],y[10010],ans=0,X[110],Y[110][110],len=0,ss,tt;
bool isx[110],isy[110];
vector <pii> need;
int main()
{
cin>>n;
rep(i,n) scanf("%d%d",&a[i],&b[i]);
cin>>m;
rep(i,m)
{
scanf("%d%d",&x[i],&y[i]);--x[i];--y[i];
int least=min(b[x[i]],b[y[i]]);
if(a[x[i]]<least) ans+=least-a[x[i]],a[x[i]]=least;
if(a[y[i]]<least) ans+=least-a[y[i]],a[y[i]]=least;
}
rep(i,m)
{
if(max(a[x[i]],a[y[i]])>=max(b[x[i]],b[y[i]])&&min(a[x[i]],a[y[i]])>=min(b[x[i]],b[y[i]])) continue;
if(b[y[i]]>a[y[i]]) swap(x[i],y[i]);
need.pb(mpr(x[i],y[i]));
isx[x[i]]=true;isy[y[i]]=true;
}
ss=++len;tt=++len;
rep(i,n)
{
if(isx[i])
{
X[i]=++len;
maxFlow::add_edge(ss,X[i],b[i]-a[i]);
}
else if(isy[i])
{
repn(j,100) Y[i][j]=++len;
repn(j,100)
{
if(j+1<=100) maxFlow::add_edge(Y[i][j+1],Y[i][j],1e12);
maxFlow::add_edge(Y[i][j],tt,(j<=a[i] ? 0:1));
}
}
}
rep(i,need.size()) maxFlow::add_edge(X[need[i].fi],Y[need[i].se][b[need[i].fi]],1e12);
maxFlow::n=len;
ans+=maxFlow::max_flow(ss,tt);
cout<<ans<<endl;
return 0;
}
[题解] Atcoder ARC 142 E Pairing Wizards 最小割的更多相关文章
- [题解] Atcoder ARC 142 D Deterministic Placing 结论,DP
题目 (可能有点长,但是请耐心看完,个人认为比官方题解好懂:P) 首先需要注意,对于任意节点i上的一个棋子,如果在一种走法中它走到了节点j,另一种走法中它走到了节点k,那么这两种走法进行完后,棋子占据 ...
- 【题解】 bzoj3894: 文理分科 (网络流/最小割)
bzoj3894,懒得复制题面,戳我戳我 Solution: 首先这是一个网络流,应该还比较好想,主要就是考虑建图了. 我们来分析下题面,因为一个人要么选文科要么选理科,相当于两条流里面割掉一条(怎么 ...
- AtCoder AGC038F Two Permutations (网络流、最小割)
题目链接 https://atcoder.jp/contests/agc038/tasks/agc038_f 题解 好题. 首先观察到一个性质,对于排列\(P\), 其所形成的每个轮换中的点\(A_i ...
- bzoj2132: 圈地计划(无比强大的最小割)
2132: 圈地计划 题目:传送门 简要题意: 给出一个矩阵,一共n*m个点,并给出三个收益矩阵.A矩阵表示这个点建A的可取收益,B矩阵表示这个点建B的可取收益,C矩阵表示如果相邻(有且仅有一条公共边 ...
- Atcoder Regular Contest 125 E - Snack(最小割转化+贪心)
Preface: 这是生平第一道现场 AC 的 arc E,也生平第一次经历了 performance \(\ge 2800\),甚至还生平第一次被 hb 拉到会议里讲题,讲的就是这个题,然鹅比较尬 ...
- Atcoder Grand Contest 038 F - Two Permutations(集合划分模型+最小割)
洛谷题面传送门 & Atcoder 题面传送门 好久前做的题了--今天偶然想起来要补个题解 首先考虑排列 \(A_i\) 要么等于 \(i\),要么等于 \(P_i\) 这个条件有什么用.我们 ...
- BZOJ-2561-最小生成树 题解(最小割)
2561: 最小生成树(题解) Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1628 Solved: 786 传送门:http://www.lyd ...
- 洛谷P1345 [USACO5.4]奶牛的电信Telecowmunication【最小割】分析+题解代码
洛谷P1345 [USACO5.4]奶牛的电信Telecowmunication[最小割]分析+题解代码 题目描述 农夫约翰的奶牛们喜欢通过电邮保持联系,于是她们建立了一个奶牛电脑网络,以便互相交流. ...
- 【bzoj1001】【最短路】【对偶图】【最大流转最小割】狼抓兔子题解
[BZOJ1001]狼抓兔子 1001: [BeiJing2006]狼抓兔子 Time Limit: 15 Sec Memory Limit: 162 MBSubmit: 18872 Solved ...
随机推荐
- 平衡树——splay 一
splay 一种平衡树,同时也是二叉排序树,与treap不同,它不需要维护堆的性质,它由Daniel Sleator和Robert Tarjan(没错,tarjan,又是他)创造,伸展树是一种自调整二 ...
- 在阿里云Centos7.6上利用docker搭建Jenkins来自动化部署Django项目
原文转载自「刘悦的技术博客」https://v3u.cn/a_id_84 一般情况下,将一个项目部署到生产环境的流程如下: 需求分析-原型设计-开发代码-内网部署-提交测试-确认上线-备份数据-外网更 ...
- Linux、Ubuntu常用命令
# 文件解压缩 # zip压缩目录(附带目录权限) zip -q -r html.zip /home/html 压缩目录 tar -zcvf pack.tar.gz pack/ #打包压缩为一个.gz ...
- SpringMvc 如何同时支持 Jsp 和 Json 接口?
后端同学基本都会见过这种场景:在同一个工程中,有些页面使用jsp模版渲染,同时还有其他接口提供Json格式的返回值.为了同时支持这两种场景,我们一般是如何处理的呢? 其实非常简单: 1.在项目中为 S ...
- Canvas 非线性图形(一):文本
基础 画布除了绘制图形以外还可以绘制文本,画布中的文本可以设置字体大小.字体格式.对齐方式(横向和纵向对齐方式),并且还可以制作很炫酷的文本,比如渐变文字. 文本有以下三个属性,控制文本的字体大小.字 ...
- 分库分表之ShardingSphere
目录 分库分表诞生的前景 分库分表的方式(垂直拆分,水平复制) 1.垂直拆分 1.1 垂直分库 1.2 垂直分表 2.水平拆分 2.1 水平分库 2.2 水平分表 分库分库中间件 ShardingSp ...
- HTML初学者小知识
引用js <script src="链接/js代码位置" type="text/javascript"></script> 引用css ...
- hadoop项目之求出每年二月的最高气温(Combiner优化)
hadoop项目之求出每年二月的最高气温(Combiner优化) 一.项目结构 一.java实现随机生成日期和气温 package com.shujia.weather; import java.io ...
- 播放器之争:VLC还是martPlayer
好多开发者跟我们交流的时候提到,为什么有了VLC这种开源播放器,大牛直播SDK还要开发SmartPlayer?以下就针对VLC和SmartPlayer功能支持和涉及侧重,做个大概的比较: VLC VL ...
- KingbaseES V8R6C5单实例sys_backup.sh备份案例
案例说明: KingbaseES V8R6C5版本中使用了securecmdd工具,用于主机节点间的通讯,默认端口8890.备份工具sys_backup.sh默认使用了securecmdd工具,对 ...