博弈论、SG定理和SG函数

巴什博奕

巴什博弈：一堆物品有n个，两个人轮流从这堆物品中取物，规定每次至少取一个，最多取m个。最后取光者得胜。

结论：

如果n%(m+1)==0 ,那么先手必败，否则后手必败

首先我们引入PN点概念
P-position：在该点时，先手必败
N-position：在该点时，先手必胜

还有关于PN点的性质：

• 终局（无法进行下一步操作的点）是P-position
• 只能进入P-position的点是P-position
• 可以进入P-position的点是N-position

证明：

如果n是0个，那么是P-position，结合上面的PN点的性质3，可以得到1~ m都是N-position，而m+1只能变成1~m，也即只能进入P-position，因此是N-position，以此类推，可以得到，当n是(m+1)的倍数时，先手必败，否则都可以一步进入P-position

SG定理与SG函数

上面的巴什博奕只有一堆石子，比较容易直接分析，那么如果是n堆石子，并且第i堆石子每次可以取1~mi个呢？

Sprague-Grundy定理：游戏和的SG函数等于各个游戏SG函数的Nim和（异或和）。

利用SG定理，我们可以将一个游戏分解成若干个子游戏来解决。

引入一个mex(minimal excludant)运算： 这是一个针对集合的运算；返回这个集合中不包含的最小的非负整数，例如mex({1,2,3})=0,mex({0,2,3})=1,mex({0,1,3})=2,

在这之前，你应该清楚，我们可以从终局往后推理，得到后面所有点的PN态，例如在上面的证明中，n=0是P-position，那么我们就可以推理得到1~ m都是N-position，这种情况下，我们就可以称n=(1~m)为n=0后继状态

sg函数：
sg(x)=mex(S)，其中S为x所有后继状态的 sg函数值 集合
也即 sg(x)=mex({sg(y1),sg(y2),sg(y3)…})

具体使用看下面的例题与代码：

练习题

代码：

#include <bits/stdc++.h>
using namespace std;
#define ll long long

ll f[1010]; //一次可以取的数量
ll sg[1010];
ll S[1010];
ll n = 1000, ct;
void getsg()
{
    fill(sg, sg + 1001, 0);
    for (int i = 1; i <= n; i++)
    {
        fill(S, S + 1001, 0); //清空集合
        for (int j = 1; f[j] <= 1000 && f[j] <= i; j++)
            S[sg[i - f[j]]] = 1; //记录所有后继状态的sg[]值
        for (int j = 0;; j++)    //模拟mex操作
            if (!S[j])
            {
                sg[i] = j;
                break;
            }
    }
}
void getfibo() //预处理菲波那切数列
{
    f[1] = 1;
    f[2] = 2;
    for (int i = 3; f[i - 1] <= 1000; i++)
        f[i] = f[i - 1] + f[i - 2];
}
int main()
{
    getfibo();
    getsg(); //sg打表
    ll a, b, c;
    while (cin >> a >> b >> c && (a || b || c))
    {
        if (sg[a] ^ sg[b] ^ sg[c]) //SG[ans]=sg[g1]+sg[g2]+sg[g3]
            cout << "Fibo" << endl;
        else
            cout << "Nacci" << endl;
    }
    return 0;
}

参考资料

1. 组合游戏 - SG函数和SG定理
2. 博弈论 SG函数