AtCoder AGC038F Two Permutations (网络流、最小割)
题目链接
https://atcoder.jp/contests/agc038/tasks/agc038_f
题解
好题。
首先观察到一个性质,对于排列\(P\), 其所形成的每个轮换中的点\(A_i\)是选\(i\)还是选\(P_i\)的状态必须相同。\(Q_i\)同理。
然后转化成最小化\(A_i=B_i\)的位置\(i\)数量。
考虑\(A_i=B_i\)的条件:
(1) \(P_i=Q_i=i\), 则此位置无用,\(A_i=B_i\)一定满足。
(2) \(P_i=i, Q_i\ne i\), 则\(A_i=B_i\)等价于\(B_i=i\).
(3) \(P_i\ne i, Q_i=i\), 则\(A_i=B_i\)等价于\(A_i=i\).
(4) \(P_i=Q_i\ne i\), 则\(A_i=B_i\)等价于\(A_i\)和\(B_i\)选的状态(即是\(i\)还是\(P_i\)或\(Q_i\))相同。
(5) \(P_i\ne Q_i\ne i\), 则\(A_i=B_i\)等价于\(A_i=i\)且\(B_i=i\).
那么我们可以从中观察到一个集合划分模型: 把所有点划分为\(S,T\)两个集合,设\(A_i=i\)表示\(A_i\)在\(S\)集,\(A_i=P_i\)表示\(A_i\)在\(T\)集,\(B_i=i\)表示\(i\)在\(T\)集,\(B_i=Q_i\)表示\(i\)在\(S\)集。那么上述条件就可以转化为:
对每个位置\(i\):
(1) \(P_i=Q_i\), 则此位置无用,一定要花费\(1\)的代价。
(2) \(P_i=i, Q_i\ne i\), 则如果\(B_i\)在\(T\)集需要花费\(1\)的代价。(\(S\)连到\(Q_i\))
(3) \(P_i\ne i, Q_i=i\), 则如果\(A_i\)在\(S\)集需要花费\(1\)的代价。(\(P_i\)连到\(T\))
(4) \(P_i=Q_i\ne i\), 则如果\(A_i\)和\(B_i\)在不同的集合要花费\(1\)的代价。(\(P_i\)与\(Q_i\)连双向)
(5) \(P_i\ne Q_i\ne i\), 则如果\(A_i\)在\(S\)集且\(B_i\)在\(T\)集需要花费\(1\)的代价。(\(P_i\)连到\(Q_i\))
于是,给每个轮换建立一个点,再对每个\(i\)分类讨论,将其所对应的轮换连边即可。
由于建出来的图是二分图,且边权都是\(1\), 故时间复杂度\(O(n\sqrt n)\).
我写完之后一直TLE, 最后发现居然是我写了两年的当前弧优化一直是个假的……挂了38发真的是很无语……
代码
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cstring>
using namespace std;
const int INF = 1e7;
namespace NetFlow
{
const int N = 2e5+2;
const int M = 4e5;
struct Edge
{
int v,w,nxt,rev;
} e[(M<<1)+3];
int fe[N+3];
int te[N+3];
int dep[N+3];
int que[N+3];
int n,en,s,t;
void addedge(int u,int v,int w)
{
// printf("addedge %d %d %d\n",u,v,w);
en++; e[en].v = v; e[en].w = w;
e[en].nxt = fe[u]; fe[u] = en; e[en].rev = en+1;
en++; e[en].v = u; e[en].w = 0;
e[en].nxt = fe[v]; fe[v] = en; e[en].rev = en-1;
}
bool bfs()
{
for(int i=1; i<=n; i++) dep[i] = 0;
int head = 1,tail = 1; que[1] = s; dep[s] = 1;
while(head<=tail)
{
int u = que[head]; head++;
for(int i=fe[u]; i; i=e[i].nxt)
{
int v = e[i].v;
if(e[i].w>0 && dep[v]==0)
{
dep[v] = dep[u]+1;
if(v==t)return true;
tail++; que[tail] = v;
}
}
}
return false;
}
int dfs(int u,int cur)
{
if(u==t||cur==0) {return cur;}
int rst = cur;
for(int &i=te[u]; i; i=e[i].nxt)
{
int v = e[i].v;
if(e[i].w>0 && rst>0 && dep[v]==dep[u]+1)
{
int flow = dfs(v,min(rst,e[i].w));
if(flow>0)
{
e[i].w -= flow;
rst -= flow;
e[e[i].rev].w += flow;
if(rst==0) {return cur;}
}
}
}
if(rst==cur) {dep[u] = -2;}
return cur-rst;
}
int dinic(int _n,int _s,int _t)
{
n = _n,s = _s,t = _t;
int ret = 0;
while(bfs())
{
for(int i=1; i<=n; i++) te[i] = fe[i];
memcpy(te,fe,sizeof(int)*(n+1));
ret += dfs(s,INF);
}
return ret;
}
}
using NetFlow::addedge;
using NetFlow::dinic;
const int N = 1e5;
int a[N+3],b[N+3];
int ca[N+3],cb[N+3];
int n,cnta,cntb;
int main()
{
scanf("%d",&n);
for(int i=1; i<=n; i++) scanf("%d",&a[i]),a[i]++;
for(int i=1; i<=n; i++) scanf("%d",&b[i]),b[i]++;
for(int i=1; i<=n; i++)
{
if(!ca[i])
{
cnta++;
int x = i;
do
{
ca[x] = cnta;
x = a[x];
} while(x!=i);
}
}
for(int i=1; i<=n; i++)
{
if(!cb[i])
{
cntb++;
int x = i;
do
{
cb[x] = cntb;
x = b[x];
} while(x!=i);
}
}
// printf("ca: "); for(int i=1; i<=n; i++) printf("%d ",ca[i]); puts("");
// printf("cb: "); for(int i=1; i<=n; i++) printf("%d ",cb[i]); puts("");
int ans = n;
for(int i=1; i<=n; i++)
{
if(a[i]==i && b[i]==i) {ans--;}
else if(a[i]==i) {addedge(1,cb[i]+cnta+2,1);}
else if(b[i]==i) {addedge(ca[i]+2,2,1);}
else
{
addedge(ca[i]+2,cb[i]+cnta+2,1);
if(a[i]==b[i]) {addedge(cb[i]+cnta+2,ca[i]+2,1);}
}
}
ans -= dinic(2+cnta+cntb,1,2);
printf("%d\n",ans);
return 0;
}
AtCoder AGC038F Two Permutations (网络流、最小割)的更多相关文章
- 【题解】 bzoj3894: 文理分科 (网络流/最小割)
bzoj3894,懒得复制题面,戳我戳我 Solution: 首先这是一个网络流,应该还比较好想,主要就是考虑建图了. 我们来分析下题面,因为一个人要么选文科要么选理科,相当于两条流里面割掉一条(怎么 ...
- 【bzoj3774】最优选择 网络流最小割
题目描述 小N手上有一个N*M的方格图,控制某一个点要付出Aij的代价,然后某个点如果被控制了,或者他周围的所有点(上下左右)都被控制了,那么他就算是被选择了的.一个点如果被选择了,那么可以得到Bij ...
- 【bzoj1143】[CTSC2008]祭祀river Floyd+网络流最小割
题目描述 在遥远的东方,有一个神秘的民族,自称Y族.他们世代居住在水面上,奉龙王为神.每逢重大庆典, Y族都会在水面上举办盛大的祭祀活动.我们可以把Y族居住地水系看成一个由岔口和河道组成的网络.每条河 ...
- 【bzoj1797】[Ahoi2009]Mincut 最小割 网络流最小割+Tarjan
题目描述 给定一张图,对于每一条边询问:(1)是否存在割断该边的s-t最小割 (2)是否所有s-t最小割都割断该边 输入 第一行有4个正整数,依次为N,M,s和t.第2行到第(M+1)行每行3个正 整 ...
- 【bzoj1976】[BeiJing2010组队]能量魔方 Cube 网络流最小割
题目描述 一个n*n*n的立方体,每个位置为0或1.有些位置已经确定,还有一些需要待填入.问最后可以得到的 相邻且填入的数不同的点对 的数目最大. 输入 第一行包含一个数N,表示魔方的大小. 接下来 ...
- 【bzoj4177】Mike的农场 网络流最小割
题目描述 Mike有一个农场,这个农场n个牲畜围栏,现在他想在每个牲畜围栏中养一只动物,每只动物可以是牛或羊,并且每个牲畜围栏中的饲养条件都不同,其中第i个牲畜围栏中的动物长大后,每只牛可以卖a[i] ...
- 【bzoj3438】小M的作物 网络流最小割
原文地址:http://www.cnblogs.com/GXZlegend/p/6801522.html 题目描述 小M在MC里开辟了两块巨大的耕地A和B(你可以认为容量是无穷),现在,小P有n中作物 ...
- 【bzoj3144】[Hnoi2013]切糕 网络流最小割
题目描述 输入 第一行是三个正整数P,Q,R,表示切糕的长P. 宽Q.高R.第二行有一个非负整数D,表示光滑性要求.接下来是R个P行Q列的矩阵,第z个 矩阵的第x行第y列是v(x,y,z) (1≤x≤ ...
- 【bzoj3894】文理分科 网络流最小割
原文地址:http://www.cnblogs.com/GXZlegend 题目描述 文理分科是一件很纠结的事情!(虽然看到这个题目的人肯定都没有纠结过) 小P所在的班级要进行文理分科.他的班级可以用 ...
随机推荐
- LeetCode面试常见100题( TOP 100 Liked Questions)
LeetCode面试常见100题( TOP 100 Liked Questions) 置顶 2018年07月16日 11:25:22 lanyu_01 阅读数 9704更多 分类专栏: 面试编程题真题 ...
- Zend 3.3.0安装 ZendOptimizer 3.3.0 for Windows 稳定版 下载
用的某php网站系统今天打开时乱码了(zend 200407...),但phpmyadmin能正常使用: 搜索下,重新安装zend可以解决,系统上原来的版本是Zend 3.3.0:下了个,安装后果然把 ...
- 查找和杀掉占用GPU显存的进程
用只有2个G的显卡跑数据就需要在训练之前先把无关进程杀掉,防止跑到一半显存满了 nvidia-smi:显示当前GPU中的线程 kill -9 PID:输入PID以结束线程
- 实现双向绑定Proxy比defineproperty优劣如何?
前言 双向绑定其实已经是一个老掉牙的问题了,只要涉及到MVVM框架就不得不谈的知识点,但它毕竟是Vue的三要素之一. Vue三要素 响应式: 例如如何监听数据变化,其中的实现方法就是我们提到的双向绑定 ...
- css高度居中
1.已知元素高度 // 子盒子 #div1{ width:200px; height:200px; position: absolute; //父元素需要相对定位 top: 50%; left: 50 ...
- 通过javascript得到当前的日期和计算出该班级的平均分
某班的成绩出来了,现在老师要把班级的成绩打印出来. 格式要求: 1.显示打印的日期. 格式为类似“XXXX年XX月XX日 星期X” 的当前的时间. 2.计算出该班级的平均分(保留整数). 同学成绩数据 ...
- mysql一些语句
<!-- 报警量排行按创建时间每月来排行 --> <select id="alarmDaySort" resultType="alarm"&g ...
- JVM学习笔记(一):Java虚拟机和虚拟机内存区域
为什么Java程序需要运行在虚拟机上 因为Java在设计之初的跨平台特性,我们知道Java程序是运行在Java虚拟机上的.如果你要问为什么Java程序要运行在虚拟机上,我可以反问你几个问题. 为什么买 ...
- scroll js 原生
1.当前位置滚动: document.documentElement.scrollTop 当前位置: 有可能是0 window.scrollTo(,document.documentElement.s ...
- java 权限控制
网上或参考书中,对于java权限控制大多给出一张看似很整齐很好记实则不好理解的表格,我整理了一个2.0升级版,自认为会好理解很多,希望可以有所帮助. 同一包内 不同包内 修饰符 当前类 非当前类(含子 ...