[bzoj4347][DP]Nim z utrudnieniem

Description

A和B两个人玩游戏，一共有m颗石子，A把它们分成了n堆，每堆石子数分别为a[1],a[2],…,a[n]，每轮可以选择一堆石子，取掉任意颗石子，但不能不取。谁先不能操作，谁就输了。在游戏开始前，B可以扔掉若干堆石子，但是必须保证扔掉的堆数是d的倍数，且不能扔掉所有石子。A先手，请问B有多少种扔的方式，使得B能够获胜。

Input

第一行包含两个正整数n,d(1<=n<=500000,1<=d<=10)。
第二行包含n个正整数a[1],a[2],…,an。
本题中m不直接给出，但是保证m<=10000000。

Output

输出一行一个整数，即方案数对10^9+7取模的结果。

Sample Input

5 2

1 3 4 1 2

Sample Output

2

题解

劲
显然题意就是要你求有多少种取数的方案 让这些数的异或和为全部数的异或和
并且取出来的数的个数是D的倍数
大力DP
设$f[i][j][k]$表示第i个数 个数%D=j 异或和为K的方案数
我们把序列从小到大排序 这样每次枚举的最大异或值不会超过2*A[i]
发现j只会从j-1转移 可以不用滚动了…直接转
j=0的地方特判一下…
这样的复杂度其实只有$O(DM)$的

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<ctime>
#define mod 1000000007
#define LL long long
#define mp(x,y) make_pair(x,y)
using namespace std;
inline int read()
{
	int f=1,x=0;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
inline void write(int x)
{
	if(x<0)putchar('-'),x=-x;
	if(x>9)write(x/10);
	putchar(x%10+'0');
}
inline void print(int x){write(x);printf(" ");}
int a[510000],n,D,S;
int f[10][1048576],g[1048576];
void ad(int &x,int y){x+=y;if(x>mod)x-=mod;}
int main()
{
	n=read();D=read();
	for(int i=1;i<=n;i++)S^=(a[i]=read());
	sort(a+1,a+1+n);
	f[0][0]=1;
	for(int i=1;i<=n;i++)
	{
		int pos=min(a[i]<<1,1048575);
		for(int j=0;j<=pos;j++)g[j]=f[D-1][j];
		for(int k=D-1;k>=1;k--)
			for(int j=0;j<=pos;j++)ad(f[k][j],f[k-1][j^a[i]]);
		for(int j=0;j<=pos;j++)ad(f[0][j],g[j^a[i]]);
	}
	LL ans=f[0][S];
	if(!(n%D)){ans=(ans-1+mod)%mod;}
	printf("%d\n",ans);
	return 0;
}