bzoj2281&&jzoj2308[Sdoi2011]黑白棋 博弈论+动态规划

题意略。

明显可以发现如果黑白都是相邻的那么肯定先手必输，根据这个策略就可以做到30分。
还有40分水法但是我不会，很简单的一个式子，好像只有在一定情况下才能符合。
能考场上或者自己切满分算法的我只能跪下来膜拜。
既然已经知道必败态，那么我们可以假设第i个白色棋子和第i个黑色棋子互相匹配，那么一连串下来其实就可以把模型化为：把他们中间的格子看成石子堆，变成一个每次可以取d堆的nim游戏

那么问题就变为了，有n堆石子，每次可以从1~k堆中取任意数量的石子（不能不取），什么情况是必胜局面
有个结论就是d堆的石子数xor起来%(d+1)==0就是必败的，证明看
http://blog.csdn.net/weixinding/article/details/7321139
然后设f[i][j]表示前i位，放置了j个石头的必败方案数，
那么有f[i+1][j]=f[i][j+x*(d+1)2^i]，注意有C(K/2,x(d+1))种组合方式。
x枚举一下。
然后用总答案减去必败的就可以辣。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=1e5+5;
const int mo=1e9+7;
typedef long long ll;
int fac[N],n,m,d;
ll tot,ans;
ll mi[25];
ll c[10005][205];
ll f[25][10005];
inline void add(ll &x,ll y){x=(x+y)%mo;}
inline int pow(int a,int b)
{
    int r=1;
    while (b)
    {
        if (b&1)r=1ll*r*a%mo;
        a=1ll*a*a%mo;
        b>>=1;
    }
    return r;
}
inline void pre()
{
    fo(i,0,n)c[i][0]=1;
    fo(i,1,n)
    fo(j,1,min(2*m,i))
    c[i][j]=(c[i-1][j]+c[i-1][j-1])%mo;
}
inline int C(int n,int m)
{
    if (m>n-m)m=n-m;
    return c[n][m];
}
int main()
{
    mi[0]=1;
    fo(i,1,16)mi[i]=mi[i-1]*2%mo;
    scanf("%d%d%d",&n,&m,&d);
    m/=2;
    pre();
    f[0][0]=1;
    fo(i,0,14)
    fo(j,0,n-2*m)
    for(int p=0;p*(d+1)<=m&&j+p*(d+1)*mi[i]<=n-2*m;p++)
    {
        add(f[i+1][j+p*(d+1)*mi[i]],f[i][j]*C(m,p*(d+1)));
    }
    fo(i,0,n-2*m)
    add(ans,f[15][i]*C(n-i-m,m));
    tot=C(n,m*2);
    printf("%lld\n",(tot-ans+mo)%mo);
}