描述 Description
威佐夫博奕(Wythoff Game)，有2堆各不同数目个石子，两个人轮流从某一堆或同时从两堆中取同样多的物品，规定每次至少取1个，多者不限。取走最后石子的人获胜。POJ1067
变种玩法：“皇后登山”游戏，在空的围棋棋盘上放一个棋子，该棋子每次只能向上或向右或沿对角线向右上方向移动（相似国际象棋），可以移动任意格，但不能不移动，两人轮流移动棋子，先将棋子移动到右上角者赢，问先移棋者的必胜策略。《智力游戏中的数学方法》。
拓展问题：Moore’s Nimk。有n堆石子，每次最少选一堆，最多选m堆，每堆都可以拿任意正整数个石子。我们把每堆石子都写成二进制的形式，然后从低位到高位求每一位1的和。如果每一位1的和sum%(m+1)=0此时为必败态。显然，当m=1时就变成了我们最开始所讨论的最基础的Nim游戏了。
分析 Analysis
奇异局势的判定
前几个奇异局势是：（0，0）、（1，2）、（3，5）、（4，7）、（6，10）、（8，13）、（9，15）、（11，18）、（12，20）。。。
可以看出,a0=b0=0,ak是未在前面出现过的最小自然数,而 bk= ak + k，奇异局势有如下三条性质：
1、任何自然数都包含在一个且仅有一个奇异局势中。
由于ak是未在前面出现过的最小自然数，所以有ak > ak-1 ，而 bk= ak + k > ak-1 + k-1 = bk-1 > ak-1 。所以性质1。成立。
2、任意操作都可将奇异局势变为非奇异局势。
事实上，若只改变奇异局势（ak，bk）的某一个分量，那么另一个分量不可能在其他奇异局势中，所以必然是非奇异局势。如果使（ak，bk）的两个分量同时减少，则由于其差不变，且不可能是其他奇异局势的差，因此也是非奇异局势。
3、采用适当的方法，可以将非奇异局势变为奇异局势。假设面对的局势是（a,b），若 b = a，则同时从两堆中取走 a 个物体，就变为了奇异局势（0，0）；如果a = ak ，b > bk，那么，取走b - bk个物体，即变为奇异局势；如果 a = ak ， b < bk ,则同时从两堆中拿走 ak - ab - ak个物体,变为奇异局势（ ab - ak , ab - ak+ b - ak）；如果a > ak ，b= ak + k,则从第一堆中拿走多余的数量a - ak 即可；如果a < ak ，b= ak + k,分两种情况，第一种，a=aj （j < k）,从第二堆里面拿走 b - bj 即可；第二种，a=bj （j < k）,从第二堆里面拿走 b - aj 即可。
从如上性质可知，两个人如果都采用正确操作，那么面对非奇异局势，先拿者必胜；反之，则后拿者取胜。
那么任给一个局势（a，b），怎样判断它是不是奇异局势呢？
这个奇异局势的序列的通项公式可以表示为：（推导见附录）
Ak = [k*(1+sqrt(5.0)/2]
Bk = Ak + k
其中k=0，1，2，...,n ，方括号表示int取整函数。
有了这个通项式子，逆向的，对于某一个局势，只需要判断其A是否是黄金分割数的某个k的倍数，然后再确认B是否等于A+k即可。
奇妙的是其中出现了黄金分割数（1+√5）/2 = 1.618...,因此,由ak，bk组成的矩形近似为黄金矩形，由于2/（1+√5）=（√5-1）/2，可以先求出j=[a（√5-1）/2]，若a=[j（1+√5）/2]，那么a = aj，bj = aj + j，若不等于，那么a = aj+1，bj+1 = aj+1+ j + 1，若都不是，那么就不是奇异局势。然后再按照上述法则进行，一定会遇到奇异局势。
我们会发现这个序列的规律，设序列第k个奇异局势元素为(Ak,Bk)，k为自然数。那么，初始条件k=0时是，A0=B0=0，递推关系为下一个奇异局势的Ak是未在前面出现过的最小自然数，且Ak = Bk + k。
附录
Wythoff’s Game (威佐夫博弈)的解释
[类似解释亦可见于编程之美NIM（3）]Wythoff's Nim
有两堆石子，不妨先认为一堆有10，另一堆有15个，双方轮流取走一些石子，合法的取法有如下两种：
1)在一堆石子中取走任意多颗；
2)在两堆石子中取走相同多的任意颗；
约定取走最后一颗石子的人为赢家，求必败态(必胜策略)。
这个可以说是MR.Wythoff(于1907年提出此游戏)一生全部的贡献。不知道高斯取整函数与Beatty定理，所做的只能是找规律而已。可以先在http://www.cut-the-knot.org/pythagoras/withoff.shtml玩几局。
简单分析一下，容易知道两堆石头地位是一样的，我们用余下的石子数(a,b)来表示状态，并画在平面直角坐标系上。
用定理：有限个结点的无回路有向图有唯一的核中所述的方法寻找必败态。先标出(0,0)，然后划去所有(0,k),(k,0),(k,k)的格点；然后找y=x上方未被划去的格点，标出(1,2)，然后划去(1,k),(k,2),(1+k,2+k)，同时标出对称点(2,1)，划去(2,k),(1,k),(2+k,1+k)；然后在未被划去的点中在y=x上方再找出(3,5)。。。按照这样的方法做下去，如果只列出a<=b的必败态的话，前面的一些是(0,0),(1,2),(3,5),(4,7),(6,10),…
接下来就是找规律的过程了，忽略(0,0)，记第n组必败态为(a[n],b[n])
命题一：a[n+1]=前n组必败态中未出现过的最小正整数
[分析]：如果a[n+1]不是未出现的数中最小的，那么可以从a[n+1]的状态走到一个使a[n+1]更小的状态，和我们的寻找方法矛盾。
命题二：b[n]=a[n]+n
[分析]：归纳法：若前k个必败态分别为 ，下证：第k+1个必败态为
从该第k+1个必败态出发，一共可能走向三类状态，从左边堆拿走一些，从右边堆拿走一些，或者从两堆中拿走一些．下面证明这三类都是胜态．
情况一：由命题一，任意一个比a[k+1]小的数都在之前的必败态中出现过，一旦把左边堆拿少了，我们只要再拿成那个数相应的必败态即可。
情况二（从右边堆拿走不太多）：这使得两堆之间的差变小了，比如拿成了 ，则可再拿成
；
情况二（从右边堆拿走很多）：使得右边一堆比左边一堆更少，这时类似于情况一，比如拿成了 (其中a[m] ；
情况三：比如拿成 ，则可再拿成
．
综上所述，任何从 出发走向的状态都可以走回核中．故原命题成立．
以上两个命题对于确定(a[n],b[n])是完备的了，给定(0,0)然后按照这两个命题，就可以写出(1,2),(3,5),(4,7),…
这样我们得到了这个数列的递推式，以下我们把这两个命题当成是(a[n],b[n])的定义。
先证明两个性质：
性质一：核中的a[n],b[n]遍历所有正整数。
[分析]：由命题一，二可得a[n],b[n]是递增的，且由a[n]的定义显然。
性质二：A={a[n]:n=1,2,3,…},B={b[n]:n=1,2,3,…}，则集合A,B不交。
[分析]：由核是内固集，显然。
实际上a[n]和b[n]就是一个Beatty序列。
Betty 定理
如果存在正无理数 A, B 满足 1/A + 1/B = 1，那么集合 P = { [A*t], t ∈ Z+}、Q = { [B*t], t ∈ Z+} 恰为集合 Z+ 的一个划分，即：P ∪ Q = Z+，P ∩ Q = ø。证明见附录2。
考虑到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。 ，到此，我们找到了该必败态的通项公式。
实际上这组Beatty序列还有一些别的性质，比如当一个数是Fibonacci数的时候，另一个数也是Fibonacci数；而且两者的比值也越来越接近黄金比，这些性质在得到通项公式之后不难证明。
启示：首先用定理所说的方法找核，然后给出核的规律（递推，或是通项）并且证明。
附上一张对应的必败态图
2.Betti theorem：
设a、b是正无理数且1/a+1/b=1。记P={[na]|n为任意的正整数}，Q={[nb]|n为任意的正整数}，([x]'指的是取x的整数部分)则P与Q是Z+的一个划分，即P∩Q为空集且P∪Q为正整数集合Z+。
证明：
因为a、b为正且1/a+1/b=1，则a、b>1，所以对于不同的整数n，[na]各不相同，类似对b有相同的结果。因此任一个整数至多在集合P或Q中出现一次。
* 现证明P∩Q为空集(反证法)：假设k为P∩Q的一个整数，则存在正整数m、n使得[ma]=[nb]=k。即k < ma、nb<k+1，等价地改写不等式为
m/(k+1)< 1/a < m/k及n/(k+1)< 1/b < n/k。相加起来得 (m+n)/(k+1) < 1 < (m+n)/k，即 k < m+n < k+1。这与m、n为整数有矛盾，所以P∩Q为空集。
*现证明Z+=P∪Q：已知P∪Q是Z+的子集，剩下来只要证明Z+是P∪Q的子集。(反证法)假设Z+\(P∪Q)有一个元素k，则存在正整数m、n使得[ma]< k <[(m+1)a]、[nb]< k <[(n+1)b]。 由此得ma < k ≦[ (m+1)a]-1<(m+1)a -1（因为a是无理数），类似地有nb < k ≦[ (n+1)b]-1<(n+1)b -1。等价地改写为 m/k < 1/a < (m+1)/(k+1)及n/k < 1/b < (n+1)/(k+1)。两式加起来，得
(m+n)/k < 1 < (m+n+2)/(k+1)，即m+n < k < k+1 < m+n+2。这与m, n, k皆为正整数矛盾。
所以Z+=P∪Q。[贝蒂定理]
代码 Code
“皇后登山”游戏
/****************************************************************************/
/* 编程之美 - NIM捡石子问题 皮皮 2014-9-10 */
/****************************************************************************/
#include <assert.h>
#include <stdio.h>
#include <math.h>
/* 黄金分割数最优算法 */
static int nim(int x, int y){
double a = (sqrt(5.0) + 1) / 2; //goldenNum 1.618...
int n = abs(x - y); //bn - an
x = x < y? x : y; //取小者
return ( x != int(a * n) );
}
int main(){
int x, y;
while(scanf("%d%d", &x, &y) == 2)
printf("%d\n", nim(x, y));
return 0;
}