bzoj3576[HNOI2014]江南乐

首先整个局面的SG值等于各个独立子局面的SG值异或和，于是只要求SG(1~100000)了。考虑将一堆i个石头分成j堆时，它的后继状态是 {\((j-i \ mod\ j)\) 个 \((i/j)\)，\((i\ mod \ j)\) 个 \((i/j+1)\)}，这个后继状态的SG值同样等于这j个局面的SG异或和，于是要计算SG(i)就只需枚举将它分成j=1~i堆(j>i是没有意义了的，在j=i时已经考虑到了)，然后把后继状态的SG取mex即可。
于是就有\(O(n^2)\)的算法:

for(int i=1;i<F;++i) SG[i]=0;
for(int i=F;i<=100000;++i)
{
    for(int j=2;j<=i;++j)
    {
        int num=(j-i%j)&1,num1=(i%j)&1;
        int sg=(num*SG[i/j])^(num1*SG[i/j+1]);
        vis[sg]=1; stk[++top]=sg;
    }
    for(int j=0;;++j) if(!vis[j]) {SG[i]=j;break;}
    while(top) vis[stk[top]]=0,top--;
}
for(int cas=1;cas<=T;++cas)
{
    n=read(); int ans=0;
    for(int i=1;i<=n;++i) a[i]=read(),ans^=SG[a[i]];
    printf("%d ",!(!ans));
}

然后容易想到这个可以分块优化，对于一段\(i/j\)相同的j，\(j-i\ mod \ j\) 和 \(i\ mod \ j\)有一个是公差为偶数的等差数列，一个是公差为奇数的等差数列，所以只需要考虑一段\(i/j\)相同的j的其中两个即可。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<vector>
#define P puts("lala")
#define cp cerr<<"lala"<<endl
#define fi first
#define se second
#define pb push_back
#define ln putchar('\n')
using namespace std;
inline int read()
{
    char ch=getchar();int g=1,re=0;
    while(ch<'0'||ch>'9') {if(ch=='-')g=-1;ch=getchar();}
    while(ch<='9'&&ch>='0') re=(re<<1)+(re<<3)+(ch^48),ch=getchar();
    return re*g;
}
typedef long long ll;
typedef pair<int,int> pii;

const int N=125;
const int M=100050;
int a[N],F,n;
int SG[M],stk[M],top=0;
bool vis[M];

int main()
{
#ifndef ONLINE_JUDGE
    freopen("1.in","r",stdin);freopen("1.out","w",stdout);
#endif
    int T=read(); F=read();
    for(int i=1;i<F;++i) SG[i]=0;
    for(int i=F;i<=100000;++i)
    {
        for(int l=2,r;l<=i;l=r+1)
        {
            r=i/(i/l);
            int num=0,num1=0,sg=0;
            num=(l-i%l)&1; num1=(i%l)&1;
            sg=(num*SG[i/l])^(num1*SG[i/l+1]);
            vis[sg]=1;
            stk[++top]=sg;

            if(l!=r)
            {
                l++;
                num=(l-i%l)&1; num1=(i%l)&1;
                sg=(num*SG[i/l])^(num1*SG[i/l+1]);
                vis[sg]=1;
                stk[++top]=sg;
            }
        }
        for(int j=0;;++j) if(!vis[j]) {SG[i]=j;break;}
        while(top) vis[stk[top]]=0,top--;
    }
    for(int cas=1;cas<=T;++cas)
    {
        n=read();
        int ans=0;
        for(int i=1;i<=n;++i) a[i]=read(),ans^=SG[a[i]];
        printf("%d ",!(!ans));
    }
    return 0;
}

转载于:https://www.cnblogs.com/thkkk/p/8677671.html