【Luogu】P2490黑白棋（博弈DP）

　　题目链接

　　题解链接

　　

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<cctype>
#define mod 1000000007
#define maxn 10020
using namespace std;
inline long long read(){
    long long num=,f=;
    char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-')    f=-;
        ch=getchar();
    }
    while(isdigit(ch)){
        num=num*+ch-'';
        ch=getchar();
    }
    return num*f;
}

long long mul[maxn];
long long inv[maxn];
long long f[][maxn];

inline long long Pow(long long a,long long b){
    long long ret=;
    while(b){
        if(b&)    ret=ret*a%mod;
        a=a*a%mod;
        b>>=;
    }
    return ret;
}

inline long long C(long long n,long long m){
    return (mul[n]*inv[m]%mod)*inv[n-m]%mod;
}

int main(){
    int n=read(),e=read()>>,w=read();    mul[]=inv[]=;
    for(int i=;i<=n;++i){
        mul[i]=mul[i-]*i%mod;
        inv[i]=Pow(mul[i],mod-);
    }
    f[][]=;
    for(int i=;i<=;++i){
        for(int j=;j<=n-(e<<);++j){
            //printf("%d %d\n",i,j);
            for(int k=;k*(w+)<=e&&j+k*(w+)*(<<i)<=n-(e<<);++k){
                long long &now=f[i+][j+k*(w+)*(<<i)];
                now=(now+f[i][j]*C(e,k*(w+)))%mod;
            }
        }
    }
    long long ans=;
    for(int i=;i<=n-(e<<);++i)    ans=(ans+f[][i]*C(n-i-e,e))%mod;
    printf("%lld\n",(C(n,e<<)-ans+mod)%mod);
    return ;
}