题目链接

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 (网络流、最小割)的更多相关文章

  1. 【题解】 bzoj3894: 文理分科 (网络流/最小割)

    bzoj3894,懒得复制题面,戳我戳我 Solution: 首先这是一个网络流,应该还比较好想,主要就是考虑建图了. 我们来分析下题面,因为一个人要么选文科要么选理科,相当于两条流里面割掉一条(怎么 ...

  2. 【bzoj3774】最优选择 网络流最小割

    题目描述 小N手上有一个N*M的方格图,控制某一个点要付出Aij的代价,然后某个点如果被控制了,或者他周围的所有点(上下左右)都被控制了,那么他就算是被选择了的.一个点如果被选择了,那么可以得到Bij ...

  3. 【bzoj1143】[CTSC2008]祭祀river Floyd+网络流最小割

    题目描述 在遥远的东方,有一个神秘的民族,自称Y族.他们世代居住在水面上,奉龙王为神.每逢重大庆典, Y族都会在水面上举办盛大的祭祀活动.我们可以把Y族居住地水系看成一个由岔口和河道组成的网络.每条河 ...

  4. 【bzoj1797】[Ahoi2009]Mincut 最小割 网络流最小割+Tarjan

    题目描述 给定一张图,对于每一条边询问:(1)是否存在割断该边的s-t最小割 (2)是否所有s-t最小割都割断该边 输入 第一行有4个正整数,依次为N,M,s和t.第2行到第(M+1)行每行3个正 整 ...

  5. 【bzoj1976】[BeiJing2010组队]能量魔方 Cube 网络流最小割

    题目描述 一个n*n*n的立方体,每个位置为0或1.有些位置已经确定,还有一些需要待填入.问最后可以得到的 相邻且填入的数不同的点对 的数目最大. 输入 第一行包含一个数N,表示魔方的大小. 接下来 ...

  6. 【bzoj4177】Mike的农场 网络流最小割

    题目描述 Mike有一个农场,这个农场n个牲畜围栏,现在他想在每个牲畜围栏中养一只动物,每只动物可以是牛或羊,并且每个牲畜围栏中的饲养条件都不同,其中第i个牲畜围栏中的动物长大后,每只牛可以卖a[i] ...

  7. 【bzoj3438】小M的作物 网络流最小割

    原文地址:http://www.cnblogs.com/GXZlegend/p/6801522.html 题目描述 小M在MC里开辟了两块巨大的耕地A和B(你可以认为容量是无穷),现在,小P有n中作物 ...

  8. 【bzoj3144】[Hnoi2013]切糕 网络流最小割

    题目描述 输入 第一行是三个正整数P,Q,R,表示切糕的长P. 宽Q.高R.第二行有一个非负整数D,表示光滑性要求.接下来是R个P行Q列的矩阵,第z个 矩阵的第x行第y列是v(x,y,z) (1≤x≤ ...

  9. 【bzoj3894】文理分科 网络流最小割

    原文地址:http://www.cnblogs.com/GXZlegend 题目描述 文理分科是一件很纠结的事情!(虽然看到这个题目的人肯定都没有纠结过) 小P所在的班级要进行文理分科.他的班级可以用 ...

随机推荐

  1. 深入理解计算机系统 第十一章 网络编程 part1 第二遍

    客户端-服务器编程模型 每个网络应用都是基于客户端-服务器模型的.采用这个模型,一个应用是由一个服务器进程和一个或者多个客户端进程组成.服务器管理某种资源,并且通过操作这种资源来为它的客户端提供某种服 ...

  2. maraidb忘记数据密码

    一.概述 服务器上安装了maraidb 数据库,但是很久未使用过它,需要使用时,忘记了密码, 此时可以给它重新设置密码. 二.操作 修改密码 修改 /etc/my.cnf,修改下图红色区域位置,修改成 ...

  3. Sharepoint 开启App 配置App

    如果没有Enable app,打开app store的时候出出现错误: Sorry, apps are turned off. If you know who runs the server, tel ...

  4. 修改NPM默认全局安装路径

    场景: 最近在新电脑上鼓捣完环境后,打算切换下源,结果使用全局安装的nrm时提示找不到命令,之前都是这么用现在怎么不行了呢? 排查过程: 于是各种折腾,发现- g安装的插件目录在C盘中的某个路径中,后 ...

  5. 1 sql server中添加链接服务器

    1  链接到另一个sql server 的实例 exec sp_addlinkedserver @server= '服务器的地址',@srvproduct='SQL Server' go 分布式查询中 ...

  6. TensorFlow教程使用RNN生成唐诗

    本教程转载至:TensorFlow练习7: 基于RNN生成古诗词 使用的数据集是全唐诗,首先提供一下数据集的下载链接:https://pan.baidu.com/s/13pNWfffr5HSN79WN ...

  7. spring framework体系结构及模块jar依赖关系

    本文对于Spring的JAR包使用和配置,结合网友材料以spring 4.3.6.RELEASE版本为例,介绍spring框架结构和各模块对应JAR包以及模块间JAR依赖关系. 注:不同版本JAR包依 ...

  8. 基于Linux解决登录ssh客户端失败问题—sshd error: could not load host key

    当你ssh远程登录时,会发现ssh登录界面刚打开就会闪退,在查看主机日志消息中,就会看到如下错误: 然而问题的根源也就是这三个文件,无法正常加载ssh主机密钥. 而我们只需要将有问题的文件删除,然后重 ...

  9. redis 的启动、关闭 判断其是否在运行中

    #检查后台进程是否正在运行 ps -ef |grep redis ps aux | grep redis #检测6379端口是否在监听 netstat -lntp | grep 6379 #使用配置文 ...

  10. 执行sudo pip3 ...报错 Traceback (most recent call last): File "/usr/bin/pip3", line 9, in <module> from pip import main ImportError: cannot import name 'main'

    对于普通pip,把pip3改成pip即可,其他的修改一样 1.执行命令 sudo gedit /usr/bin/pip3 2.改成下面的形式 from pip import __main__ # 需要 ...