[CC]Bored of Nim

题意简述

给出一个$n$个元素的数列，从中拿出连续$m$个，做Nim游戏，问每一种长度的$m$有多少种情况先手必输。

数据范围

$1 \leq n \leq 10^5$
$1 \leq a_i \leq 10^9$

思路

做异或前缀和，区间xor为0就是前缀和有两种相同的数。
对数的个数进行讨论，如果$\leq B$那就暴力。
如果$>B$那就做FFT。
$j,i(j>i)$贡献到$j-i$的位置上。
这样不符合卷积的形式。我们把$A_i$和$B_{n-j}$进行卷积。
$B=\sqrt {n\log n}$时取得最优复杂度。
因为FFT的大常数，我们适当扩大 B <script type="math/tex" id="MathJax-Element-502">B</script>。

代码

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<vector>
#include<cstring>
using namespace std;
#define MAXN 100010
#define MAXL 262150
const double pi=acos(-1);
pair<int,int> DS[MAXN];
struct C{
    double x,y;
    C (double _x=0,double _y=0)
    {
        x=_x,y=_y;
    }
    C operator * (const C &n1)
    {
        return C(x*n1.x-y*n1.y,x*n1.y+y*n1.x);
    }
    C operator *= (const C &n1)
    {
        double a=x,b=y;
        x=a*n1.x-b*n1.y;
        y=b*n1.x+a*n1.y;
    }
    C operator /= (const int &n1)
    {
        x/=n1,y/=n1;
    }
    C operator + (const C &n1)
    {
        return C(x+n1.x,y+n1.y);
    }
    C operator - (const C &n1)
    {
        return C(x-n1.x,y-n1.y);
    }
    void operator += (const int &n1)
    {
        x+=n1;
    }
    double real()
    {
        return x;
    }
}A[MAXL],B[MAXL];
int len,ti,tmp,block;
int r[MAXL];
int T,n,now;
int seq[MAXN];
long long ans[MAXN];
vector<int> vec[MAXN];
void fft(C *a,int f)
{
    for (int i=0;i<len;i++)
        if (i<r[i])
            swap(a[i],a[r[i]]);
    for (int i=1;i<len;i<<=1)
    {
        C wn(cos(pi/i),f*sin(pi/i));
        for (int j=0;j<len;j+=(i<<1))
        {
            C w=1;
            for (int k=0;k<i;k++,w*=wn)
            {
                C x=a[j+k],y=w*a[j+k+i];
                a[j+k]=x+y;a[j+k+i]=x-y;
            }
        }
    }
}
int main()
{
    scanf("%d",&T);
    while (T--)
    {
        scanf("%d",&n);
        block=10000;
        for (ti=0,len=1,tmp=n<<1;len<=tmp;ti++,len<<=1);
        for (int i=0;i<len;i++)
            r[i]=(r[i>>1]>>1)|((i&1)<<ti-1);
        for (int i=1;i<=n;i++)
            scanf("%d",&seq[i]);
        for (int i=1;i<=n;i++)
            seq[i]^=seq[i-1];
        for (int i=0;i<=n;i++)
            DS[i+1]=make_pair(seq[i],i);
        sort(DS+1,DS+n+2);
        now=0;
        DS[0].first=-1;
        for (int i=1;i<=n+1;i++)
        {
            if (DS[i].first!=DS[i-1].first)
                vec[++now].clear();
            vec[now].push_back(DS[i].second);
        }
        memset(ans,0,sizeof(ans));
        for (int i=1;i<=now;i++)
        {
            tmp=vec[i].size();
            if (tmp<=block)
                for (int a1=0;a1<tmp;a1++)
                    for (int a2=a1+1;a2<tmp;a2++)
                        ans[vec[i][a2]-vec[i][a1]]++;
            else
            {
                memset(A,0,sizeof(A));
                memset(B,0,sizeof(B));
                for (int a1=0;a1<tmp;a1++)
                    A[vec[i][a1]]+=1,B[n-vec[i][a1]]+=1;
                fft(A,1);fft(B,1);
                for (int j=0;j<len;j++)
                    A[j]*=B[j];
                fft(A,-1);
                for (int j=0;j<len;j++)
                    A[j]/=len;
                for (int j=n+1;j<=2*n;j++)
                    ans[j-n]+=(long long)(A[j].real()+0.5);
            }
        }
        for (int i=1;i<=n;i++)
            printf("%lld%c",ans[i]," \n"[i==n]);
    }
    return 0;
}