博弈论及算法实现（一）

最新推荐文章于 2024-10-08 11:17:34 发布

weixin_33804582

最新推荐文章于 2024-10-08 11:17:34 发布

阅读量259

点赞数

原文链接：http://www.cnblogs.com/keam37/p/3403695.html

版权

http://www.cnblogs.com/keam37 keam所有转载请注明出处

在生活中五子棋也是一种先手有必赢策略的游戏，有人会说五子棋先手我也会输啊，所以

博弈论问题都有个类似如“参与者足够聪明”，“两人都不犯错"的前提。

在此前提下，讨论几种常见的博弈情形。

{====================}

一、巴什博弈（Bash Game）

只有一堆n个物品，两个人从轮流中取出（1~m）个；最后取光者胜。

考虑到若n=m+1 那么第一个人不论如何取都不能取胜。

进一步我们发现若 n=k*(m+1)+r; 先取者拿走 r 个，那么后者再拿（1~m）个

n=（k-1）*（m+1）+s；先取者再拿走s 个最后总能造成剩下n=m+1 的局面。

因此，此时先手有必赢策略。

相对应的，若n=k*(m+1) 那么先取者必输。

因此我们可以写出对应的程序（默认 n m都大于0）

1 int Bash_Game(int n,int m)//是否先手有必赢策略
2 {
3     if (n%(m+1)!=0) return 1;
4     return 0;
5 }

到了这，不如尝试一下当有N堆，每堆有Mi>0个物品，依旧是两个人来取该怎么判断？

先考虑取的最大数目无上限即可以把一堆全部取完的情形

{====================================}

二、尼姆博弈(Nimm Game)

正如上述.

从巴什博弈我们知道一个当情形对应一种状态,而由一个状态只能变为另一种状态时能很轻易地判断是否先手有必赢策略

那么怎么样才能在尼姆博弈里找到这样的状态呢?

如果把n堆抽象为n个非负整数,再将n个整数转化为二进制,然后对n个二进制数按位相加(不进位),若每一位相加都为偶数,

那么称这个状态为偶状态,否则称它为奇状态.

可以证明:任何一个偶状态在其中一个数变小后一定成为奇状态,而一个奇状态一定可以通过改变一个数变成偶状态.

前一点很显然,因为一个数变小至少有一位发生改变,这一位就改变了原来的偶状态.

对于后一点,对于一个从高位到低位某一位和为奇的奇状态,必定有一个数的二进制表示在此位为1,对于后面的较低位和为奇的情况,只要把这个数对应位取反即可得到一个偶状态.

到此,成功的构造了两个可以转换的状态!!!

那么对于n堆物品,只要判断它是否是奇状态就可以判断是否先手有必赢策略.

但是求每个数的二进制表示略显麻烦,可以用更好的办法,也是我偏爱的位运算.

XOR 和判断:

如果有奇数个二进制数在第K位为1 那么在这一位上的和为奇,同样的，偶数个1和为偶.

很明显位运算xor满足我们的要求,偶数个1异或和为0,奇数个为1;

由此,终于可以给出算法

1 int Nimm_Game(int n)//假设n个数存在数组f[]中,有必胜策略返回1
2 {
3     int flag=0;
4     for(int i=1;i<=n;i++)
5     flag^=f[i];
6     if(flag) return 1;
7     return 0;
8 }

{====================================}

但是我遇到过n非常大，且每一堆的物品数为连续的整数的情况

对此我们要考虑连续非负整数的异或和问题

记 f(x, y) 为x到y的所有整数的异或值。

f[1,n]=f[0,n];

当n为2^k-1(2的K次方减一)时;

0 到 2^k-1 共2^k个数等于∑C（n，i）

可以看做在k个位置中放入i个0，最后求和

同时可以看做在空格位置中放入i个1；最后求和

即在每一位上1个0的个数都相等，每个位上有2^(k-1)个1,当k>=2时 1的个数为偶数;

而我们已经知道偶数个1的异或和为0

所以 f[0, 2^k - 1] = 0 (k >= 2)

对 f[0, n] (n>=4) 设n的最高位1是在第k位(k >= 2)，

f[0, n] = f[0, 2^k - 1] xor f[2^k, n] = f[2^k, n]

对2^k到n这n+1-2^k个数，最高位(第k位)共有 m = n+1-2^k 个1,

2^k总是偶数,因此,当n为奇数时，m是偶数，f[0, n] = f[2^k, n] = f[0, n - 2^k]

当n为偶数(n - 2^k)总是偶数，所以：f[0,n] = f[0,n-2^k-2^(k-1)...-2^2](只到3是因为k>=2)

此时只剩下两位是我们需要的我们可以用（n & 3）很快得到后两位

由于n是偶数所以（n & 3）只可能得到 1 或 3；

1对应二进制数（01）所以是奇数个1 此时f [0,n]=1;

3对应二进制数（11）此时f[0,n]=0;

当n为偶数时，m是奇数，因而 f[0, n] = f[2^k, n] = f[0, n - 2^k] xor 2^k

可得：f[0, n] = f(0, n & 3) xor 2^k xor n[k]*2^(k-1) xor ....n[2]*2^2 n[k] 为 n的二进制数的第k位；

很明显当n为偶数时 f[0,n]的二进制从最高位到第3位（如果不止3位）跟n的二进制数从高位到第三位相同；

此时只需要判断第二位

n & 3=0对应后二位为（00）此时 f[0,n]=n;

n & 3=2对应后二位为（10）此时 f[0,n]=n+1;

综上所述：

n n % 4 == 0

f[1, n] = f[0, n] = 1 n % 4 == 1

n +1 n % 4 == 2

0 n % 4 == 3

f[x,y] = f[0,y] xor f[0,x-1] .(x>0)

对应的代码在这

 1 //读入n，表示有从物品数分别1到n的n堆物品，假设n个数存在数组f[]中
 2 int xor_n(int n)//从1到n的异或和
 3 {
 4      int t = n & 3;
 5      if (t & 1) return t / 2 ^ 1;
 6      return t / 2 ^ n;
 7 }
 8 int Nimm_Game(int n)//有必胜策略返回1
 9 {
10     int flag=0;
11     for(int i=1;i<=n;i++)
12     flag^=xor_n(f[i]);
13     if(flag) return 1;
14     return 0;
15 }