http://aeghi.blog.hexun.com/33559875_d.html
看 POI Source 应该是某年 NOI 的题目
想了三节数学课终于想了出来,贴上相关的思考过程与大家分享,同时也算对这类问题做个总结。
题目大意:有两堆石子,各为 N、M 个,甲乙两人轮流取石子,每次可以从任意一堆里取走任意多个石子,也可以从两堆中取走相同数量的石子,但不能不取。取走最后一个石子的人为赢家。给出初始状态,问先手有无必胜策略。
数据范围:
RQNOJ 256:1 <= N、M <= 10000。
POJ 1067:1 <= N、M <= 1,000,000,000。
大致看完题目,想当然就知道这是一道博弈论的问题,最容易想的就是直接用博弈论的必败、必胜态进行动态规划求解。但是朴素的动态规划是 O(N * M) 的,如果做一些优化可能可以过掉 RQNOJ 的题目,但是对于 POJ 1067 来说就完全无能为力。所以我们尝试分析数据,看看有没有什么规律(以下用 (a, b) 表示两堆石子的个数,即游戏中的一个状态)。
列举了几个状态之后容易发现,必胜态的数目比必败态要多很多,所以我们先手工求出前几个必败态:
(1, 2)、(2, 1)、(3, 5)、(5, 3)、(4, 7)、(7, 4)、(6, 10)、(10, 6)……
首先回顾必胜态和必败态的朴素求法:
定理 0:一个状态是必败态,当且仅当它的所有后继状态都是必胜态;而一个状态是必胜态,只要它的后继状态有一个以上的必败态即可。
证明略去。
容易发现下面的定理:
定理 1:(a,b) 和 (b, a) 的胜负性是相同的(a <> b)。
证明:如果 (a, b) 是必胜态,那么将必胜策略中所有的操作,对第一堆的变为第二堆,对第二堆的变为第一堆,就构成 (b, a) 的必胜策略
定理 2:若 (a, b) 是必败态,则对于所有的 x <> a 和 y <> b,(x, b) 和 (a, y) 是必胜态。
证明:
对于 x > a 和 y > b,不管是哪一种情况,总可以从 x 堆或 y 堆中取出一定量的石子使当前状态变为必败态 (a, b),由定理 1,(x, b) 和 (a, y) 为必胜态。
对于 x < a 和 y < b,不管是哪一种情况,如果 (x, b) 或 (a, y) 是必败态的话,由上述可得 (a, b) 是必胜态,矛盾。故 (x, b) 和 (a, y) 均为为必胜态。
定理 3: 若 (a, b) 是必败态,则对于所有的 d > 0,(a + d, b + d) 是必胜态。
证明:
与定理 2 类似。
定理 4:在所有的必败态中,每个数字恰巧出现一次。
证明:
有了定理 1,对于对称的状态我们只需要处理其中一个,而两个数不会相同(相同的状态必然是必胜态),于是我们把每个状态中较小的数字放在前面,每行写一个状态,去掉括号并按照升序排列每行的第一个数,就构成了如下的矩阵:
1 2
3 5
4 7
6 10
……
观察这个矩阵,我们又可以得到新的定理:
定理 5:矩阵中每行第一个数恰巧是前面每一行中没有出现过的最小正整数。
证明:
由定理 4,矩阵中每个数字恰巧出现一次,而按照这个矩阵的定义,第二列的数总比同行第一列大,第一列又按照升序排列,所以每一行的第一个数正好为前面每一行中没有出现过的最小正整数。
定理 6:矩阵第 i 行的第二个数正好为第一个数加上 i
证明:
用数学归纳法。
1) 对于第一行显然成立
2) 若对于前 i - 1 行均成立,则所有的 (a[p], a[p] + p) (a[p] 为第 p 行第一个数,p < i) 均为必败态,那么考察第 i 行的状态 (a[i], a[i] + delta)。容易看出 delta >= i,因为如果 delta < i,一定可以通过一次操作变为前面出现过的必败态,那么这个状态就是必胜态。下面由 delta >= i,我们来说明 delta = i。
首先,我们考虑从第一堆中取出 p 个石子,得到状态 (a[i] - p, a[i] - p + delta),由定理 5,比 a[i] 小的数都在之前出现过,若 a[i] - p 出现在某一行的第一列,由于存在必败态 (a[i] - p, a[i] - p + d) (d < delta),故 (a[i] - p, a[i] - p + delta) 一定为必胜态(定理 2);若 a[i] - p 出现在某一行的第二列,由于第一列是单增的,因而其对应的第一列数必小于 a[i] + delta,故而也可推出其状态为必胜态。
对于从两堆石子中取出相同数目的情况与之类似,容易看出一定为必胜态。
于是,(a[i], a[i] + delta) 状态的胜负性只与状态 (a[i], a[i] + d) (d < delta) 有关。不难看出,delta = i 时恰为必败态,因为不论从第二堆中取出多少个石子,作为另一堆的第一堆石子并没有在之前出现过,所以得到的一定是一个必胜态,因而 (a[i], a[i] + delta) 为必败态,由定理 2 及定理 4 可得,原命题成立。即矩阵中第 i 行第二列的数等于同行第一列的数加上 i。
这时,我们所有的问题都转化到了矩阵上,只要能通过合适的方法表示出这个矩阵,我们就可以很好地解决原问题。
下面的过程可能需要比较高的数学技巧,首先给出我们需要的一个重要定理([x] 表示 x 的整数部分,{x} 表示 x 的小数部分,即 {x} = x - [x]):
定理 7(Betty 定理):如果存在正无理数 A, B 满足 1/A + 1/B = 1,那么集合 P = { [At], t ∈ Z+}、Q = { [Bt], t ∈ Z+} 恰为集合 Z+ 的一个划分,即:P ∪ Q = Z+,P ∩ Q = ø。
证明:暂时略去,将来补充。
考虑到 Betty 定理中“恰为 Z+ 的划分”这一说,这意味着,Z+ 中的每个数都恰好出现一次,这与上述矩阵的性质十分吻合。于是我们猜想每一行第一列的数满足 [Φi] 的形式。
于是我们得到每一行第二列的数为 [Φi] + i = [Φi + i] = [(Φ + 1)i]
我们的目的是要让 Z+ 中每个数都在这个矩阵中出现,于是考虑到 Betty 定理的条件,Φ 和 (Φ + 1) 应满足 1/Φ + 1/(Φ + 1) = 1。解这个方程,我们得到 Φ = (sqrt(5) + 1) / 2,于是 Φ + 1 = (sqrt(5) + 3) / 2。
Φ 恰为黄金分割比,这是多么令人惊奇的结论!
于是应用 Betty 定理,我们得到最终我们需要的定理:
定理 8:上述矩阵中每一行第一列的数为 [Φi],第二列的数为 [(Φ + 1)i],其中 Φ = (sqrt(5) + 1) / 2 为黄金分割比。
证明:由 Betty 定理显然得证。
有了定理 8,代码的实现就十分简单了,由于是数学算法,总复杂度为 O(1)。至此,本题完美解决。
总结:遇到这样困难的题目时,我们不应该轻言放弃。而应该仔细分析题目隐含的信息,学会分析和转化问题,从而找到问题的突破口,一举歼灭问题。
另附一道类似的题目,安徽 2004 年省选第二题,数字迷阵简述:
一个 M * N 的矩阵,第一行为 Fibonacci 数列,第二行开始,第 i 行的第一个数为前面各行中没有出现过的最小正整数,第二个数为第一个数 * 2 - (i - 1),第三个数开始满足 Fibonacci 关系。简列如下:
1 2 3 5 8 ……
4 7 11 18 29 ……
6 10 16 26 42 ……
8 13 21 34 55 ……
9 14 23 37 60 ……
……
……
问题是求出第 i 行第 j 列的数 mod P 的结果(即忽略高精度)。
数据规模记不清了,大概是 1 <= N, M <= 几十万吧。。