博弈论---Nimk问题。 dp再搞搞。
很容易看出,该游戏的终态是每两个棋子都紧靠着。当一颗棋子移动,另一方与该棋子对应的那一刻可以立即追上,使得仍旧紧靠,最终棋子动弹不得,游戏结束。
还能看出,对于白色棋子(先手),往左走没有意义。因为黑子(后手)可以紧随其上使得两者距离不变。同理黑子只往左走。(黄学长貌似提出了反例?)
所以,问题可以抽象为Nim,与传统Nim只能选1堆不同,你可以选1-d堆。
这个拓展问题叫做Nimk问题。对于这种问题,我们可以证明:当将n堆棋子化为二进制,每一位上如果1的个数mod(k+1)==0 则为必败态。
详细证明,大传送术!http://blog.csdn.net/weixinding/article/details/7321139
最后只需要计算方案数。使用dp,dp[i][j]表示当前在二进制第i位上,总计用了j石头的方案。转移方程为:
dp[i+1][j+a*(k+1)*bin[i]]+=dp[i][j]*C(n,a*(k+1));
注意组合数处理,取%等细节即可。
1 #include<bits/stdc++.h>
2 using namespace std;
3 #define mo 1000000007
4 #define N 10005
5 #define LL long long
6 LL c[N][205],dp[20][N],ans;
7 int n,k,d,K,bin[20];
8 void pre(){
9 bin[0]=1; for(int i=1;i<=15;i++) bin[i]=bin[i-1]<<1;
10 for(int i=0;i<=n;i++) c[i][0]=1;
11 for(int i=1;i<=n;i++)
12 for(int j=1;j<=min(i,K);j++)
13 c[i][j]=(c[i-1][j]+c[i-1][j-1])%mo;
14 }
15 LL C(int n,int m){
16 if(n-m<m) m=n-m;
17 return c[n][m];
18 }
19 LL cal(LL& x,LL y){
20 x=(x+y)%mo;
21 }
22 int main(){
23 scanf("%d%d%d",&n,&K,&d); pre(); dp[0][0]=1;
24 for(int i=0;i<15;i++)
25 for(int j=0;j<=n-K;j++)
26 for(int k=0;k*(d+1)<=K/2 && j+k*(d+1)*bin[i]<=n-K;k++)
27 cal(dp[i+1][j+k*(d+1)*bin[i]],dp[i][j]*C(K/2,k*(d+1)));
28 for(int i=0;i<=n-K;i++)
29 cal(ans,dp[15][i]*C(n-i-K/2,K/2));
30 LL tot=C(n,K);
31 cout<<(tot-ans+mo)%mo;
32 return 0;
33 }