题目大意: 一个长度为 n n n 的序列上有 m m m 个金币,两个人轮流操作,一个人可以将一个金币向左移动任意格,但是不能越过别的金币,问有多少种局面先手必胜。
题解
转化一下这个博弈:相当于有 m + 1 m+1 m+1 堆石子,每次可以将一堆石子中的若干个移到下一堆内,第 m + 1 m+1 m+1 堆的不可以被操作。
对于第 i i i 堆石子,假如 m + 1 − i m+1-i m+1−i 为偶数,那么这堆石子其实是没用的,假如 Alice \text{Alice} Alice 将其中若干个移到下一堆,那么 Bob \text{Bob} Bob 接着将这若干个再移到下一堆,那么这若干个依然在 m + 1 − i m+1-i m+1−i 为偶数的堆内,一直搞下去那么这些石子最后就会停在第 m + 1 m+1 m+1 堆内,此时依然是 Alice \text{Alice} Alice 操作。
那么如果将 m + 1 − i m+1-i m+1−i 为奇数的堆内的石子移到下一堆,那么这些石子就变成了没用的石子了,则 m + 1 − i m+1-i m+1−i 为奇数的堆就构成了一个 N i m Nim Nim 游戏,于是 dp \text{dp} dp 一下就做完了, dp \text{dp} dp 部分和上面链接里第二题是一样的,做过的话就不用再看我讲一次了。
现在的子问题是 n n n 个石子分成 m m m 堆构成的 N i m Nim Nim 游戏有多少种情况先手必胜(注意这里 n , m n,m n,m 和上面的不是同一个),这不太好求,正难则反,考虑必败的情况数,令 f [ i ] f[i] f[i] 表示 i i i 个石子分 m m m 堆先手必败的方案数。先手必败需要满足:将每堆石子数量转成二进制数后,对于任意 j j j,第 j j j 位为 1 1 1 的数的个数为偶数。考虑枚举每一位 1 1 1 的个数,用组合数将他们分配到 m m m 堆内,就可以做 d p dp dp 了,时间复杂度 O ( 17 n m ) O(17nm) 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);
}