首先定义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}, //我们可以将 f[] , 从小到大排序
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；
注意：
f[]需要从小到大排序
1.可选步数为1~m的连续整数，直接取模即可，SG(x) = x % (m+1);
2.可选步数为任意步，SG(x) = x; ， Nim问题
3.可选步数为一系列不连续的数，用GetSG()计算
模板1如下：
// SG打表， 利用动态规划
// f[] 为堆中可取的石子个数
// x为堆中石子的总数， 如果有多个堆， 取堆中最多数为x.
int f[N],sg[N]; void getSG(int x) { sg[0]=0; for(int i=1; i<=x ; i++) // 这个x 是所有要用到的x 的最大值 { set<int>S; // 集合中无重复的元素， 且元素是按从小到大排序 for(int j=0 ; f[j]<= i; j++) // f[]数组 的下标是 从0开始的 S.insert(sg[i-f[j]]); int g=0; while(S.count(g) != 0) g++; sg[i] = g; } }
模板2如下：（dfs）
sg[1] =mex{sg[0]} = mex{0} sg[1] =1
sg[2] = mex{sg[1]} = mex{1} sg[2] = 0
sg[3] = mex{sg[2], sg[0]} = mex{0, 0} = 1
sg[4] = mex{sg[3], sg[1], sg[0]} = mex{1, 1, 0} sg[4] = 2
sg[5] = mex{sg[4], sg[2], sg[1]} = mex{2, 0, 1} sg[5] = 3
代码如下：
//注意 f数组要按从小到大排序
//SG函数要初始化为-1 对于每个集合只需初始化1遍
//n是集合f的大小 f[i]是定义的特殊取法规则的数组
int f[N_f],vis[N_f]; int sg[N_x],n; int SG_dfs(int x) { int i; if(sg[x]!=-1) return sg[x]; memset(vis,0,sizeof(vis)); for(i=0;i<n;i++) // n是集合f[]的大小 { if(x>=f[i]) { SG_dfs(x-f[i]); vis[sg[x-f[i]]]=1; } } int e; for(i=0;;i++) if(!vis[i]) { e=i; break; } return sg[x]=e; }
hdu 1848
三堆石子， 每堆只能去 fabonacci 数列。 直接对 三堆石子的总数 求 sg , 然后异或 ，即 sg[m]^sg[n]^sg[p] !=0. 。first 赢。
代码如下：
#include<iostream> #include<stdlib.h> #include<stdio.h> #include<math.h> #include<string.h> #include<string> #include<queue> #include<algorithm> #include<map> #include<set> #define N 1010 using namespace std; int f[17]; int sg[1005]; // F[16] = 1597,故Fibonacci 打表只需要达到16 void Fibo_init() { f[1] =1; f[2]= 2; for(int i=3;i<=16; i++) f[i] = f[i-1] +f[i -2]; } // SG函数打表 void init_SG(int x) { sg[0] =0 ; for(int i=1; i<=x ; i++) { set<int>S; for(int j=1; f[j]<=i ; j++) //f[]下标从1开始 S.insert(sg[i-f[j]]); int g=0; while(S.count(g)!= 0) g++; sg[i]= g; } } int main(){ Fibo_init(); init_SG(1000); int n,m,p,x; while(scanf("%d%d%d",&m,&n,&p)!=EOF) { if(m==0 && n==0 && p==0) break; x= sg[m]^sg[n]^sg[p]; if(x) puts("Fibo"); else puts("Nacci"); } return 0; }
到此，我们需要更加深刻的理解多堆，表示多个游戏。是否具有像积性函数 这样的性质？
假设有N堆， 每堆数为x1,x2,...,xn,从中按照一定的数量规则f[]取石子。
可以转换为， N 个棋子， 每个棋子为xi, 作为有向无环图的顶点， 分别计算 sg[xi].
我们有个结论： 即当多个棋子的有向图游戏的局面是first赢， 当且仅当 所有棋子所在位置的sg函数值 异或 为 0.
所以我们可以定义有向图游戏的和（sum of Graph Games）： 设G1，G2,、、、Gn 是n 个有向图游戏， 定义游戏G 是 G1,G2,..., Gn 的和， 游戏G 的移动规则为， 任选一个子游戏 Gi 并移动上面的棋子. Sprague- Grundy Theorem :就是： sg[G] = sg[G1] ^sg[G2] ^...^sg[Gn]. 也就是说， 游戏的和 的SG函数值 是它的 所有子游戏 的sg 函数值的 异或。
可以参考：http://www.docin.com/p-557540496.html
hdu 1847 , 一堆石子， 可以取的数组为 2的幂次（10次就可以了 ， 2^10=1024>1000）, 直接求sg[x]，如果sg[x] !=0 first 赢
Good Luck in CET-4 Everybody!
题目来源：http://acm.hdu.edu.cn/showproblem.php?pid=1847
代码如下：
#include<iostream> #include<stdlib.h> #include<stdio.h> #include<math.h> #include<string.h> #include<string> #include<queue> #include<algorithm> #include<map> #include<set> #define Inf 0x7fffffff // 0x 是数字0,而不是字母o #define N 1010 using namespace std; int f[11]; int sg[1005]; // SG函数打表 void init_SG(int x) { sg[0] =0 ; for(int i=1; i<=x ; i++) { set<int>S; for(int j=0; f[j]<=i ; j++) //f[]下标从1开始 S.insert(sg[i-f[j]]); int g=0; while(S.count(g)!= 0) g++; sg[i]= g; } } int main(){ int x; for(int i=0 ; i<=10; i++) f[i]=pow(2.0,(double)i); init_SG(1000); while(cin>>x) { if(sg[x]) puts("Kiki"); else puts("Cici"); } return 0; }
hdu 1849 Nim, m个棋子，m堆， 每堆可以取任意棋子， 直到为空。 ， 直接对输入的堆的个数（相当于位移长度）取 异或 即可。
题目来源：http://acm.hdu.edu.cn/showproblem.php?pid=1849
代码如下：
int main(){ int m,temp,x; while(cin>>m&&m) { x=0; for(int i=0;i<m ;i++) { cin>>temp; x^=temp; } if(x) puts("Rabbit Win!"); else puts("Grass Win!"); } return 0; }
hdu 1850 Nim博弈， 如果异或值x不为0， 表示first赢， first 可以走的方案数， 为 使异或值 x为0的temp[i] 数的个数， 即（temp[i] > (temp^x)）
题目来源：
http://acm.hdu.edu.cn/showproblem.php?pid=1850
代码如下：
#include<iostream> #include<stdlib.h> #include<stdio.h> #include<math.h> #include<string.h> #include<string> #include<queue> #include<algorithm> #include<map> #include<set> #define Inf 0x7fffffff // 0x 是数字0,而不是字母o #define N 105 using namespace std; int main(){ int m,temp[105],x,ans; while(cin>>m&&m) { x=0,ans=0; for(int i=0;i<m ;i++) { cin>>temp[i]; x^=temp[i]; } for(int i=0; i<m; i++) if(temp[i]> (temp[i]^x)) ans++; cout<<ans<<endl; } return 0; }
hdu 1851 题目来源：http://acm.hdu.edu.cn/showproblem.php?pid=1851
共有n堆， 每堆石子数为mi, 最多可取Li , 则 sg[mi] = mi %(Li+1), 然后对每堆 取异或x，x不为0， first赢。
代码如下：
#include<iostream> #include<stdlib.h> #include<stdio.h> #include<math.h> #include<string.h> #include<string> #include<queue> #include<algorithm> #include<map> #include<set> #define Inf 0x7fffffff // 0x 是数字0,而不是字母o #define N 15 using namespace std; int sg[N]; int main(){ int t,n,m,l; cin>>t; while(t--) { cin>>n; int x=0; for(int i=0;i<n;i++) { cin>>m>>l; x^=m%(l+1); } if(x) puts("No"); else puts("Yes"); } return 0; }