SDOI 2019 移动金币 题解

题目传送门

题目大意： 一个长度为 $n$ 的序列上有 $m$ 个金币，两个人轮流操作，一个人可以将一个金币向左移动任意格，但是不能越过别的金币，问有多少种局面先手必胜。

题解

感觉像这两题（1，2）的合体，不过合的也很巧妙。

转化一下这个博弈：相当于有 $m+1$ 堆石子，每次可以将一堆石子中的若干个移到下一堆内，第 $m+1$ 堆的不可以被操作。

对于第 $i$ 堆石子，假如 $m+1-i$ 为偶数，那么这堆石子其实是没用的，假如 $\text{Alice}$ 将其中若干个移到下一堆，那么 $\text{Bob}$ 接着将这若干个再移到下一堆，那么这若干个依然在 $m+1-i$ 为偶数的堆内，一直搞下去那么这些石子最后就会停在第 $m+1$ 堆内，此时依然是 $\text{Alice}$ 操作。

那么如果将 $m+1-i$ 为奇数的堆内的石子移到下一堆，那么这些石子就变成了没用的石子了，则 $m+1-i$ 为奇数的堆就构成了一个 $Nim$ 游戏，于是 $\text{dp}$ 一下就做完了，$\text{dp}$ 部分和上面链接里第二题是一样的，做过的话就不用再看我讲一次了。

现在的子问题是 $n$ 个石子分成 $m$ 堆构成的 $Nim$ 游戏有多少种情况先手必胜（注意这里 $n,m$ 和上面的不是同一个），这不太好求，正难则反，考虑必败的情况数，令 $f[i]$ 表示 $i$ 个石子分 $m$ 堆先手必败的方案数。先手必败需要满足：将每堆石子数量转成二进制数后，对于任意 $j$，第 $j$ 位为 $1$ 的数的个数为偶数。考虑枚举每一位 $1$ 的个数，用组合数将他们分配到 $m$ 堆内，就可以做 $dp$ 了，时间复杂度 $O(17nm)$。

代码如下：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define maxn 150010
#define mod 1000000009

int n,m,k;
int Binom[maxn][61];
void add(int &x,int y){x=(x+y>=mod?x+y-mod:x+y);}
int add(int x){return x>=mod?x-mod:x;}
void BinomInit(){
	Binom[0][0]=1;
	for(int i=1;i<=n;i++){
		Binom[i][0]=1;
		for(int j=1;j<=i&&j<=m;j++){
			Binom[i][j]=add(Binom[i-1][j]+Binom[i-1][j-1]);
		}
	}
}
int f[maxn],tmp[maxn];

int main()
{
	scanf("%d %d",&n,&m);k=(m+1)/2;
	BinomInit();f[0]=1;
	for(int i=0;i<=17;i++){
		int bin=1<<i;
		for(int j=0;j<=k&&j*bin<=n-m;j+=2){
			for(int l=0;l+j*bin<=n-m;l++){
				add(tmp[l+j*bin],1ll*f[l]*Binom[k][j]%mod);
			}
		}
		for(int j=0;j<=n-m;j++)f[j]=tmp[j],tmp[j]=0;
	}
	int ans=0;
	for(int i=0;i<=n-m;i++)add(ans,1ll*f[i]*Binom[n-m-i+(m+1-k-1)][m+1-k-1]%mod);
	ans=(Binom[n-m+(m+1-1)][m+1-1]-ans+mod)%mod;
	printf("%d",ans);
}