刚刚接触到sg函数突然感觉到原来可以这么好用,sg函数应该算是博弈论中比较经典的东西了。下面来说说sg函数:
从网上搜集资料终于能看懂了下面解释来自http://www.cnblogs.com/cj695/archive/2012/07/31/2617378.html,自己写不出来收藏了大神的思想。
他几乎可以解决博弈论中的所有问题。你可以将sg函数看作是一个深搜的的过程。而每一堆的石子就相当于图中间的节点。所以说整个sg函数的过程就是在对一个有向无环图进行dfs的过程。
sg函数的具体内容可以用一个公式来表示(虽然我最不喜欢公式,不过我还是得写。不然这没法说清楚):
sg(x) =mex{sg(y) : y ∈ F(x)}。其中{}内的是一个集合(只要上过高中都应该知道吧),在:右边的是该集合元素所满足的条件。sg(y)为该元素的值(其实就是一个递归的过程),F(x)貌似是该状态可以达到的状态。重点来了。。mex()函数表示是不在该集合中的最小的非负整数的值。比如mex({1,2,3})=0...mex({0,2,4})=1...mex({})=0..最后得出来的结果中。为0为必败点,不为零必胜点。。
接下来才是sg函数精妙之处了。假如说是在一个游戏中有多个石子堆该怎么办了。我们只需要把对每个石子堆进行sg函数的调用,将得到的所有的值进行异或。得出来的结果为0则该情形为必败态。否则为必胜态。。
此处贴上sg的一个模板:
int sg[N]; bool hash[N]; void sg_solve(int *s,int t,int N) //N求解范围代表一堆石子中的个数 S[]数组是可以每次取的值,t是s的长度。 { int i,j; memset(sg,0,sizeof(sg)); for(i=1;i<=N;i++) //枚举石子的个数 { memset(hash,0,sizeof(hash)); for(j=0;j<t;j++) if(i - s[j] >= 0) //枚举每次拿走的个数并标记 hash[sg[i-s[j]]] = 1; for(j=0;j<=N;j++) if(!hash[j]) break; sg[i] = j; //找到这个F[](该状态可以达到的状态)中不存在的最小的数 } }
下面写一个模板的应用:
题意:
是有1堆石子,两个人轮流从其中一堆拿,每次只能从一堆拿,可以拿任意个(大于0),也可以把一堆石子分成三堆(每堆大于0),问谁先赢。
void st() { int i,j,k; sg[0]=0; for(i=1;i<100;i++) { memset(vis,0,sizeof(vis)); for(j=0;j<i;j++) vis[sg[i-j]]=1; if(i>=3) { for(j=1;j<i-2;j++) for(k=1;k<i-2;k++) if(j+k<i) vis[sg[j]^sg[k]^sg[i-k-j]]=1; } j=0; while(vis[j]) j++; sg[i]=j; printf("%d ",sg[i]); } }