取石子游戏之威佐夫博弈

下面博文中有一部分是转载自ACM and JAVA的文章。

威佐夫博弈：有两堆各若干个物品，两个人轮流从某一堆或同时从两堆中取同样多的物品，规定每次至少取一个，多者不限，最后取光者得胜。
我们用(ak,bk)(ak<=bk,k=0,1,2,3...n)表示两堆物品的数量，并将这种形式记为局势。如果甲面对局势(0,0)，说明甲输了，我们把这种情况叫做奇异局势。

现在的问题我给你一个局势(a,b)，我们怎么判断它是不是奇异局势呢？

判断公式为：ak=[k(1+√5)/2]，bk=ak+k（k=0，1，2，…,n，[]:取整-高斯函数).可以令:k=b-a;那么如果:a=[(b-a)(1+√5)/2],那么这个局势为奇异局势。

如果面对非奇异局势，先拿者必胜；反之，则后拿者取胜。

例题：POJ 1067，直接用公式即可，下面为模版：

#include<iostream>
#include<cmath>
using namespace std;
#define S sqrt(5.0)
int main()
{   int a,b;
    while(cin>>a>>b)
    {   
        if(a>b) swap(a,b);
        if(a==int((b-a)*(S+1)/2)) cout<<0<<endl;
        else cout<<1<<endl; //非奇异局势,先拿着必胜 
    }
    return 0;
}

我们知道a0=b0=0,ak是没有在前面出现过的最小自然数，bk=ak+k，由上面的代码可以知道前几组的奇异局势：(0，0)、(1，2)、(3，5)、(4，7)、(6，10)、(8，1)、(9，15)、(11，18)、(12，20)。

奇异局势的三个性质为：

(1)、任何自然数都包含在一个且仅有一个奇异局势中。

原理：ak是没有在前面出现过的最小自然数,所以ak>ak-1，bk=ak+k>ak-1+(k-1)=bk-1>ak-1。

(2)、任意操作都可将奇异局势变为非奇异局势。

原理：若是只改变奇异局势(ak,bk)中的一个分量，由于另一个分量不可能存在于其他的奇异局势中，所以改变后一定是非奇异局势，如果使(ak，bk)的两个分量同时减少，则由于其差不变，且不可能是其他奇异局势的差，因此也是非奇异局势。

(3)、采用适当的方法，可以将非奇异局势变为奇异局势。

必败态的公式求解过程，需用到贝亚蒂定理：

若p,q∈R+(正实数集),且p,q为无理数，使得1/p+1/q=1，定义集合(贝亚蒂序列)P={ [np]:n∈Z+ },Q={ [nq]:n∈Z+ }，则P 和Q 构成正整数集的一个分划：P∩Q=空集，PUQ=Z+，

也就是说：如果两个正无理数的倒数之和是1，则任何正整数都可刚好以一种形式表示为不大于其中一个无理数的正整数倍的最大整数。

由性质1知：ak,bk刚好为贝亚蒂序列， 那么我可以令ak=[kp],bk=[kq]，且1/p+1/q=1。由于bk=ak+k=[kp]+k=[k(p+1)]=[kq],所以q=p+1,则有：1/p+1/(p+1)=1，可以求出p=(√5+1)/2。这就是必败态的通向公式：ak=[k(1+√5)/2]。

下面是一张必败态的表，很Nice的~