[SDOI2019]移动金币(博弈论+阶梯Nim+按位DP)

首先可以把问题转化一下:m堆石子,一共石子数不超过(n-m)颗,每次可以将一堆中一些石子推向前一堆,无法操作则失败,问有多少种方法使得先手必胜?

然后这个显然是个阶梯Nim,然后有这样的结论:奇数层异或和为0。具体证明:参考这篇博客,当然不是我写的。如果不知道结论,里面有例题POJ1704可以做一下。

然后直接DP显然会T飞,考虑一个按位DP的技巧,f[i][j]表示确定前i位异或起来为0,剩下j个棋子的方案数。组合数相乘转移,注意一些细节即可。

[SDOI2019]移动金币(博弈论+阶梯Nim+按位DP)
#include<bits/stdc++.h>
using namespace std;
const int N=151000,mod=1e9+9;
int n,m,ans,f[19][N],fac[N],inv[N];
int C(int a,int b){return 1ll*fac[a]*inv[b]%mod*inv[a-b]%mod;}
int main()
{
    cin>>n>>m;
    fac[0]=inv[0]=inv[1]=1;for(int i=2;i<=n+m;i++)inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;
    for(int i=1;i<=n+m;i++)fac[i]=1ll*fac[i-1]*i%mod,inv[i]=1ll*inv[i-1]*inv[i]%mod;
    ans=C(n,m),n-=m;
    f[18][n]=1;
    for(int i=17;~i;i--)
    for(int j=0;j<=n;j++)
    for(int k=0;j+(2*k<<i)<=n&&k<=(m+1)/4;k++)
    f[i][j]=(f[i][j]+1ll*f[i+1][j+(2*k<<i)]*C((m+1)/2,2*k))%mod;
    for(int i=0;i<=n;i++)ans=(ans-1ll*f[0][i]*C(i+m/2,m/2)%mod+mod)%mod;
    cout<<ans;
}
View Code

 

[SDOI2019]移动金币(博弈论+阶梯Nim+按位DP)

上一篇:Uni-app实战项目注意事项


下一篇:mui真机调试时无法查找到手机