[BZOJ 3576][HNOI 2014]江南乐(SG博弈)

题目链接

http://www.lydsy.com/JudgeOnline/problem.php?id=3576

思路

显然是个SG博弈。
以下内容转自http://blog.csdn.net/gromah/article/details/27326991：

对于10分的数据：明显当石子数大于1时，先手胜。因为只有一堆，先手只要把石子全拆成只有1个的若干堆就赢了，但如果本来就只有1个石子，就是后手赢，复杂度O(1)。
对于70分的数据：直接暴力求在当前的 f
值下每个石堆对的的sg值即可。每堆石子有O(n)中分裂方式，每个分裂方式的sg求值是O(1)的，有O(n)个石子，复杂度O(n ^ 2)。
对于100分的数据：我们考虑两个性质

1、（n / i）（商向下取整），至多有 2 * sqrt(n) 个不同的结果。 2、如果 n / i == n / (i + 2)
，那么通过将石子分裂成 i 份得到的 i 堆石子的sg值异或和与将石子分裂成 i + 2 份石子所得到的结果相同。

如果上述两个性质是对的，我们就可以只枚举 n / i 来求解，因为有很多分裂方式在本质上是相同的，我们完全可以不考虑。
现在我们来证明这些性质：

证明性质1：在1 - sqrt(n) 之间，因为只有 sqrt(n) 个数，所以 (n / i) 在1 - sqrt(n)
之间当然最多只有sqrt(n) 个结果。
在 sqrt(n) - n 之间，因为 n / n = 1 , n / sqrt(n) = sqrt(n)，n 除以该范围内的数的结果必然在1 - sqrt(n) 之间 。所以 (n / i) 在 sqrt(n)
- n 之间也只会出现sqrt(n) 个不同的结果。

证明性质2：因为把 n 个石子按照题目要求分成 i 份，会产生 n % i 个石子数为 n / i + 1 的石头堆，i - n % i
个石子数为 n / i 的石头堆。（这个挺显然的）
而把 n 个石子按照题目要求分成 i + 1 份，就相当于在分成了 i 份的基础上，从 n / i 个石子数为 n / i + 1 的石头堆中，各拿出1个石子，组成一个新的石子数为n / i 的石头堆。
而把 n 个石子按题目要求分成 i + 2 份的话，就相当于做了两次组新石子堆的工作，其石子堆的奇偶性不变，因为你操作了两次，所以各种石子堆的改变量就必然是2的倍数，既然这样，石子堆的sg值异或和就不会变。
所以性质 2 是正确的。
好的，现在我们保证了这个算法的正确性，来说说具体做法吧。 对于每个石子数n，我们枚举它的 n / i 值，然后取分裂成 n / i 份和 n / (i +1) 份的两种情况就可以了。 实现过程中，我们还可以进行记忆化搜索来提高程序的效率。 每个石子数的sg求值是
O(sqrt(n)) 的，有 O(n) 个石子，所以时间复杂度是 O(n * sqrt(n)) 的。

嗯，上面的题解讲得非常清楚，我就不用多说了。。。对于每堆石子，也就是每个游戏而言，它的每个后继都是若干个子游戏的集合，是一堆SG值的亦或和，那么它的SG值就是mex{每个后继的SG值}
70分的暴力还是比较好写的，然后就是各种优化了
在枚举分成的堆数i时，只需要枚举大概sqrt(n)规模的i，这个在题目中有证明，然后我们还需要关注i的奇偶性，因此是要同时求i和i+1的SG值的
坑爹的是这个题目的比较方式是不过滤掉行末最后一个空格和文末最后一个换行符的，导致我贡献3发PE，教训啊

代码

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <algorithm>

#define MAXN 101000

using namespace std;

int T,F;
int mex[MAXN],sg[MAXN];
bool visit[MAXN];
int id; //mex数组中为了防止冲突，给每次DFS留的标记

int GetSG(int n) //求一个大小为n的堆的SG值
{
    if(n<F) return 0; //SG=0:堆的大小小于F，无法操作
    if(visit[n]) return sg[n]; //已经求过sg(n)，就不重复劳动了
    visit[n]=true;
    for(int i=2;i<=n;i=n/(n/i)+1) //枚举分成i堆
        for(int j=i;j<=i+1&&j<=n;j++) //枚举分成j堆,j=i或i+1
        {
            GetSG(n/j); //一种是分成的堆大小为n/j
            GetSG(n/j+1); //另一种是分成的堆大小为n/j+1
        }
    id++;
    for(int i=2;i<=n;i=n/(n/i)+1)
        for(int j=i;j<=i+1&&j<=n;j++)
        {
            int tmp=0;
            if((n%j)%2==1) tmp^=sg[n/j+1]; //奇数个大小为(n/j+1)的堆就需要亦或它的SG值，偶数个就相当于没亦或
            if((j-n%j)%2==1) tmp^=sg[n/j]; //奇数个大小为(n/j)的堆就需要亦或它的SG值，偶数个就相当于没亦或
            mex[tmp]=id;
        }
    for(sg[n]=0;mex[sg[n]]==id;sg[n]++);
    return sg[n];
}

int main()
{
    scanf("%d%d",&T,&F);
    while(T--)
    {
        int n,xorsum=0; //n堆石子,xorsum是SG值的亦或和
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
            int x;
            scanf("%d",&x);
            xorsum^=GetSG(x);
        }
        if(xorsum==0) //先手必败
            printf("0%c",T?' ':'\n');
        else
            printf("1%c",T?' ':'\n');
    }
    //printf("\n");
    return 0;
}