Position:
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-10 using 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(" "); } return 0; }