HOJ 2645 WNim（sg函数找规律以及求解）

最新推荐文章于 2020-09-07 23:59:51 发布

沙雕.

最新推荐文章于 2020-09-07 23:59:51 发布

阅读量225

点赞数

分类专栏：博弈论

本文链接：https://blog.csdn.net/weixin_43768644/article/details/89280332

版权

博弈论专栏收录该内容

7 篇文章 0 订阅

订阅专栏

题目链接：http://acm.hit.edu.cn/problemset/2645

参考了这篇题解：https://blog.csdn.net/y1196645376/article/details/52165245

解题思路：

符合nim游戏的博弈，重点就是找sg函数的规律以及对应求法。

找规律请参考那位大佬的博客，我仅描述一下怎样求sg值

举个栗子：

假设有10个数，那么我们就有1个10位的二进制数，某一位为1表示那个数参与此次求多个数的最小公倍数

（多个数最小公倍数求法：ans=a1*a2/gcd（a1，a2），之后每次ans = ans*a[i]/gcd（ans，a[i]））

所以外循环 from 1 to 2^10-1

内循环判断 1-10位上是否为1，（判断方法 i & (1<<j-1)）///j表示第几位，可以直接将内层循环写为0-9，更方便

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define debug(x) printf("Line:%s----\n",#x)
#define ll long long
using namespace std;

int think[11];

ll gcd(ll a,ll b){return b==0?a:gcd(b,a%b);}

int getsg(int x,int tot)
{
    int ans=x;
    for (int i=1;i<(1<<tot);i++){//debug(1);
        int cnt=0;
        ll lcm=1;
        for (int j=0;j<tot;j++){
            if (i & (1<<j)) lcm = lcm*think[j]/gcd(lcm,think[j]),cnt++;
            if (lcm>x) break;
        }//debug(2);
        if (cnt%2) ans -= x/lcm;
        else ans += x/lcm;//debug(3);
    }//debug(4);
    return ans;
}

int main()
{
    int t,n,m;
    scanf("%d",&t);
    while (t--){
        scanf("%d %d",&n,&m);
        for (int i=0;i<n;i++) scanf("%d",think+i);
        int x,flag=0,ans=0;
        for (int i=0;i<m;i++){
            scanf("%d",&x);
            if (flag) continue;
            for (int j=0;j<n;j++) if (x%think[j]==0) flag = 1;
            //printf("sg[%d]=%d\n",x,getsg(x,n));
            ans ^= getsg(x,n);
        }
        //printf("ans=%d\n",ans);
        if (flag) printf("Xiao Hong\n");
        else printf("%s\n",ans ? "Xiao Hong":"Xiao Gang");
    }
    return 0;
}