博弈论基础与习题(未完)

转自博弈论SG函数

一.巴什博奕(Bash Game):

首先我们来玩一个比较古老的报数游戏。A和B一起报数,每个人每次最少报一个,最多报4个。轮流报数,看谁先报到30.

如果不知道巴什博弈的可能会觉得这个是个有运气成分的问题,但是如果知道的人一定知道怎样一定可以赢。

比如A先报数的话,那么B一定可以赢(这里假定B知道怎么正确的报数)

B可以这样报数,每次报5-k(A)个数,其中k(A)是A报数的个数这样的话没一次

两人报完数之后会变成5 10 15 20 25 30这样是不是B一定会赢呢?是不是有一种被欺骗的感觉呢?好吧下面我们来看看这个原理。我们先看下一个一眼就能看出答案的例子 比如说我们报到5(4+1),每次报最多报4个,最少报1个.那么是不是后者一定可以赢呢?答案是肯定的。好了到这巴什博弈的精髓基本就OK了。

那么如果我们要报到n+1,每次最多报n个,最少报1个的话,后者一定能够赢。

现在我们需要报数到n,而每次最多报数m个,最少报数1个.我们可以化成这样

n = k*(1+m)+r(0 <= r <= m)这样的话如果r不等于0那么先手一定会赢,为什么呢?首先先手报r个,那么剩下k倍(1+m)个数,那么我们每次报数1+m-k(B)个数就一定能保证最后剩下1+m个,那么就到了上面我们说的那个了,先手就一定会赢,如果r=0那么后手一定会赢,道理一样的。

到这巴什博弈也就介绍完了,知道这个道理之后我们也可以去骗小朋友了。-_-//

#include<iostream>  
#include<string>  
#include<cstring>  
#include<cstdio>  
#include<algorithm>  
#define CLR(arr, val) memset(arr, val, sizeof(arr))  
using namespace std;  
  
int main()  
{  
    //freopen("Input.txt", "r", stdin);  
    int N, num, limit;  
    scanf("%d", &N);  
    while(N--)  
    {  
        scanf("%d%d", &num, &limit);  
        if(num % (limit + 1) != 0) //必胜局面  
            printf("Win
");  
        else  
            printf("Lose
");  
    }  
    return 0;  
}

二.威佐夫博奕(Wythoff Game):

   这种博弈比前面一种要稍微复杂一点。我们来看下下面这个游戏。

   有两堆火柴棍,每次可以从某一堆取至少1根火柴棍(无上限),或者从两堆取相同的火柴棍数。最后取完的是胜利者。好了,如果你不知道这个博弈定理,对于小数目的火柴棍数,可能还能推出来,但是如果火柴棍数一多,就不行了。看了下面的这个介绍,你也会有一种被骗的感觉。

   首先我们知道两堆火柴是没有差别的,也就是说第一堆有a根,第二堆有b根和第一堆有b根,第二堆有a根是一样的结果。

   我们用一个二维的状态(a,b)来记录当前剩下的火柴数,表示第一堆剩下a根火柴,第二堆剩下b根火柴。同样我们假设两个人的编号是A和B,且A先取。

那么如果某个人遇到了这样的状态(0,0)那么也就是说这个人输了。这样的状态我们叫做奇异状态,也可以叫做失败态。

那么接下来的几个失败态为(1,2),(3,5),(4,7),(6,10),(8,13)……

我们用a[i]表示失败态中的第一个,b[i]表示失败态中的第二个.(i从0开始).

那么我们可以看到b[i] = a[i]+i;(i >= 0),a[i]是前面的失败态中没有出现过的最小的整数

下面我们可以得到三个基本的结论。

  1.每个数仅包含在一个失败态中

  首先我们知道a[k]是不可能和前面的失败态中的a[i],b[i]重复的(这点由a[i]的得到可以知道)

b[k] = a[k]+k > a[k-1]+k>a[k-1]+k-1+1>a[k-1]+(k-1) = b[k-1]>a[k-1]这样我们知道每个数仅在一个失败态中。

  2.每个失败态可以转到非失败态。

 加入当前的失败态为(a,b),那么如果我们只在一堆中取的话,肯定会变成非失败态(这点由第一点可以保证),如果从两堆同时取的话,由于每个失败态的差是不一样的,所以也不可能得到一个失败态。也就是说一个失败态不管你怎么取,都会得到一个非失败态。

   3.每个非失败态都可以转到一个失败态

对于这个结论,首先我们要知到每个状态(a,b)要么a = a[i],要么b = b[i].(每个数都出现在一个失败态中),下面我们分两种情况来讨论

   I.a = a[i].如果b = a的话那么一次取完就变成了(0,0).如果b > b[i]的话,那么我们从第二堆中取走b-b[i]就变成了一个失败态。如果b < b[i].那么我们从两堆中同时取走a-a[b-a[i]]这样得到失败态(a[b-a[i]],a[b-a[i]]+b-a[i])(a[i] = a)

   II.b = b[i].如果a > a[i]那么我们从第一堆中取走a-a[i]根火柴.

              如果a < a[i].这里又分两种情况。第一是a = a[k](k < i)

那么我们从第二堆取走b - b[k]就行了。

第二是a = b[k]这样的话由于两堆火柴是没有区别的,所以我们把b变成a[k]就行了,也即是从第二堆火柴中取走b - a[k]就变成了失败态

至于怎么判断一个状态是否是失败态.我们可以用下面的方法来判断(本人暂时还不会证明)

  a[i] = [i*(1+√5)/2](这里的中括号表示向下取整)   b[i] = a[i]+i;

  那么这就是一个失败态

代码如下:
 

#include<iostream>
#include<string>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define CLR(arr, val) memset(arr, val, sizeof(arr))
using namespace std;
 
int main()
{
	//freopen("Input.txt", "r", stdin);
	int num1, num2, tmp; //第一堆剩的数量为num1,第二堆剩num2
	while(scanf("%d%d", &num1, &num2) != EOF)
	{
		if(num1 > num2)
			swap(num1, num2); 
		tmp = floor((num2 - num1) * (1 + sqrt(5.0)) / 2.0); //黄金分割
		if(tmp == num1)	printf("Lose
"); //奇异局势必败
		else	printf("Win
");
	}
	return 0;
}        

三.尼姆博奕(Nimm Game):

指的是这样的一个博弈游戏,目前有任意堆石子,每堆石子个数也是任意的,双方轮流从中取出石子,规则如下:
1)每一步应取走至少一枚石子;每一步只能从某一堆中取走部分或全部石子;
2)如果谁取到最后一枚石子就胜。
也就是尼姆博弈(Nimm Game)。
必败局面:也叫奇异局势。无论做出何出操作,最终结果都是输的局面。必败局面经过2次操作后,可以达到另一个必败局面。
必胜局面:经过1次操作后可以达到必败局面。
即当前局面不是必败局面就是必胜局面,而必胜局面可以一步转变成必败局面。
最终状态:
(1)最后剩下一堆石子;(必胜局面)
(2)剩下两堆,每堆一个;(必败局面)
(3)当石子剩下两堆,其中一堆只剩下1颗,另一堆剩下多于n颗石子时,当前取的人只需将多于1颗的那一堆取出n-1颗,则局面变为刚才提到的必败局面。(必胜局面)
判断当前局势是否为必胜(必败)局势:
1)把所有堆的石子数目用二进制数表示出来,当全部这些数按位异或结果为0时当前局面为必败局面,否则为必胜局面;
2)在必胜局面下,因为所有数按位异或的结果是大于零的,那么通过一次取,将这个(大于其它所有数按位异或的结果的)数下降到其它所有数按位异或的结果,这时局面就变为必败局面了。
定理:一组自然数中必然存在一个数,它大于等于其它所有数按位异或的结果。
证明:原命题等价于,设a1^a2^... ^an=p,p≠0时,必存在k,使得ak^p<ak(当p=0时,对于任意的k,有ak^p=ak)。
设p的最高位是第q位,则至少存在一个k,使得ak的第q位也是1,而ak^p的第q位为0,所以ak^p<ak
    补缀一点,(a^b)^b=a^(b^b)=a^0=a,所以ak^p相当于“其它所有数按位异或的结果”。

例1:2 45 45


45^45=0,45和45的异或等于0。
例 2:3 3 6 9

局势(3,6,9)因为3^6^9不等于0,所以这是一个必胜局势。
 3 011

^6 110

 5 101 
即从第3堆中的9个中取走9-5=4个,则(3,6,9)->(3,6,5),3^6^5=0,故(3,6,5)为奇异局势,即从必胜局势转变成必败局势。

代码如下:
 

#include<iostream>
using namespace std;
int temp[ 20 ]; //火柴的堆数
 
int main()
{
	int i, n, min;
	while( cin >> n )
	{
		for( i = 0; i < n; i++ )
			cin >> temp[ i ]; //第i个火柴堆的数量
		min = temp[ 0 ];
		for( i = 1; i < n ; i++ )
			min = min^temp[ i ]; //按位异或
		if( min == 0 )
			cout << "Lose" << endl; //输
		else
			cout << "Win" << endl; //赢
	}
	return 0;
}

讲完三种博弈,接下来就是精髓

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 的函数值。

void getsg(int n)
{
    bool mex[11111];
    for(int i=0;i<n;i++)
    {
        memset(mex,false,sizeof(mex));
        		for (int j = 1; j <= k; ++j)//k为fa数组元素个数
        		{
        			if (i - fa[j] >= 0)
        			{
        				mex[sg[i - fa[j]]] = true;
        			}
        		}
        for(int j=0;j<=n;j++)
            if(mex[j]==false)
            {
                sg[i]=j;
                break;
            }
    }
}

例题

A:LightOJ-1186 Incredible Chess:题意:棋盘上每一纵列有黑白两个棋子, 黑色一定在白色上方,两个人一人各执一色, 白色先走,每次操作可以把自己的颜色向上或者向下移动任意格, 但是不能超出范围也不能越过另一色的棋, 无法操作的一方败。思路:nim博弈, 将两个棋子之间的距离看作是堆中的石子的数量, 然后就转化成一个简单的nim博弈。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <string>
#define maxn 10005
#define INF 0x3f3f3f3f
typedef long long ll;
using namespace std;
int main()
{
    int n;
    int t,T=1;
    cin>>t;
    int a[110],b[110];
    while(t--){
        cin>>n;
        int ans=0;
        for(int i=0;i<n;i++){
            cin>>a[i];
        }
        for(int i=0;i<n;i++){
            cin>>b[i];
        }
        for(int i=0;i<n;i++)
            ans^=(b[i]-a[i]-1);
        if(ans==0)
            printf("Case %d: black wins
",T++);
        else
            printf("Case %d: white wins
",T++);
    }
    return 0;
}

B:LightOJ-1192 Left Right:给定n对点,每一对点中第一个点向右走,第二个点向左右,不能跨越...谁不能走谁输...问最后谁赢。

和上一题同样的思路;

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <string>
#define maxn 10005
#define INF 0x3f3f3f3f
typedef long long ll;
using namespace std;
int main()
{
    int n;
    int t,T=1;
    cin>>t;
    while(t--){
        cin>>n;
        int ans=0;
        for(int i=0;i<n;i++){
            int x,y;
            cin>>x>>y;
            ans^=(y-x-1);
        }
        if(ans==0)
            printf("Case %d: Bob
",T++);
        else
            printf("Case %d: Alice
",T++);
    }
    return 0;
}

C:LightOJ-1253 Misere Nim:在K堆石子的任意堆中任取>0个石子,Alice为先手,Bob为后手,先取完者输掉比赛(注意这里)。

思路没什么说的就是改成先取完输。具体怎么写模板改一下就行了。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <string>
#define maxn 10005
#define INF 0x3f3f3f3f
typedef long long ll;
using namespace std;
int main()
{
    int n;
    int t,T=1;
    cin>>t;
    while(t--){
        cin>>n;
        int ans=0,f=0;
        for(int i=0;i<n;i++){
            int x;
            cin>>x;
            ans^=x;
            if(x>1)
                f=1;
        }
        if(f==0){
            if(n%2)
                printf("Case %d: Bob
",T++);
            else
                printf("Case %d: Alice
",T++);
        }
        else{
            if(ans==0)
                printf("Case %d: Bob
",T++);
            else
                printf("Case %d: Alice
",T++);
        }
    }
    return 0;
}

D:LightOJ-1393 Crazy Calendar:题意:给定n*m的一个方格,每个方格含x个石子,Pinky 先手,每次向右或者下移动,     谁不能移动谁输。思路:如果该方格中的一个点的坐标与n+m的奇偶性相同,则不需要考虑该点。因为若该点与n+m奇偶性相同,设该点坐标为(i,j)(0<=i<=n,0<=j<=m)则(n+m)-(i+j),即到达终点的步数为偶数,也就是说先手怎么移动后手就怎么移动,也就是说等于没移动。所以我们只需要考虑与n+m不同的,说白了也就是意思n+m和i+j一个奇数一个偶数,也就是说,这些点距离终点只需要也只能移动一次,每次可以从一个方格中取任意个石子移动,这就转化为了尼姆博弈中的取石子问题,因为尼姆博弈中可以随意取石子。异或为0为必败态一定会转化为必胜态,异或为1为必胜态可以转化为必败态,具体详见下最基础的尼姆博弈--取石子问题,就ok了。
 

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <string>
#define maxn 10005
#define INF 0x3f3f3f3f
typedef long long ll;
using namespace std;
int main()
{
    int n,m;
    int t,T=1;
    cin>>t;
    while(t--){
        cin>>n>>m;
        int ans=0;
        int tem=(n+m)%2;
        for(int i=0;i<n;i++){
            for(int j=0;j<m;j++){
                int x;
                cin>>x;
                if((i+j)%2!=tem)
                    ans^=x;
            }
        }
        if(ans==0)
            printf("Case %d: lose
",T++);
        else
            printf("Case %d: win
",T++);
    }
    return 0;
}

E:HDU-2516 取石子游戏:题意:1堆石子有n个,两人轮流取.先取者第1次可以取任意多个,但不能全部取完.以后每次取的石子数不能超过上次取子数的2倍。取完者胜.思路: 斐波那契博弈,n为斐波那契数时先手必败,记录一下斐波那契数即可。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <string>
#define maxn 10005
#define INF 0x3f3f3f3f
typedef long long ll;
using namespace std;
int main()
{
    int n,fib[50];
    fib[0]=2;
    fib[1]=3;
    for(int i=2;i<=50;i++)
        fib[i]=fib[i-1]+fib[i-2];
    while(cin>>n&&n){
        int flag=1;
        for(int i=0;i<=50;i++){
            if(fib[i]==n)
                flag=0;
            if(fib[i]>n)
                break;
        }
        if(flag)
            cout<<"First win"<<endl;
        else
            cout<<"Second win"<<endl;
    }
    return 0;
}
原文地址:https://www.cnblogs.com/shmilky/p/14089032.html