[51nod1169] 石子游戏

题目大意

给定n堆石子的个数。两人轮流操作，每个人每次操作可以从一堆里取出至少一个石子。不能拿的人输。
现在可以在游戏开始前，选择任意堆石子（不能全选），并从中分别取出任意颗石子。求有多少种取法使先手必败。

n≤100 每堆石子个数不超过$10^9$，答案对$10^9+7$取模

分析

首先一堆石子的sg值为石子个数，然后选手必败需要满足所有sg值异或和为0。
然后考虑到不能全选这一条件，可以考虑容斥：首先算出随便取的答案，然后给每堆石子数量都减一，然后算出此时的答案，相减即可得到正确答案。
剩下的问题就是求任意取的答案了。

首先从高位开始枚举二进制位，假设当前枚举到k，a[]是原数组，b[]是取完后的数组，枚举k表示所有b[i]比k高的位都等于a[i]，然后存在i，b[i]与a[i]的第k位不同。低位则不管。

接着枚举一个i，表示a[i]与b[i]的第k位不同（为了满足题目条件，可以发现a[i]的第k位为1，b[i]的第k位为0）。然后就要保证二进制第k位的异或和为0。考虑低位，由于b[i]的低位可以任意取，所以不用考虑其它数的低位异或和是多少了。
为了避免算重，要保证i是第一个不同的数。

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

const int N=105,M=30,mo=1e9+7;

typedef long long LL;

int n,a[N],ans,g[N][2],b[N];

int calc(int *a)
{
    ans=0;
    int k;
    for (k=M;k>=0;k--)
    {
        memset(g,0,sizeof(g));
        g[n+1][0]=1;
        int p=1<<k;
        for (int i=n;i;i--)
        {
            b[i]=((a[i]&p)>0);
            if (!b[i])
            {
                g[i][0]=(LL)g[i+1][0]*(a[i]%p+1)%mo;
                g[i][1]=(LL)g[i+1][1]*(a[i]%p+1)%mo;
            }else
            {
                g[i][0]=((LL)g[i+1][1]*(a[i]%p+1)+(LL)g[i+1][0]*p)%mo;
                g[i][1]=((LL)g[i+1][0]*(a[i]%p+1)+(LL)g[i+1][1]*p)%mo;
            }
        }
        int j=1,t=0;
        for (int i=1;i<=n;i++)
        {
            if (b[i]) ans=(ans+(LL)j*g[i+1][t])%mo;
            j=(LL)j*(a[i]%p+1)%mo; t^=b[i];
        }
        if (t) break;
    }
    if (k<0) ans++;
    return ans;
}

int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++) scanf("%d",&a[i]);
    int s=calc(a);
    for (int i=1;i<=n;i++) a[i]--;
    s=(s+mo-calc(a))%mo;
    printf("%d\n",s);
    return 0;
}