Poj2960 S-Nim

Position：

• http://poj.org/problem?id=2960

---

List

Description

• 大意：有n堆石子，每堆石子个数已知，两人轮流从中取石子，每次可取的石子数x满足x属于集合S(k) = {s1,s2,s3…sk-1},问先拿者是否有必胜策略？
• 普通Nim取石子游戏但加了一些限制条件，比如每次只能取S={s1,s2,s3……}，就把前驱的条件改一下就行。

Knowledge

Sprague-Grundy Function-SG函数–博弈论
博弈论也是最近新学的知识，上面是一个写得很好的博客。
简单脑补：对于公平博弈（一般是NIM游戏），我们有一个重要的工具————就是SG函数。
SG函数的定义：
必败态的sg值为0，其余态的sg值为其后继状态的sg值的mex和。
其中mex和操作（mex{a1,a2,a3,…,ar}）的含义是a1,a2,a3,…,ar中最小的没出现过的自然数。
而对于组合游戏（就是由若干个子游戏组合而成，每个子游戏之间状态独立，每次操作任选一个子游戏操作），其sg值为所有子游戏sg值的异或和。如果一个状态求得sg值为0，则为必败态，否则为必胜态。（证明略，大致是因为先手总能通过一步使sg不为0的状态变为0，而sg为0的状态只能变成sg不为0的状态，最后不能操作的状态sg也为0）而一般sg都是打表找规律，常用分析方法：(1) 等差分析(2) 等比分析(3) 特定数值位置分析(4) 奇偶位置分析。对于多组数据都不同就要暴力求，如本题。

Solution

分析：
1.可将问题转化为n个子问题，每个子问题分别为：
从一堆x颗石子中取石子，每次可取的石子数为集合S(k)中的一个数
2.分析(1)中的每个子问题，易得:SG(x)=mex(SG[(x-s[i]>0)])(k>=i>=1)
3.后面就是SG函数的应用，根据Sprague-Grundy Therem：g(G)=g(G1)^g(G2)^g(G3)^…^g(Gn)即游戏的和的SG函数值是它的所有子游戏的SG函数值的异或，即SG(G) = SG(x1)^SG(x2)^…^SG(xn),故若SG(G)=0那么必输。

Notice

memset:①比for快②#include - cstring
map复杂度加一个log，对于加入的数少的情况用，else flag数组。

Code

// <S-Nim.cpp> - 08/03/16 20:18:06// This file is made by YJinpeng，created by XuYike's black technology automatically.// Copyright (C) 2016 ChangJun High School, Inc.// I don't know what this program is.#include <iostream>#include <vector>#include <algorithm>#include <cstring>#include <cstdio>#include <map>#include <cstdlib>#include <cmath>#define MOD 1000000007#define INF 1e9#define EPS 1e-10using namespace std;typedef long long LL;const int MAXN=10010;const int MAXM=100010;inline int max(int &x,int &y) {return x>y?x:y;}inline int min(int &x,int &y) {return x<y?x:y;}inline int getint() {    register int w=0,q=0;register char ch=getchar();    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();    if(ch=='-')q=1,ch=getchar();    while(ch>='0'&&ch<='9')w=w*10+ch-'0',ch=getchar();    return q?-w:w;}int n,T,sg[MAXN],m,ans,s[MAXN];inline int mex(int x){    //map<int,bool>f;    bool f[MAXN];    memset(f,0,sizeof(f));//fast <cstring>    int temp;    for(int i=1;i<=n;i++){        temp=x-s[i];        if(temp>=0){            if(sg[temp]==-1)sg[temp]=mex(temp);            f[sg[temp]]=1;        }    }    for(int i=0;;i++)if(!f[i])return i;}inline int SG(int x){    if(sg[x]==-1)sg[x]=mex(x);    return sg[x];}int main(){    freopen("S-Nim.in","r",stdin);    freopen("S-Nim.out","w",stdout);    while(n=getint(),n){        for(int i=1;i<=n;i++)s[i]=getint();        sg[0]=0;        for(int i=1;i<MAXN;i++)sg[i]=-1;        T=getint();        while(T--){            m=getint();ans=0;            while(m--)ans^=SG(getint());            if(ans)printf("W");else printf("L");        }        printf("\n");    }    return 0;}