通常的Nim游戏的定义是这样的:有若干堆石子,每堆石子的数量都是有限的,合法的移动是“选择一堆石子并拿走若干颗(不能不拿)”,如果轮到某个人时所有的石子堆都已经被拿空了,则判负(因为他此刻没有任何合法的移动)。

结论:(不给证明,可去搜证明)

对于一个Nim游戏的局面(a1,a2,...,an),它是必败态当且仅当a1^a2^...^an=0,其中^表示异或(xor)运算。怎么样,是不是很神奇?我看到它的时候也觉得很神奇,完全没有道理的和异或运算扯上了关系。但这个定理的证明却也不复杂,基本上就是按照两种position的证明来的。

int N,A[maxn]
void solve( )
{
int x = ;
for(int i= ; i<N ; i++)
x^=A[i];
if(x!=)
puts("A");
else
puts("B");
return ;
}

现在来实战分析Nim

poj 1704

题意】

从左到右有一排石子,给出石子所在的位置。规定每个石子只能向左移动,且不能跨过前面的石子。最左边的石子最多只能移动到1位置。每次选择一个石子按规则向左移动,问先手是否能赢。

分析:如果我们将棋子两两看成一对,之间的距离表示Nim里面小石头堆的数量,那这就是个Nim.那这里就有个疑问了,那奇数的时候怎么办了呐? 很简单拿第一个棋子的距离与 1 当成是一个堆其他的不变就好了,那问题又来了,为什么,可以这样转换呢? 这个道理很简单我棋子向左移动,是不是距离就减少了,那不就是相当于拿走石头吗,等等!聪明仔可能发现了问题了,一个棋子左移动,那肯定是有一堆的石头的数量增加了呀 , 这不符合Nim呀!想想也是,不过如果我们将其看做是对手不讲道义,破坏规则,进行加石头 。 那我们只要将对手增加石头的部分减回去,不就是原来的状态了吗。神奇不,注意需要排序,超坑

#include<stdio.h>
#include<algorithm>
using namespace std;
const int maxn = ;
int n,p[maxn],val[maxn];
int main( )
{
int t;
scanf("%d",&t);
while(t--)
{ scanf("%d",&n);
for(int i= ; i<=n ; i++)
scanf("%d",&p[i]) ;
sort(p+,p++n);
int cnt=;
if(n%==)
{
for(int i= ; i<=n ; i+=)
{
val[cnt++]=p[i]-p[i-]-;
}
}
else
{
val[cnt++]=p[]-;
for(int i= ; i<=n ; i+=)
{
val[cnt++]=p[i]-p[i-]-;
}
}
int x=;
for(int i= ; i<cnt ; i++)
{
x^=val[i];
//printf("%d ",val[i]);
} if(x!=)
puts("Georgia will win");
else
puts("Bob will win");
}
return ;
}

题目:

A 和 B 在玩游戏,给k个数字a1 , a2 , ... , ak ; 一开始有n堆硬币 , 每一堆有 xi 枚 。 A 和 B 轮流选出一堆硬币 ,从中取出硬币 , 每次所取硬币的枚数一定在a1 , a2 , a3 ,...., ak ; 里面,A先取,取光硬币的一方获胜 。

分析:这里引入一个概念 Grundy , 利用这个东西,很多游戏都可以转换为Nim>

只有一堆石头的情况

int grundy(int x){
集合S={};
for(j=:k){
if(a_j<=x) 将grundy(x-a_j)加到S集合中
}
return 最小的不属于S的非负整数
}

rundy值:除(任意一步所能转移到 的状态  的Grundy值 )以外的最小非负整数,这样的Grundy值,和Nim中的一个石子堆类似,有如下性质:

mex{0,1,2}=3;mex{ 1, 2}=0 ; mex{ 2, 3}=1

1.Nim中有x颗石子的石子堆,能转移成有0,1,……,x-1堆石子的石子堆

2.从Grundy值为x的状态出发,可以转移到Grundy值为0,1,……,x-1的状态。

Nim:

所有石子堆的石子数xi的XOR

x1 XOR x2 XOR …XOR xk

为0必败,为1必胜

Grundy值等价于Nim中石子数,所以对于Grundy的情况:

所有硬币堆的Grundy的值

grundy(x1)  XOR  grundy(x2)  XOR …… XOR grundy(xn)

为0必败,为1必胜。

#include <cstdio>
#include <set>
#include <algorithm>
#define maxn 105
#define maxk 105
using namespace std;
int N,K,X[maxn],A[maxk];
int grundy[maxn+]; void solve(){
//轮到自己剩0的时候必败
grundy[]=; //计算grundy
int max_x=*max_element(X,X+N);//找最大值
for(int j=;j<=max_x;j++){
set<int> s; //存储当前的状态
for(int i=;i<K;i++){
if(A[i]<=j) s.insert(grundy[j-A[i]]);
}
int g=; //寻找当前状态的最小排斥
while(s.count(g)!=) g++;
grundy[j]=g;
}
int ans=;
for(int i=;i<N;i++) ans^grundy[X[i]]; if(ans!=) puts("Alice");
else puts("Bob"); } int main(){
scanf("%d%d",&N,&K);
for(int i=;i<N;i++){
scanf("%d",&A[i]);
}
for(int j=;j<K;j++){
scanf("%d",&X[j]);
}
solve(); return ;
}

poj2311

题意:给定一张n*m的矩形纸片,定义一个决策为:从所有的纸片中任选一张(初始时仅有1张),沿平行于矩形边剪切,将原矩形分割为2个较小的矩形,双方轮流决策,先剪出1*1纸片者获胜,问先手是否必胜.

经典的组合游戏问题,首先一张纸片经过一次裁剪变为两张,即一个游戏变为两个子游戏,那么子游戏和的sg值=子游戏A的sg值异或子游戏B的sg值,即sg[A∪B] = sg[A] xor sg[B].

很容易想出将(1,1)作为终止条件.但是考虑组合游戏本质为在一张有向图中从S到T双方轮流移动,但此题中只要出现一张(1,1)游戏便结束,以(1,1)作为终止条件可能会产生有向图并没有走完但游戏已经结束的情况.

稍加分析可知,某一方某一次决策之后产生了一张(1,k)的纸片,那么这一方必败(因为另一方将(1,k)裁剪为(1,1)&(1,k-1)可以直接获胜),那么对于(i,j)的纸片,最优决策一定会避免剪出(1,k)的纸片,当纸片无论怎么剪都会出现(1,k),称这张纸片为"不可剪",当某纸片不可剪时,则会选择其他纸片进行裁剪,易证,存在某一个状态,所有的纸片都是"不可剪".

不可剪状态的纸片有:(2,2), (2,3), (3,2), (3,3).

那么只需要将这些纸片作为终止条件问题就解决了.

#include<stdio.h>
#include<set>
#include<string.h>
using namespace std ;
int mem[][];
int grundy(int w , int h)
{
if(mem[w][h] != -)
return mem[w][h] ; set<int> s;
for(int i= ; w-i>= ; i++)
s.insert(grundy(i,h)^grundy(w-i,h));
for(int i= ; h-i>= ; i++)
s.insert(grundy(w,i)^grundy(w,h-i));
int res = ;
while(s.count(res))
res++;
return mem[w][h]=res ;
}
int main( )
{
int w,h;
memset(mem,-,sizeof(mem));
while(scanf("%d%d",&w,&h)!=EOF)
{ if(grundy(w,h) != )
puts("WIN");
else
puts("LOSE");
}
}
<span style="font-size:18px;">#include <iostream>
#include <cstdio>
#include <cstring> using namespace std; #define rep(i, l, r) for (int i = l; i <= r; i++)
#define REP(i, l, r) for (int i = l; i >= r; i--)
#define MAXN 1010 int n, m, sg[MAXN][MAXN];
bool xx[MAXN]; inline int mex() {int i; for (i = ; xx[i]; ++i); return i;} inline int getsg(int n, int m) {
if (~sg[n][m]) return sg[n][m];
memset(xx, , sizeof(xx));
rep(i, , n-) xx[getsg(i, m) ^ getsg(n-i, m)] = ;
rep(i, , m-) xx[getsg(n, i) ^ getsg(n, m-i)] = ;
return sg[n][m] = mex();
} int main() {
memset(sg, -, sizeof(sg));
sg[][] = sg[][] = sg[][] = sg[][] = sg[][] = ;
while (scanf("%d%d", &n, &m) != EOF)
cout << (getsg(n, m) ? "WIN" : "LOSE") << endl; return ;
}
</span>

深度了解问题

最新文章

  1. Standard C 语言标准函数库介绍
  2. 悲剧啊!Mysql的上古BUG!!!
  3. Git 安装
  4. 源码方式安装mysql5.5
  5. Oracle SQL基本操作
  6. LeetCode &quot;Binary Tree Level Order Traversal II&quot; using DFS
  7. redis使用日志(一) 安装,调试
  8. 关于Liferay所有的能够进行自定义和扩展的东西的总结
  9. Cocos2d-JS坐标系
  10. 【6】和作为连续序列s
  11. C++获取本机IP等信息
  12. using Newtonsoft.Json;
  13. 机器学习入门11 - 逻辑回归 (Logistic Regression)
  14. Software Engineer Title Ladder
  15. 解决Could not open Hibernate Session for transaction; nested exception is java.lang.NoClassDefFoundError: org/hibernate/engine/transaction/spi/TransactionContext
  16. JDBC 编程步骤
  17. Qt addStretch()详解
  18. Shell编程-05-Shell中条件测试与比较
  19. LeetCode-63.不同路径Ⅱ
  20. Codeforces 493 E.Devu and Birthday Celebration

热门文章

  1. Oracle 12C 新特性之级联truncate
  2. bzoj 3996 线性代数 —— 最大权闭合子图
  3. org.springframework.web.client.HttpClientErrorException: 400 null
  4. Hyperledger fablic 0.6 在centos7环境下的安装与部署
  5. SQL SERVER 中的*=和=*
  6. JINKENS
  7. MyBatis总结(1)
  8. ORM查询相关
  9. 3.JasperReports学习笔记3-在浏览器生成PDF文件
  10. Centos6.5安装JDK8教程(一)