首先定义mex(minimal excludant)运算,这是施加于一个集合的运算,表示最小的不属于这个集合的非负整数。例如mex{0,1,2,4}=3、mex{2,3,5}=0、mex{}=0。
对于一个给定的有向无环图,定义关于图的每个顶点的Sprague-Grundy函数g如下:g(x)=mex{ g(y) | y是x的后继 },这里的g(x)即sg[x]
例如:取石子问题,有1堆n个的石子,每次只能取{1,3,4}个石子,先取完石子者胜利,那么各个数的SG值为多少?
sg[0]=0,f[]={1,3,4},
x=1时,可以取走1-f{1}个石子,剩余{0}个,mex{sg[0]}={0},故sg[1]=1;
x=2时,可以取走2-f{1}个石子,剩余{1}个,mex{sg[1]}={1},故sg[2]=0;
x=3时,可以取走3-f{1,3}个石子,剩余{2,0}个,mex{sg[2],sg[0]}={0,0},故sg[3]=1;
x=4时,可以取走4-f{1,3,4}个石子,剩余{3,1,0}个,mex{sg[3],sg[1],sg[0]}={1,1,0},故sg[4]=2;
x=5时,可以取走5-f{1,3,4}个石子,剩余{4,2,1}个,mex{sg[4],sg[2],sg[1]}={2,0,1},故sg[5]=3;
以此类推…
x 0 1 2 3 4 5 6 7 8…
sg[x] 0 1 0 1 2 3 2 0 1…
计算从1-n范围内的SG值。
f(存储可以走的步数,f[0]表示可以有多少种走法)
f[]需要从小到大排序,这个模版f是从1开始的。hash数组大小跟f[]大小差不多
1.可选步数为1~m的连续整数,直接取模即可,SG(x) = x % (m+1);
2.可选步数为任意步,SG(x) = x;
3.可选步数为一系列不连续的数,用GetSG()计算
-
//f[]:可以取走的石子个数
-
//sg[]:0~n的SG函数值
-
//hash[]:mex{}
-
int f[N],sg[N],hash[N];
-
void getSG(int n)
-
{
-
int i,j;
-
memset(sg,
0,
sizeof(sg));
-
for(i=
1;i<=n;i++)
-
{
-
memset(hash,
0,
sizeof(hash));
-
for(j=
1;f[j]<=i;j++)
-
hash[sg[i-f[j]]]=
1;
-
for(j=
0;j<=n;j++)
//求mes{}中未出现的最小的非负整数
-
{
-
if(hash[j]==
0)
-
{
-
sg[i]=j;
-
break;
-
}
-
}
-
}
-
}
转自:http://www.cnblogs.com/gongxijun/p/3297683.html
下边补充一点东西。
上边那个模版求hash时候并没有考虑f[j]有效长度,在某些题目中可以通过,比如这个。因为在求斐波那契数列时候肯定多求了一个,而就是因为这个会在求hash时候打破循环。
其实总的来说还是这个函数并不严密。因为有的时候f[]是有有效长度的,如果多出了这个长度就会出现错误,如果你的初值都是0,那么就会取到0,如果是-1,那么就会取到-1,肯定不对。比如这个题。
下面这两个模版应该就比较严密了,这个里边的f[]是从零开始的。转自:http://blog.csdn.net/primoblog/article/details/13376057
1、sg打表
-
//f[]:可以取走的石子个数
-
//sg[]:0~n的SG函数值
-
//hash[]:mex{}
-
int f[K],sg[N],hash[N];
-
void getSG(int n)
-
{
-
memset(sg,
0,
sizeof(sg));
-
for(
int i=
1; i<=n; i++) {
-
memset(hash,
0,
sizeof(hash));
-
for(
int j=
0; f[j]<=i && j < k; j++)
//k是f[]的有效长度
-
hash[sg[i-f[j]]]=
1;
-
for(
int j=
0; ; j++) {
//求mes{}中未出现的最小的非负整数
-
if(hash[j]==
0) {
-
sg[i]=j;
-
break;
-
}
-
}
-
}
-
}
2、dfs
-
//注意 S数组要按从小到大排序 SG函数要初始化为-1 对于每个集合只需初始化1遍
-
//n是集合s的大小 S[i]是定义的特殊取法规则的数组
-
int s[N],sg[N],n;
-
int getSG(int x)
-
{
-
if(sg[x]!=
-1)
-
return sg[x];
-
bool vis[M];
-
memset(vis,
0,
sizeof(vis));
-
for(
int i=
0; i<n; i++) {
-
if(x>=s[i])
-
vis[getSG(x-s[i])]=
1;
-
}
-
for(i=
0;; i++)
-
if(!vis[i]) {
-
sg[x]=i;
-
break;
-
}
-
return sg[x];
-
}