HDOJ2177 取(2堆)石子游戏

有两堆各若干个物品，两个人轮流从某一堆或同时从两堆中取同样多的物品，规定每次至少取一个，多者不限，最后取光者得胜。相关的经典例题：

http://acm.hdu.edu.cn/showproblem.php?pid=2177这种情况下是颇为复杂的。我们用（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）。

假设面对的局势是（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，若都不是，那么就不是奇异局势。然后再按照上述法则进行，一定会遇到奇异局势。

所以只要if(m == (int)(t*(1 + d) / 2))  printf("0/n"); 就OK了。

 

#include <iostream> #include<cmath> using namespace std; int main() { int i,j,k,n,m,t,a,b,flag,tag; int a1,b1,a2,b2; double d = sqrt(5.0); while (scanf("%d %d",&a,&b)!=EOF) { a1 = b1 = a2 = b2 = -1; flag = 0; tag = 0; if (a==0 && b==0) { break; } if (a == (int)(abs(b-a)*(1 + d) / 2)) { printf("0/n"); continue; } for (i=1;i<=b;i++) { if ( (a-i>=0) && (flag==0))// && ((a-i)!=(b-i)))//两个同取 { // printf(" %d %d --> %d ------ %d /n/n",a-i,b-i,(a-i)^(b-i),a^b); if ((a-i) == (int)(abs(b-a)*(1 + d) / 2)) { a1 = a-i; b1 = b-i; flag = 1; } } if ( (a-i>=0) && (tag==0) )//&& (a-i)!=b) { if((a-i) == (int)(abs(b-a+i)*(1 + d) / 2)) { a2 = a-i; b2 = b; tag = 1; } } if ( (tag==0))// && a!=(b-i)) { if(((b-i)<a?(b-i):a) == (int)(abs((b-i)-a)*(1 + d) / 2)) { // printf("a=%d b-i=%d %d/n",a,(b-i)<a?(b-i):a, (int)(abs((b-i)-a)*(1 + d) / 2 +0.5)) ; a2 = a; b2 = b-i; tag = 1; } } if (tag && flag) { break; } } if (flag || tag) { printf("1/n"); if (flag) { printf("%d %d/n",a1,b1); } if (tag) { if (a2 > b2) { printf("%d %d/n",b2,a2); } else printf("%d %d/n",a2,b2); } } } return 0; }