取石子游戏

取石子游戏

Time Limit: 1000MS		Memory Limit: 10000K
Total Submissions: 25176		Accepted: 7961

Description

有两堆石子，数量任意，可以不同。游戏开始由两个人轮流取石子。游戏规定，每次有两种不同的取法，一是可以在任意的一堆中取走任意多的石子；二是可以在两堆中同时取走相同数量的石子。最后把石子全部取完者为胜者。现在给出初始的两堆石子的数目，如果轮到你先取，假设双方都采取最好的策略，问最后你是胜者还是败者。

Input

输入包含若干行，表示若干种石子的初始情况，其中每一行包含两个非负整数a和b，表示两堆石子的数目，a和b都不大于1,000,000,000。

Output

输出对应也有若干行，每行包含一个数字1或0，如果最后你是胜者，则为1，反之，则为0。

Sample Input

2 1
8 4
4 7

Sample Output

0
1
0

我当时做的思路：

此题如果用常规的博弈算法，比如极大极小、估值函数之类的方法来做，会举步维艰。因为内容限制为10MB，但是数据值可大至1,000,000,000（约等于1G），所以不可能建一个如此大的数组。我开始观察必赢/输的数值规律：

发现在轮到你时对于给定的石子数（x,y），可以确定是必赢还是必输。以1表示必赢，0表示必输。

（0,1）=1，........，（0，n）=1

（1,1）=1，（1,2）=0，（1，3）=1，........ ，（1，n）=1     //我可以拿其中一堆中n-2个石子，给对手留下（1,2）的局面

（2,2）=1，（2,3）=1，.........，（2，n）=1      //我可以拿其中一堆n-1个石子，给对手留下（2,1）的局面

（3,3）=1，（3,4）=1，（3,5）=0，（3,6）=1，.........，（3,n）=1   //拿n-5个石子，留下（3,5）的局面

..................... 

（4,7）=0；（6,10）=0；（8,13）=0；（9,15）=0 。。。。。。。。。。。

其中必输局面的规律可以得出来：

（1,2）、（3,5）、（4,7）、（6,10）、（8,13）、（9,15）、（11,18）、（12,20）。。。

（x，y）中 y-x 的值是递增1的，而且x和y是不重复的自然数。

到这里就有些眉目了，只要用一个stone[i]数组保存（x，y）中的x值（因为y=x+i），再用一个布尔数组记录自然数序列中的值是否已用。这样看来，问题变得很简单了。

我的代码如下：

#include <iostream>
using namespace std;

#define LEN 1000000

unsigned int *s;

bool exist[LEN];
int maxDelta=0;

int findNext(int x)
{
	x++;
	while(1)
	{
		if(exist[x])x++;
		else return x;
	}
}

void Extend(int del)
{
	int i,next;
	next=s[maxDelta];
	for(i=maxDelta+1;i<=del;i++)
	{
		next = findNext(next);
		exist[next]=true; 
		exist[next+i]=true;
		s[i]=next;
	}
}

int main()
{
	unsigned int x,y;
	int del;

	memset(exist,0,LEN);
	exist[0]=true;
	s=(unsigned int *)malloc(sizeof(unsigned int)*LEN);
	memset(s,0,sizeof(unsigned int)*LEN);
	while(cin>>x && cin>>y)
	{
		del=abs(x-y);
		if(del==0){cout<<1<<endl;continue;}
		if(del>maxDelta)
		{
			if(x<LEN && y<LEN)
			{
				if(exist[x] || exist[y])
				{	cout<<1<<endl;continue;}
			}
			Extend(del);//向上处理
			maxDelta=del;
		}
		if((s[del]==x && s[del]+del==y) ||(s[del]==y && s[del]+del==x))cout<<0<<endl;
		else cout<<1<<endl;	
		//for(int i=0;i<=maxDelta;i++)fout<<s[i]<<","<<s[i]+i<<endl;
	}
	return 0;
}

提交上去，出现Runtime Error错误，回头看代码意识到应该是bool exist[]的大小的问题，因为要记录下0~大约min(x,y)/2个数的使用状态，还是容易内存超量使用。

在网上搜了下解题报告，发现问题可以变得更简单，这是典型的威佐夫博奕算法。现将算法原理抄录如下：

威佐夫博奕（Wythoff Game）：

有两堆各若干个物品，两个人轮流从某一堆或同时从两堆中取同样多的物品，规定每次至少取一个，多者不限，最后取光者得胜。
    这种情况下是颇为复杂的。我们用（ak，bk）（ak ≤ bk ,k=0，1，2，...,n)表示两堆物品的数量并称其为局势，如果甲面对（0，0），那么甲已经输了，这种局势我们称为奇异局势。前几个奇异局势是：（0，0）、（1，2）、（3，5）、（4，7）、（6，10）、（8，13）、（9，15）、（11，18）、（12，20）。
    可以看出,a0=b0=0,ak是未在前面出现过的最小自然数,而 bk= ak + k，奇异局势有
如下三条性质：

    1。任何自然数都包含在一个且仅有一个奇异局势中。
    由于ak是未在前面出现过的最小自然数，所以有ak > ak-1 ，而 bk= ak + k > ak-1 + k-1 = bk-1 > ak-1 。所以性质1。成立。
    2。任意操作都可将奇异局势变为非奇异局势。
    事实上，若只改变奇异局势（ak，bk）的某一个分量，那么另一个分量不可能在其他奇异局势中，所以必然是非奇异局势。如果使（ak，bk）的两个分量同时减少，则由于其差不变，且不可能是其他奇异局势的差，因此也是非奇异局势。
    3。采用适当的方法，可以将非奇异局势变为奇异局势。

    假设面对的局势是（a,b），若 b = a，则同时从两堆中取走 a 个物体，就变为了奇异局势（0，0）；

如果a = ak ，b > bk，那么，取走b  - bk个物体，即变为奇异局势；

如果 a = ak ，  b < bk ,则同时从两堆中拿走 ak - ab - ak个物体,变为奇异局势（ ab - ak , ab - ak+ b - ak）；

如果a > ak ，b= ak + k,则从第一堆中拿走多余的数量a - ak 即可；

如果a < ak ，b= ak + k,分两种情况，第一种，a=aj （j < k）,从第二堆里面拿走 b - bj 即可；第二种，a=bj （j < k）,从第二堆里面拿走 b - aj 即可。

    从如上性质可知，两个人如果都采用正确操作，那么面对非奇异局势，先拿者必胜；反之，则后拿者取胜。

    那么任给一个局势（a，b），怎样判断它是不是奇异局势呢？我们有如下公式：
    ak =[k（1+√5）/2]，bk= ak + k  （k=0，1，2，...,n 方括号表示取整函数)
奇妙的是其中出现了黄金分割数（1+√5）/2 = 1。618...,因此,由ak，bk组成的矩形近似为黄金矩形，由于

2/（1+√5）=（√5-1）/2，可以先求出j=[a（√5-1）/2]，若a=[j（1+√5）/2]，那么a = aj，bj = aj + j，若不等于，那么a = aj+1，bj+1 = aj+1+ j + 1，若都不是，那么就不是奇异局势。然后再按照上述法则进行，一定会遇到奇异局势。

如此以来，就变成对于局势（x,y）利用公式判断是否满足奇异局势。

判断 x ?= ((y-x)*(sqrt(5)+1)) / 2 向下取整

代码如下：

#include <iostream>
#include <math.h>
using namespace std;

int main()
{
	int x,y,t;
	while(cin>>x && cin>>y)
	{
		if(x>y)
		{
			t=x;
			x=y;
			y=t;
		}
		t=floor((y-x)*(sqrt(5.0)+1.0))/2;
		if(t==x)cout<<0<<endl;
		else cout<<1<<endl;
	}
	return 0;
}