SG 大法(Sprague-Grundy函数)

SG函数的定义:

g(x) = mex ( sg(y) |y是x的后继结点 )

其中mex(x)(x是一个自然是集合)函数是x关于自然数集合的补集中的最小值,比如x={0,1,2,4,6} 则mex(x)=3;

什么是后继结点?

所谓后继结点就是当前结点经过一个操作可以变成的状态。比如对于取4石子游戏,假如每次可以取的数目是1,2,4,当前的石子数目也就是当前状态是5,那么5的后继结点就是{5-1, 5-2, 5-4}={4,3,1};

如果5的三个后继结点的SG函数值分别为0,1,3,那么5的SG值就是集合{0,1,3}的补集的最小元素,也就是2。

关于整个游戏的sg值之和sum,定义sum=sg1 ^ sg2 ^ sg3 ^ ……sgn.  其中^表示按位异或运算。

结论:一个游戏的初始局面是必败态当且仅当sum=0。


插一段大神的论文



上一期的文章里我们仔细研究了Nim游戏,并且了解了找出必胜策略的方法。但如果把Nim的规则略加改变,你还能很快找出必胜策略吗?比如说:有n堆石子,每次可以从第1堆石子里取1颗、2颗或3颗,可以从第2堆石子里取奇数颗,可以从第3堆及以后石子里取任意颗……这时看上去问题复杂了很多,但相信你如果掌握了本节的内容,类似的千变万化的问题都是不成问题的。
现在我们来研究一个看上去似乎更为一般的游戏:给定一个有向无环图和一个起始顶点上的一枚棋子,两名选手交替的将这枚棋子沿有向边进行移动,无法移动者判负。事实上,这个游戏可以认为是所有Impartial Combinatorial Games的抽象模型。也就是说,任何一个ICG都可以通过把每个局面看成一个顶点,对每个局面和它的子局面连一条有向边来抽象成这个“有向图游戏”。下面我们就在有向无环图的顶点上定义Sprague-Garundy函数。

首先定义mex(minimal excludant)运算,这是施加于一个集合的运算,表示最小的不属于这个集合的非负整数。例如mex{0,1,2,4}=3、mex{2,3,5}=0、mex{}=0。

对于一个给定的有向无环图,定义关于图的每个顶点的Sprague-Garundy函数g如下:g(x)=mex{ g(y) | y是x的后继 }。

来看一下SG函数的性质。首先,所有的terminal position所对应的顶点,也就是没有出边的顶点,其SG值为0,因为它的后继集合是空集。然后对于一个g(x)=0的顶点x,它的所有后继y都满足g(y)!=0。对于一个g(x)!=0的顶点,必定存在一个后继y满足g(y)=0。

以上这三句话表明,顶点x所代表的postion是P-position当且仅当g(x)=0(跟P-positioin/N-position的定义的那三句话是完全对应的)。我们通过计算有向无环图的每个顶点的SG值,就可以对每种局面找到必胜策略了。但SG函数的用途远没有这样简单。如果将有向图游戏变复杂一点,比如说,有向图上并不是只有一枚棋子,而是有n枚棋子,每次可以任选一颗进行移动,这时,怎样找到必胜策略呢?

让我们再来考虑一下顶点的SG值的意义。当g(x)=k时,表明对于任意一个0<=i<k,都存在x的一个后继y满足g(y)=i。也就是说,当某枚棋子的SG值是k时,我们可以把它变成0、变成1、……、变成k-1,但绝对不能保持k不变。不知道你能不能根据这个联想到Nim游戏,Nim游戏的规则就是:每次选择一堆数量为k的石子,可以把它变成0、变成1、……、变成k-1,但绝对不能保持k不变。这表明,如果将n枚棋子所在的顶点的SG值看作n堆相应数量的石子,那么这个Nim游戏的每个必胜策略都对应于原来这n枚棋子的必胜策略!

对于n个棋子,设它们对应的顶点的SG值分别为(a1,a2,...,an),再设局面(a1,a2,...,an)时的Nim游戏的一种必胜策略是把ai变成k,那么原游戏的一种必胜策略就是把第i枚棋子移动到一个SG值为k的顶点。这听上去有点过于神奇——怎么绕了一圈又回到Nim游戏上了。

其实我们还是只要证明这种多棋子的有向图游戏的局面是P-position当且仅当所有棋子所在的位置的SG函数的异或为0。这个证明与上节的Bouton‘s Theorem几乎是完全相同的,只需要适当的改几个名词就行了。

刚才,我为了使问题看上去更容易一些,认为n枚棋子是在一个有向图上移动。但如果不是在一个有向图上,而是每个棋子在一个有向图上,每次可以任选一个棋子(也就是任选一个有向图)进行移动,这样也不会给结论带来任何变化。

所以我们可以定义有向图游戏的和(Sum of Graph Games):设G1、G2、……、Gn是n个有向图游戏,定义游戏G是G1、G2、……、Gn的和(Sum),游戏G的移动规则是:任选一个子游戏Gi并移动上面的棋子。Sprague-Grundy Theorem就是:g(G)=g(G1)^g(G2)^...^g(Gn)。也就是说,游戏的和的SG函数值是它的所有子游戏的SG函数值的异或。

再考虑在本文一开头的一句话:任何一个ICG都可以抽象成一个有向图游戏。所以“SG函数”和“游戏的和”的概念就不是局限于有向图游戏。我们给每个ICG的每个position定义SG值,也可以定义n个ICG的和。所以说当我们面对由n个游戏组合成的一个游戏时,只需对于每个游戏找出求它的每个局面的SG值的方法,就可以把这些SG值全部看成Nim的石子堆,然后依照找Nim的必胜策略的方法来找这个游戏的必胜策略了!

回到本文开头的问题。有n堆石子,每次可以从第1堆石子里取1颗、2颗或3颗,可以从第2堆石子里取奇数颗,可以从第3堆及以后石子里取任意颗……我们可以把它看作3个子游戏,第1个子游戏只有一堆石子,每次可以取1、2、3颗,很容易看出x颗石子的局面的SG值是x%4。第2个子游戏也是只有一堆石子,每次可以取奇数颗,经过简单的画图可以知道这个游戏有x颗石子时的SG值是x%2。第3个游戏有n-2堆石子,就是一个Nim游戏。对于原游戏的每个局面,把三个子游戏的SG值异或一下就得到了整个游戏的SG值,然后就可以根据这个SG值判断是否有必胜策略以及做出决策了。其实看作3个子游戏还是保守了些,干脆看作n个子游戏,其中第1、2个子游戏如上所述,第3个及以后的子游戏都是“1堆石子,每次取几颗都可以”,称为“任取石子游戏”,这个超简单的游戏有x颗石子的SG值显然就是x。其实,n堆石子的Nim游戏本身不就是n个“任取石子游戏”的和吗?

所以,对于我们来说,SG函数与“游戏的和”的概念不是让我们去组合、制造稀奇古怪的游戏,而是把遇到的看上去有些复杂的游戏试图分成若干个子游戏,对于每个比原游戏简化很多的子游戏找出它的SG函数,然后全部异或起来就得到了原游戏的SG函数,就可以解决原游戏了。这种“分而治之”的思想在下一节介绍的“翻硬币游戏”中将被应用得淋漓尽致。还是敬请期待。



SG打表模板1

 1 int f[K],sg[maxn],mex[K];  //K为集合的大小   maxn一个堆的大小
 2 int k ;  //k为集合的大小
 3 void getsg(int H)
 4 {
 5     int i,j;
 6     sg[0]=0;
 7     for (i=1;i<H;i++){
 8         memset(mex,0,sizeof(mex));
 9         j=1;
10         while (j<=k && i>=f[j]){
11             mex[sg[i-f[j]]]=1;
12             j++;
13         }
14         j=0;
15         while (mex[j]) j++;
16         sg[i]=j;
17     }
18 }  

SG打表模板2

//f[]  可以取走的石子数  从f[1]开始  f[]里的值要从小到大排序
//sg[]   0~n的sg函数值
//hash[]  mex{}
int f[maxn],sg[maxn],hash[maxn];
void getsg(int n)
{
    memset(sg,0,sizeof(sg));
    for(int i=1;i<=n;i++)
    {
        memset(hash,0,sizeof(hash));
        for(int j=1;f[j]<=i;j++)
            hash[sg[i-f[j]]]=1;
        for(int j=0;j<=n;j++)
        {
            if(hash[j]==0)
            {
                sg[i]=j;
                break;
            }
        }
    }
}  

HDU 1536  和 1944

题意:
首先输入K 表示一个集合的大小 之后输入集合 表示对于这对石子只能去这个集合中的元素的个数
之后输入 一个m 表示接下来对于这个集合要进行m次询问
之后m行 每行输入一个n 表示有n个堆 每堆有n1个石子 问这一行所表示的状态是赢还是输 如果赢输入W否则L

Sample Input
2 2 5 // 集合大小 集合里的数
3 //询问次数
2 5 12 // 多少堆 堆的大小
3 2 4 7
4 2 3 7 12
5 1 2 3 4 5
3
2 5 12
3 2 4 7
4 2 3 7 12
0

Sample Output
LWW //L为先手败
WWL

 1 # include <iostream>
 2 # include <cstdio>
 3 # include <cmath>
 4 # include <cstring>
 5 # include <algorithm>
 6 using namespace std ;
 7
 8 int f[110],sg[10010],mex[110];
 9 int k ;  //k为集合的大小
10 void getsg(int H)
11 {
12     int i,j;
13     sg[0]=0;
14     for (i=1;i<H;i++){
15         memset(mex,0,sizeof(mex));
16         j=1;
17         while (j<=k && i>=f[j]){
18             mex[sg[i-f[j]]]=1;
19             j++;
20         }
21         j=0;
22         while (mex[j]) j++;
23         sg[i]=j;
24     }
25 }
26
27 int main ()
28 {
29 //    freopen("in.txt","r",stdin) ;
30
31     while (scanf("%d" , &k) , k)
32     {
33         int i ;
34         for (i = 1 ; i <= k ; i++)
35         {
36             scanf("%d" , &f[i]) ;
37         }
38         sort(f+1 , f+1+k) ;
39         getsg(10010) ;
40         int m ;
41         scanf("%d" , &m) ;
42         while (m--)
43         {
44             int n ;
45             scanf("%d" , &n) ;
46             int x ;
47             int sum = 0 ;
48             for (i = 1 ; i <= n ; i++)
49             {
50                 scanf("%d" , &x) ;
51                 sum ^= sg[x] ;
52             }
53             if (sum == 0)
54               printf("L") ;
55             else
56               printf("W") ;
57         }
58         printf("\n") ;
59     }
60
61     return 0 ;
62 }

HDU 1848

1、 这是一个二人游戏;
2、 一共有3堆石子,数量分别是m, n, p个;
3、 两人轮流走;
4、 每走一步可以选择任意一堆石子,然后取走f个;
5、 f只能是菲波那契数列中的元素(即每次只能取1,2,3,5,8…等数量);
6、 最先取光所有石子的人为胜者;

假设双方都使用最优策略,请判断先手的人会赢还是后手的人会赢。
如果先手的人能赢,请输出“Fibo”,否则请输出“Nacci”

Sample Input
1 1 1 // m n p
1 4 1
0 0 0

Sample Output
Fibo
Nacci

 1 # include <iostream>
 2 # include <cstdio>
 3 # include <cmath>
 4 # include <cstring>
 5 # include <algorithm>
 6 using namespace std ;
 7
 8 const int maxn = 1100 ;
 9
10 int f[maxn],sg[maxn],Hash[maxn];
11 void getsg(int n)
12 {
13     memset(sg,0,sizeof(sg));
14     for(int i=1;i<=n;i++)
15     {
16         memset(Hash,0,sizeof(Hash));
17         for(int j=1;f[j]<=i;j++)
18             Hash[sg[i-f[j]]]=1;
19         for(int j=0;j<=n;j++)
20         {
21             if(Hash[j]==0)
22             {
23                 sg[i]=j;
24                 break;
25             }
26         }
27     }
28 }
29
30 int main ()
31 {
32     f[1] = 1 ;
33     f[2] = 2 ;
34     int i , n ,m , p ;
35     for (i = 3 ; i <= 16 ; i++)
36       f[i] = f[i - 1] + f[i - 2] ;
37     getsg(1000) ;
38
39     while (scanf("%d %d %d" , &m , &n , &p) )
40     {
41         if (m == 0 && n == 0 && p == 0)
42            break ;
43
44         int sum = 0 ;
45         sum = sg[m] ^ sg[n] ^ sg[p] ;
46         if (sum == 0)
47           printf("Nacci\n") ;
48         else
49           printf("Fibo\n") ;
50     }
51
52     return 0 ;
53 }

HDU 1847 (如果找规律法不会 可以用SG)

 1 # include <iostream>
 2 # include <cstdio>
 3 # include <cmath>
 4 # include <cstring>
 5 # include <algorithm>
 6 using namespace std ;
 7
 8 const int maxn = 1100 ;
 9
10 int f[maxn],sg[maxn],Hash[maxn];
11 void getsg(int n)
12 {
13     memset(sg,0,sizeof(sg));
14     for(int i=1;i<=n;i++)
15     {
16         memset(Hash,0,sizeof(Hash));
17         for(int j=1;f[j]<=i;j++)
18             Hash[sg[i-f[j]]]=1;
19         for(int j=0;j<=n;j++)
20         {
21             if(Hash[j]==0)
22             {
23                 sg[i]=j;
24                 break;
25             }
26         }
27     }
28 }
29
30 int main ()
31 {
32     int n ;
33     f[1] = 1 ;
34     int i ;
35     for (i = 2 ; i <= 12 ; i++)
36        f[i] = f[i-1] * 2 ;
37     getsg(1000) ;
38     while (scanf("%d" , &n) != EOF)
39     {
40         int sum = 0 ;
41         sum ^= sg[n] ;
42         if (sum == 0)
43           printf("Cici\n") ;
44         else
45           printf("Kiki\n") ;
46     }
47
48     return 0 ;
49 }

时间: 2024-10-08 10:28:25

SG 大法(Sprague-Grundy函数)的相关文章

【小结】SG生成函数(Grundy函数)

SG生成函数(Grundy函数)小结 转移到子游戏x&y,则 sg[now]=sg[x]∧sg[y] 模板 #include <cstdio> #include <cstring> #include <algorithm> using namespace std; const int MAX = 100007; const int MAX_S = 128; const int MAX_X = MAX; int sg[MAX], stone[MAX_S]; int

博弈SG函数

转自:Sprague-Grundy Function-SG函数--博弈论(3) 公平游戏的Sprague-Grundy定理 公平游戏是一种双人游戏,在游戏中双方都有完整的信息,没有牵涉,任何状态的合法操作对双方来说都是相同的. 一个公平游戏可以抽象地用一个有向无环图来表示,这个图中每个点都对应这一个状态,每条有向边代表从一个状态到另一个状态的合法操作. 我们可以想象一个代币最初放在某个点上,然后两个玩家轮流将其从当前的点移动到它的后继点.当代币移动到汇点时游戏结束,汇点是一个没有出度的点,最后一

hdu1847 Good Luck in CET-4 Everybody! ,巴什博奕,理解SG函数

hdu1847 Good Luck in CET-4 Everybody! 题意: 总共n张牌,双方轮流抓牌,每人每次抓牌的个数只能是2的幂次(即:1,2,4,8,16-),抓完牌,胜负结果也出来了:最后抓完牌的人为胜者.给出n,问先手赢还是后手赢? PS:当然这题可以直接推出 n%3==0必败,否则必胜. //巴什博奕 下面介绍另外一种做法 SG值:一个点的SG值就是一个不等于它的后继点的SG的且大于等于零的最小整数.//同mex()函数 简单点来讲就是当前状态离最近一个必败点的距离. SG(

硬币游戏2&amp;&amp;Cutting Game——Grundy值

Grundy值 当前状态的Grundy值就是除任意一步所能转移到的状态的Grundy值以外的最小非负整数, 以硬币问题一为例,可写成: int init_grundy() { sg[0] = 0; for(int i = 1;i <= x;i++) //递推求前x个SG值 { set<int>st; for(int j = 0;j < k;j++) if(a[j] <= i) st.insert(sg[i - a[j]]); int g = 0; while(st.count

Django 缓存、信号

Reference: http://www.cnblogs.com/lianzhilei/p/6365877.html 缓存 由于Django是动态网站,所有每次请求均会去数据进行相应的操作,当程序访问量大时,耗时必然会更加明显,最简单解决方式是使用:缓存,缓存将一个某个views的返回值保存至内存或者memcache中,5分钟内再有人来访问时,则不再去执行view中的操作,而是直接从内存或者Redis中之前缓存的内容拿到,并返回 Django中提供了6种缓存方式: 开发调试 内存 文件 数据库

巴什博奕

巴什博奕(Bash Game):只有一堆n个物品,两个人轮流从这堆物品中取物,规定每次至少取一个,最多取m个.最后取光者得胜. 显然,如果n=m+1,那么由于一次最多只能取m个,所以,无论先取者拿走多少个,后取者都能够一次拿走剩余的物品,后者取胜.因此我们发现了如何取胜的法则:如果n=(m+1)*r+s,(r为任意自然数,s≤m),那么先取者要拿走s个物品,如果后取者拿走k(≤m)个,那么先取者再拿走m+1-k个,结果剩下(m+1)(r-1)个,以后保持这样的取法,那么先取者肯定获胜.总之,要保

Python学习进程(9)序列

序列是Python中最基本的数据结构.     (1)序列简介: 序列中的每个元素都分配一个数字标明它的位置或索引,第一个索引是0,第二个索引是1,依此类推.序列都可以进行的操作包括索引,切片,加,乘,检查成员. 此外,Python已经内置确定序列的长度以及确定最大和最小的元素的方法. Python有6个序列的内置类型,但最常见的是列表和元组.     (2)列表: 列表可以作为一个方括号内的逗号分隔值出现.列表的数据项不需要具有相同的类型.创建一个列表,只要把逗号分隔的不同的数据项使用方括号括

HDOJ 5724 博弈SG函数

链接: http://blog.csdn.net/tc_to_top/article/details/51958964 题意: n行20列的棋盘,对于每行,如果当前棋子右边没棋子,那可以直接放到右边,如果有就跳过放到其后面的第一个空位子,A先操作,最后谁无法操作则输,给定每行棋子状态,问先手是否必胜 题解: 组合博弈问题,直接sg函数,因为列只有20,可以状压搞,枚举每个状态,找到该状态下可行的操作然后标记 代码: 31 int sg[1 << 21]; 32 int vis[21]; 33

TYVJ 2049 魔法珠 sg函数

题意:链接 方法:sg函数 解析: tyvj的题大部分都没题解啊- - 不过这样貌似会更好?感觉做这的题都需要自己动脑啊- - 虽然嘴上说着好烦然而心里觉得好评? 回归正题 设sg[x]表示数x的sg值,这好像是废话 然后对于读入的a[i],将所有的a[i]的sg值异或起来如果不是零则先手赢反之后手 维护的时候有个坑. 每次求约数的时候,数组要在sg里开,因为如果递归下去的话,全局变量的话会被更改,会被坑死. 然后就是怎么维护了 对于x,先求约数 之后枚举哪个数不取,将其他的异或(或者先都异或起