51Nod_1185 威佐夫游戏 V2

                                          51Nod_1185 威佐夫游戏 V2

                                         http://www.51nod.com/Challenge/Problem.html#!#problemId=1185

 

 

题目

有2堆石子。A B两个人轮流拿，A先拿。每次可以从一堆中取任意个或从2堆中取相同数量的石子，但不可不取。拿到最后1颗石子的人获胜。假设A B都非常聪明，拿石子的过程中不会出现失误。给出2堆石子的数量，问最后谁能赢得比赛。例如：2堆石子分别为3颗和5颗。那么不论A怎样拿，B都有对应的方法拿到最后1颗。

输入

第1行：一个数T，表示后面用作输入测试的数的数量。（1 <= T <= 10000)；第2 - T + 1行：每行2个数分别是2堆石子的数量，中间用空格分隔。(1 <= N <= 10^18)

输出

共T行，如果A获胜输出A，如果B获胜输出B。

样例输入

3
3 5
3 4
1 9

样例输出

B
A
A

分析

此题与 51Nod_1072 威佐夫游戏 的不同之处在于N的取值范围从(1 <= N <= 2000000)变为(1 <= N <= 10^18)。如果继续使用(1+sqrt(5))/2来近似代替黄金分割数，就会产生精度问题，因此必须使用更加准确的值，也就是1.618033988749894848204586834。由于double类型的数据不能表示这么多位的小数，因此，可以把小数部分用整数进行存储，每9位表示一个整数，分别为p[0]=618033988,p[1]=749894848,p[2]=204586834。设两堆石子的数量差为temp，mod=1e9，则temp的范围为0<=temp<10^18，如果直接用temp乘p[i]会导致溢出，因此将temp分为高9位和低9位，分别用 l 和 r 表示,则有l=temp/mod, r=temp%mod 。此时temp*0.618033988749894848204586834可以表示为

(l*mod+r)*(p[0]*mod*mod+p[1]*mod+p[2])/(mod*mod)

=( r*p[2]+(l*p[2]+r*p[1])*mod+(l*p[1]+r*p[0])*mod^2+(l*p[0]) *mod^3)/(mod^2)

由于结果只取整数，因此每次只让运算结果的高9位参与下次运算

C++程序

#include<iostream>
#include<algorithm>
#include<cmath>

using namespace std;

typedef unsigned long long ull;

const ull p[3]={618033988,749894848,204586834};
const ull mod=1e9;

int main()
{
	int t;
	ull n,m;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%llu%llu",&n,&m);
		if(n>m) swap(n,m);
		ull temp=m-n;
		ull l=temp/mod;//temp的高9位 
		ull r=temp%mod;//temp的低9位 
		ull sum=r*p[2];
		sum=l*p[2]+r*p[1]+sum/mod;
		sum=l*p[1]+r*p[0]+sum/mod;
		sum=l*p[0]+sum/mod;
		sum+=temp;//加temp是由于sum现在的值为0.618....与temp的乘积，不是1.618...,所以还要加上一个temp
		if(sum==n)
		  printf("B\n");
		else
		  printf("A\n");
	}
	return 0;
}