巴什博奕
巴什博弈:一堆物品有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;
}