[0] 约定

\(r_i = \sum\limits_{j = 1}^{m}[A_{i,j}\le k]\)

\(r^{'}_i = \sum\limits_{j = 1}^{m}[B_{i,j}\le k]\)

\(c_j = \sum\limits_{i = 1}^{n}[A_{i,j}\le k]\)

\(c^{'}_j = \sum\limits_{i = 1}^{n}[B_{i,j}\le k]\)

[1] 分析思路

[1.1] 分析操作

对于 \(A\),若 \(\forall k,\{r_i|i\in[1,n]\} = \{r^{'}_i | i\in[1,n]\}\),则先行再列的排序后 \(A=B\)。

同理,若 \(\forall k,\{c_j|j\in[1,m]\} = \{c^{'}_j | j\in[1,m]\}\),则先列再行的排序后 \(A=B\)。

由此,确定 \(10\) 对排列 \(p_{1\sim n},q_{1\sim m}\),满足对于第 \(k\) 对排列,排序后 \(r_i = r^{'}_{p_i},c_i = c^{'}_{q_i}\)。

[1.2] 转化限制

对于每对排列单独考虑,需要满足 \([A_{i,j}\le k] = [B_{p_i,q_j}\le k]\)。

对于 \(k\) 相邻的两对排列,需要满足 \(B_{p_i,q_j}\le k \rightarrow B_{p\circ p^{'}_i,q\circ q^{'}_j}\le k + 1\)。

更进一步,\([B_{i,j}\le k] = [j\le r'_i] = [i\le c'_j]\)。

于是有 \(P:\forall i, j,j\le r'_i, p_i\le c'_{q_j}\)。

排序,使 \(c',r'\) 单调不增(以下写作 \(c,r\))。

\(P:\forall i, p_i\le c_{\max\limits_{j = 1}^{r_i}q_j}\)

发现只要符合 \(P\) 的排列就可以组合在一起。

[2] 如何计数

[2.1] DP 方程

发现需要计数的是满足以下条件的排列:

\(x_i = \max\limits_{j = 1}^{a_i}q_j\le m\)

\(p_i\le c_{x_{i}}\le n\)

考虑 DP 计数。

考虑与 树上异或 相似的分阶段计数。

令 \(f_{i,j}\) 代表 \(p_{1\sim i}\) 已确定,\(q_{a_{i} + 1\sim m}\) 已确定,且 \(x_i = j\) 的方案数。

发现可以在 \(f_{i - 1, j} \to f_{i,j}\) 的过程中计算 \(p_i\) 的贡献 \(j - (i - 1)\)。

可以在 \(i - 1\to i\) 的过程中计算 \(q_{a_{i} + 1\sim a_{i - 1}}\) 的贡献,用组合数算即可。

[2.2] 细节

就是无解或得到负数的情况。

[3] 后记

[3.0] 注意

以下全为作者自己口胡的,不要看,不要看!很可能是错的

[3.1] 正文

好的,让我们从头开始。

考虑排序对 \(A\) 造成的影响,发现比较难想。

发现 \(B\) 必然已经有序。

减少限制,若 \(A_{i,j}\in[0,1]\)。

[3.1.2] 约定

\(r_i\) 代表第 \(i\) 行 \(1\) 的个数。

\(c_j\) 代表第 \(j\) 列 \(1\) 的个数。

结论:

若可重集 \(r_A=r_B,c_A=c_B\) 那么排序后 \(A,B\) 相同。

证明:

考虑排序对 \(A\) 造成的影响。

先行再列,相当于对 \(r_i\) 排序。

先列再行,相当于对 \(c_j\) 排序。

回到原题,\(A_{i,j}\in[0,9]\)。

[3.1.3] 约定

\(r_{k,i}\) 代表第 \(i\) 行 \([A_{i,j}\le k] = 1\) 的个数。

\(c_{k,j}\) 代表第 \(i\) 行 \([A_{i,j}\le k] = 1\) 的个数。

再次考虑排序的意义,发现因为对于每个 \(k\),将 \(A\) 中元素视为 \([A_{i,j}\le k]\in[0,1]\),发现对于每个 \(k\) 排序得到的最终 \(A_{k}'\) 与 \(A_{i,j}\in[0,1]\) 的情况等价。

考虑通过 \(k\) 个 \(A_{k}'\) 还原 \(A'\),若 \(A_{k,i,j}'\) 为 \(1\) 说明 \(A'_{i,j}\le k\),同理 \(0\) 代表 \(A'_{i,j}>k\)。

直接计数应该也可以,但明显时间复杂度有问题。

发现 \(k\) 与 \(k + 1\) 的关系并没有那么紧密,因为只需要保证

\[\forall i\in[1,n],j\in[1,m],k\in[1,9],r_{k,i}\ge r_{k-1,i},c_{k,j}\ge c_{k-1,j}
\]

即可。

[4] 参考资料

E - RowCol/ColRow Sort Editorial by evima

AT_agc057_e 题解的更多相关文章

  1. 2016 华南师大ACM校赛 SCNUCPC 非官方题解

    我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...

  2. noip2016十连测题解

    以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...

  3. BZOJ-2561-最小生成树 题解(最小割)

    2561: 最小生成树(题解) Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1628  Solved: 786 传送门:http://www.lyd ...

  4. Codeforces Round #353 (Div. 2) ABCDE 题解 python

    Problems     # Name     A Infinite Sequence standard input/output 1 s, 256 MB    x3509 B Restoring P ...

  5. 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解

    题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...

  6. 2016ACM青岛区域赛题解

    A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Jav ...

  7. poj1399 hoj1037 Direct Visibility 题解 (宽搜)

    http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...

  8. 网络流n题 题解

    学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...

  9. CF100965C题解..

    求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...

  10. JSOI2016R3 瞎BB题解

    题意请看absi大爷的blog http://absi2011.is-programmer.com/posts/200920.html http://absi2011.is-programmer.co ...

随机推荐

  1. springboot 整合 pagehelper

    pom.xml <dependency> <groupId>com.github.pagehelper</groupId> <artifactId>pa ...

  2. Spring5.X的注解配置项目

    pom.xml <?xml version="1.0" encoding="UTF-8"?> <project xmlns="htt ...

  3. CentOS 8安装docker

    1.查看Linux内核(Docker最低支持CentOS 7 64位 内核3.10) uname -a 2.安装docker(输入yes,然后等待-) yum install docker 3.启动d ...

  4. [rCore学习笔记 014]批处理系统

    写在前面 本随笔是非常菜的菜鸡写的.如有问题请及时提出. 可以联系:1160712160@qq.com GitHhub:https://github.com/WindDevil (目前啥也没有 本章目 ...

  5. Midnight Commander (MC)

    Midnight Commander GNU Midnight Commander 是一个可视化文件管理器,根据 GNU 通用公共许可证获得许可,因此有资格成为自由软件.它是一个功能丰富的全屏文本模式 ...

  6. 靶机练习: y0usef

    靶机: y0usef 准备工作 靶机地址: https://download.vulnhub.com/y0usef/y0usef.ova MD5 校验:28c5d869b003be94b2d8ab4b ...

  7. Win10下安装LabelImg以及使用--LabelImg

    labelImg是图片标注软件,用于数据集的制作.标注等等.下面介绍labelImg的安装过程. 我用的是anaconda,所以以anaconda prompt作为终端: 在Anaconda Prom ...

  8. Docker Compose在Centos下的安装与卸载

    实践环境 Centos7.8 先决条件 已安装Docker Engine 安装Docker Compose 运行以下命令下载稳定版本Docker Compose $ sudo curl -L &quo ...

  9. MFC的消息反射

    这边我用对话框来测试mfc的消息反射 当对话框初始化完成之后,点击按钮,首先相应CDialog::OnCommand其中wpARAM(低八位nId,高八位子控件句柄) 和LPARAM(句柄) OnCo ...

  10. innodb存储引擎了解

    mysql常用的存储引擎分为innodb和myisam 其中innodb具有支持事务,执行行级锁,支持MVCC,外键,自动增长列,崩溃恢复等特性.并且mysql在5.5.5之后是数据的默认存储引擎 文 ...