[BZOJ2281][SDOI2011]黑白棋(K-Nim博弈)

2281: [Sdoi2011]黑白棋

Time Limit: 3 Sec  Memory Limit: 512 MB
Submit: 626  Solved: 390
[Submit][Status][Discuss]

Description

小A和小B又想到了一个新的游戏。

这个游戏是在一个1*n的棋盘上进行的，棋盘上有k个棋子，一半是黑色，一半是白色。

最左边是白色棋子，最右边是黑色棋子，相邻的棋子颜色不同。

小A可以移动白色棋子，小B可以移动黑色的棋子，他们每次操作可以移动1到d个棋子。

每当移动某一个棋子时，这个棋子不能跨越两边的棋子，当然也不可以出界。当谁不可以操作时，谁就失败了。

小A和小B轮流操作，现在小A先移动，有多少种初始棋子的布局会使他胜利呢？

Input

共一行，三个数，n,k,d。

Output

输出小A胜利的方案总数。答案对1000000007取模。

Sample Input

10 4 2

Sample Output

182

HINT

1<=d<=k<=n<=10000, k为偶数，k<=100。

Source

stage 2 day1

很有意思的一道博弈题，可惜HZWER学长给出了反例。

那么这一题通过手玩可以发现，最终状态必定是所有棋子全部扎堆在棋盘左端或右端，棋子之间没有间隙。不过仔细观察可以发现，可能在游戏状态中会出现所有棋子扎堆但不在棋盘一端的情况，其实那个时候就已经决定了最终的胜负。因为只要一方朝自己来的方向走了，则另一方必定能也往那边走一步，最终会步步紧逼直到走到棋盘一端。

根据这一点，感性理解一下，这个游戏就是一个把对方棋子“怼”过去的过程，谁怼赢了就是胜者。所以从一开始双方都一定拼尽全力往对面怼，所以有一个结论：先手不可能往左走，后手不可能往右走。

这样这个问题就变成了一个取石子游戏，每对相邻的白子和黑子之间的格子数是石子数（显然共有K/2堆石子），每人每次选不超过k堆取一个石子。

这个问题叫K-Nim，结论是:将所有石子数转成二进制，如果对于每一位二进制，这一位上为1的石子堆数都能被k+1整除则为必败态，否则为必胜态。

证明主要思路是: 1.最终态二进制每一位都为0必为必败态。2.只要有某位的1的个数不被k+1整除，则必然有一种走法使每一位都被整除。 3.如果每一位都被k+1整除，则无论怎么走都不可能使得每一位都仍然能被整除。

这三点分别保证了：最终态是必败态。必胜态必定能走到必败态。必败态只能走到必胜态。

详细证明:http://blog.csdn.net/weixinding/article/details/7321139

这样，我们分别用了“寻找最终态”和“模仿”的技巧将问题转化为了K-Nim问题。回到这一题，最终答案=总方案数-必败态的方案数。

设$f_{i,j}$表示前$i$个二进制位共放了$j$个石子的方案数，则$$ans=C_n^K-\sum_{i=0}^{n-K} f_{s,i}*C_{n-i-K/2}^{K/2}$$s为最高位的1，这里取15就够了。

考虑$f$的转移方程即可：$$f_{i+1,j+k*(d+1)*(1<<i)}\ \ +=\ \ f_{i,j}*C_{K/2}^{k*(d+1)}$$

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #define rep(i,l,r) for (int i=l; i<=r; i++)
 typedef long long ll;
 using namespace std;
 
 const ll N=,mod=;
 ll tot,ans,bin[],c[N][],f[][N];
 int n,K,d,p;
 
 void add(ll &x,ll y){ x=(x+y)%mod; }
 void pre(){
     rep(i,,n) c[i][]=;
     rep(i,,n) rep(j,,min(*K,i)) c[i][j]=(c[i-][j]+c[i-][j-])%mod;
 }
 int C(int n,int m){ if (m>n-m) m=n-m; return c[n][m]; }
 
 int main(){
     freopen("bzoj2281.in","r",stdin);
     freopen("bzoj2281.out","w",stdout);
     bin[]=; rep(i,,) bin[i]=bin[i-]<<;
     scanf("%d%d%d",&n,&K,&d); K>>=;
     pre(); f[][]=;
     rep(i,,) rep(j,,n-*K)
         for (int k=; k*(d+)<=K && j+k*(d+)*bin[i]<=n-*K; k++)
             add(f[i+][j+k*(d+)*bin[i]],f[i][j]*C(K,k*(d+)));
     rep(i,,n-*K) add(ans,f[][i]*C(n-i-K,K));
     tot=C(n,K*); printf("%lld\n",(tot-ans+mod)%mod);
     return ;
 }