#L3771. 「APIO2022」火星

「APIO2022」火星

题目描述

在 LibreOJ,本题仅保证使用 GNU C++1717 编译将得到预期的编译结果,使用其他编译器版本不保证编译结果一定符合预期。

你们晓得,法老们是最先去过外太空的人。他们发射过首次登陆行星图特摩斯一世(Thutmus I,现在一般叫它 火星)的飞船。行星的表面可以建模成由方形单元构成的 (2n+1)×(2n+1)(2n + 1) \times (2n + 1) 网格,其中每个单元中或者为陆地、或者为水域。对于第 ii 行第 jj 列(0i,j2n0 \leq i, j \leq 2n)的单元,如果单元中为陆地,则其状态表示为 s[i][j]=’1’s[i][j] = \text{'1'};如果单元中为水域,则表示为 s[i][j]=’0’s[i][j] = \text{'0'}

如果在两个陆地单元之间存在某条仅由陆地单元构成的路径,而且路径中每两个连续的前后单元都有公共边,则称这两个陆地单元是连通的。行星上的岛屿被定义为两两连通的陆地单元的极大集合。

飞船的任务是统计该行星上岛屿的数量。然而,考虑到飞船的上古电脑,这事儿并不容易。电脑的内存储器 hh 以一个 (2n+1)×(2n+1)(2n + 1) \times (2n + 1) 的二维数组的形式存储数据,且数组的每个位置上可以保存长度为 100100 的字符串,串中的每个字母为 ’0’\text{'0'}(ASCII 码 4848)或 ’1’\text{'1'}(ASCII 码 4949)。初始时,存储器的每个位置的第 00 位记录的是上述网格中每个单元的状态,即 h[i][j][0]=s[i][j]h[i][j][0] = s[i][j](对所有 0i,jn0 \leq i, j \leq n)。hh 中的其他位在初始时都被置为 ’0’\text{'0'}(ASCII 码 4848)。

在处理存储器中的数据时,电脑只能访问存储器中的 3×33 \times 3 区块,并且改写该区块左上角位置的值。说得更正式一点,电脑可以访问 h[ii+2][jj+2]h[i \dots i + 2][j \dots j + 2]0i,j2(n1)0 \leq i, j \leq 2(n − 1))中的值,并且改写 h[i][j]h[i][j] 中的值。在下文中,该过程被叫做处理单元 (i,j)(i, j)

为了解决电脑能力的局限,法老们搞出了下面的套路:

  • 电脑可以分成 nn 个阶段来操作存储器。
  • 在阶段 kk0kn10 \leq k \leq n − 1),令 m=2(nk1)m = 2 (n − k − 1), 电脑将对所有的 0i,jm0 \leq i, j \leq m,按照 ii 的升序以及每个 iijj 的升序,处理单元 (i,j)(i, j)。换句话说,电脑将按照如下顺序处理这些单元:$(0, 0), (0, 1), \dots, (0, m), (1, 0), (1, 1), \dots, (1, m), \dots, (m, 0), (m, 1), \dots, (m, m)$。
  • 在最后一个阶段(k=n1k = n − 1),电脑仅处理单元 (0,0)(0, 0)。该阶段结束后,写入到 h[0][0]h[0][0] 的值应该等于行星上的岛屿数量,而且该值应以字符串的形式表示成二进制,其中最低有效位对应于字符串的首字符

下图给出了电脑操作某个 5×55 \times 5n=2n = 2)存储器的方式。蓝色单元表示该单元正在被改写,而着色的单元则表示被处理的子数组。

在阶段 00,电脑将以如下顺序处理下面的子数组:

(图略)

在阶段 11,电脑将仅处理一个子数组:

(图略)

你的任务是给出一个方法,让电脑能在给定的操作方式下,统计出行星图特摩斯一世上的岛屿数量。


实现细节

请在程序开头引入库 mars.h。你需要实现下面的函数:

string process(string[][] a, int i, int j, int k, int n);
  • a:一个 3×33 \times 3 数组,表示正在被处理的子数组。特别说明,有 a=h[ii+2][jj+2]a = h[i \dots i + 2][j \dots j + 2],这里 aa 中的每个元素均为长度恰好为 100100 的字符串,而且串中的字符为 ’0’\text{'0'}(ASCII 码 4848)或 ’1’\text{'1'}(ASCII 码 4949)。
  • i, j:电脑当前正在处理的单元的行号和列号。
  • k:当前阶段的序号。
  • n:阶段总数,同时也是行星表面的大小,此时行星表面包含 (2n+1)×(2n+1)(2n + 1) \times (2n + 1) 个单元。

该函数应返回一个长度为 100100 的二进制表示字符串。返回值将保存在电脑存储器中的 h[i][j]h[i][j] 处。

k=n1k = n − 1 时,是该函数的最后一次调用。在此次调用中,函数应以字符串的形式返回行星上的岛屿数量的二进制表示,其最低有效位对应下标 00 处的字符(二进制字符串的首字符),次低有效位对应下标 11 处的字符,以此类推。

该函数必须独立于任何的静态或全局变量,且其返回值应仅依赖于传递给该函数的参数。

每个测试用例包括 TT 个独立的场景(也就是说,不同的行星表面情形)。你的函数在每个场景上的行为,必须与这些场景的顺序无关,因为对同一场景的 process 函数调用可能不是连续发生的。但是,可以确保对每个场景,会按照题面所描述的顺序来调用函数 process

此外,对每个测试用例,你的程序可能会同时运行多个实例。内存限制和 CPU 用时限制将施加在所有这些实例的总和上。任何故意在这些实例之间偷偷传递数据的行为,都将被认定为作弊,选手可能会因此被取消比赛资格。

特别说明,在调用函数 process 时保存在静态或全局变量中的信息,不保证在下次调用时可以读出。


约束条件

  • 1T101 \leq T \leq 10
  • 1n201 \leq n \leq 20
  • s[i][j]s[i][j]’0’\text{'0'}(ASCII 码 4848)或 ’1’\text{'1'}(ASCII 码 4949)(对所有 0i,j2n0 \leq i, j \leq 2n);
  • h[i][j]h[i][j] 的长度恰好为 100100(对所有 0i,j2n0 \leq i, j \leq 2n);
  • h[i][j]h[i][j] 中的每个字符均为 ’0’\text{'0'}(ASCII 码 4848)或 ’1’\text{'1'}(ASCII 码 4949)(对所有 0i,j2n0 \leq i, j \leq 2n)。

对函数 process 的每次调用,都有:

  • 0kn10 \leq k \leq n − 1
  • 0i,j2(nk1)0 \leq i, j \leq 2 (n − k − 1)

子任务

  • 66 分)n2n \leq 2
  • 88 分) n4n \leq 4
  • 77 分) n6n \leq 6
  • 88 分) n8n \leq 8
  • 77 分) n10n \leq 10
  • 88 分) n12n \leq 12
  • 1010 分) n14n \leq 14
  • 2424 分) n16n \leq 16
  • 1111 分) n18n \leq 18
  • 1111 分) n20n \leq 20

样例 1

输入

1
1
1 0 0
1 1 0
0 0 1

输出

2

考虑 n=1n = 1 的样例,其中 ss 如下所示:

'1' '0' '0'
'1' '1' '0'
'0' '0' '1'

在本例中,行星表面包括 3×33 \times 3 个单元,其中有 22 个岛屿。对函数 process 的调用至多只有 11 个阶段。

在阶段 00,评测程序将调用函数 process 恰好一次:

process([["100","000","000"],["100","100","000"],["000","000","100"]],0,0,0,1);

该函数应返回 "0100..."(省略的位全部为零),这里二进制的 ...0010 等于十进制的 22。注意,这里省略了 9696 个零并用 \dots 来代替。


样例 2

输入

1
2
1 1 0 1 1
1 1 0 0 0
1 0 1 1 1
0 1 0 0 0
0 1 1 1 1

输出

4

考虑 n=2n = 2 的样例,其中 ss 如下所示:

'1' '1' '0' '1' '1'
'1' '1' '0' '0' '0'
'1' '0' '1' '1' '1'
'0' '1' '0' '0' '0'
'0' '1' '1' '1' '1'

在本例中,行星表面包括 5×55 \times 5 个单元,其中有 44 个岛屿。对函数 process 的调用至多只有 22 个阶段。

在阶段 00,评测程序将调用函数 process 99 次:

process([["100","100","000"],["100","100","000"],["100","000","100"]],0,0,0,2);
process([["100","000","100"],["100","000","000"],["000","100","100"]],0,1,0,2);
process([["000","100","100"],["000","000","000"],["100","100","100"]],0,2,0,2);
process([["100","100","000"],["100","000","100"],["000","100","000"]],1,0,0,2);
process([["100","000","000"],["000","100","100"],["100","000","000"]],1,1,0,2);
process([["000","000","000"],["100","100","100"],["000","000","000"]],1,2,0,2);
process([["100","000","100"],["000","100","000"],["000","100","100"]],2,0,0,2);
process([["000","100","100"],["100","000","000"],["100","100","100"]],2,1,0,2);
process([["100","100","100"],["000","000","000"],["100","100","100"]],2,2,0,2);

假定上面调用得到的返回值分别为 "011""000""000""111""111""011""110""010""111",被省略的位均为零。因此,在阶段 00 结束后,hh 将保存有如下的值:

"011", "000", "000", "100", "100"
"111", "111", "011", "000", "000"
"110", "010", "111", "100", "100"
"000", "100", "000", "000", "000"
"000", "100", "100", "100", "100"

在阶段 11,评测程序将调用函数 process 一次:

process([["011","000","000"],["111","111","011"],["110","010","111"]],0,0,1,2);

最后,本次函数调用应返回 "0010000..."(被省略的位均为零),这里二进制的 ...0000100 等于十进制的 44。注意这里省略了 9393 个零并用 \dots 来代替。


评测程序示例

评测程序示例读取如下格式的输入:

  • 11 行:TT
  • ii 个区块(0iT10 \leq i \leq T − 1):该区块表示第 ii 个场景。
    • 11 行:nn
    • 2+j2 + j 行(0j2n0 \leq j \leq 2n):s[j][0],s[j][1],,s[j][2n]s[j][0], s[j][1], \dots, s[j][2n]

评测程序示例将按照如下格式打印出结果:

  • 1+i1 + i 行(0iT10 \leq i \leq T − 1):在第 ii 个场景上,函数 process 最后一次的返回值的十进制表示。