SG函数
来源:互联网 发布:托福 一个月 知乎 编辑:程序博客网 时间:2024/06/07 04:13
http://blog.csdn.net/luomingjun12315/article/details/45555495 (后面的sg函数及其实例取石子问题来源)
这篇博客讲的蛮好的。
组合游戏的和,假设有k个游戏,每次在当前回合,游戏者可以选择其中的一个游戏进行操作,然后保证其他游戏的局面不变,也就是说一次只能玩一个,不能同时玩。最后不能操作的人输。
其中nim游戏就是这种组合游戏特殊的一种,每一堆石头都可以看作一个游戏,每个游戏都相等。
Sprague-Grundy定理(SG定理)
游戏和的SG函数等于各个游戏SG函数的
nim和。所以我们可以将组合游戏中的每一个游戏单独考虑,从而简化了问题。
单堆nim游戏中的sg(x)=x;
SG函数:
首先定义mex(minimal excludant)运算,这是施加于一个集合的运算,表示最小的不属于这个集合的非负整数。例如mex{0,1,2,4}=3、mex{2,3,5}=0、mex{}=0。 对于任意状态 x , 定义 SG(x) = mex(S),其中 S 是 x 后继状态的SG函数值的集合。如 x 有三个后继状态分别为 SG(a),SG(b),SG(c),那么SG(x) = mex{SG(a),SG(b),SG(c)}。 这样 集合S 的终态必然是空集,所以SG函数的终态为 SG(x) = 0,当且仅当 x 为必败点P时。
【实例】取石子问题
有1堆n个的石子,每次只能取{ 1, 3, 4 }个石子,先取完石子者胜利,那么各个数的SG值为多少?
SG[0]=0,f[]={1,3,4},
x=1 时,可以取走1 - f{1}个石子,剩余{0}个,所以 SG[1] = mex{ SG[0] }= mex{0} = 1;
x=2 时,可以取走2 - f{1}个石子,剩余{1}个,所以 SG[2] = mex{ SG[1] }= mex{1} = 0;
x=3 时,可以取走3 - f{1,3}个石子,剩余{2,0}个,所以 SG[3] = mex{SG[2],SG[0]} = mex{0,0} =1;
x=4 时,可以取走4- f{1,3,4}个石子,剩余{3,1,0}个,所以 SG[4] = mex{SG[3],SG[1],SG[0]} = mex{1,1,0} = 2;
x=5 时,可以取走5 - f{1,3,4}个石子,剩余{4,2,1}个,所以SG[5] = mex{SG[4],SG[2],SG[1]} =mex{2,0,1} = 3;
以此类推…..
x 0 1 2 3 4 5 6 7 8….
SG[x] 0 1 0 1 2 3 2 0 1….
由上述实例我们就可以得到SG函数值求解步骤,那么计算1~n的SG函数值步骤如下:
1、使用 数组f 将 可改变当前状态 的方式记录下来。
2、然后我们使用 另一个数组 将当前状态x 的后继状态标记。
3、最后模拟mex运算,也就是我们在标记值中 搜索 未被标记值 的最小值,将其赋值给SG(x)。
4、我们不断的重复 2 - 3 的步骤,就完成了 计算1~n 的函数值。
代码实现如下:
//f[N]:可改变当前状态的方式,N为方式的种类,f[N]要在getSG之前先预处理 //SG[]:0~n的SG函数值 //S[]:为x后继状态的集合 int f[N],SG[MAXN],S[MAXN]; void getSG(int n){ int i,j; memset(SG,0,sizeof(SG)); //因为SG[0]始终等于0,所以i从1开始 for(i = 1; i <= n; i++){ //每一次都要将上一状态 的 后继集合 重置 memset(S,0,sizeof(S)); for(j = 0; f[j] <= i && j <= N; j++) S[SG[i-f[j]]] = 1; //将后继状态的SG函数值进行标记 for(j = 0;; j++) if(!S[j]){ //查询当前后继状态SG值中最小的非零值 SG[i] = j; break; } } }
HDU 1848
#include <iostream>#include <cstring>#include <algorithm>#include <cstdio>#include <cmath>#include <vector>using namespace std;const int maxn = 1005;#define LL long longint SG[maxn];int f[maxn];int S[maxn];void getfb(){ f[1]=1; f[2]=2; for(int i=3;;i++) { f[i]=f[i-1]+f[i-2]; if(f[i]>1000) break; }}void getSG(){ memset(SG,0,sizeof(SG)); for(int i=1;i<=1000;i++) { memset(S,0,sizeof(S)); for(int j=1;f[j]<=i;j++) { S[SG[i-f[j]]]=1; } for(int j=0;;j++) if(!S[j]) { SG[i]=j; break; } } //return SG[x];}int main(){ int m,n,p; getfb(); getSG(); while(scanf("%d %d %d",&m,&n,&p)!=EOF) { if(m==0&&n==0&&p==0) break; int x=SG[m]^SG[n]^SG[p]; if(x==0) printf("Nacci\n"); else printf("Fibo\n"); } return 0;}
- sg函数
- SG函数
- SG函数
- SG函数
- SG函数
- SG函数
- SG函数
- SG函数
- SG函数
- SG函数
- sg函数
- SG函数
- SG函数
- SG函数
- sg函数
- SG函数
- SG函数
- SG函数
- idr机制(32叉树)
- 每日vim系列0712
- (49)静态同步函数的锁是Class对象,例子:懒汉式--在多线程中存在问题--同步解决
- 日期和星期工具类
- ListView setAdapter has already been called 异常
- SG函数
- Android状态栏颜色修改,activity状态栏颜色,非布局写法
- SqlServerCLR聚合函数:分组连接字符串
- 轻松理解PV操作(转)
- 个简单的学生登陆系统的jsp实现之登陆身份验证
- unbuntu系统 在虚拟机下的安装
- 内部类向上转型为接口
- #423
- SAP物料主数据采购视图采购价值代码设置