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;
}