• [SDOI2019] 移动金币


    分析

    阶梯NIM模型:共有m+1堆石子,石子总数不超过n-m,求必胜的,即奇数堆石子数目异或和非零的局面数。补集转化,答案C(n,m)-奇数堆石子数目异或和位0的局面数。

    可以想到按位dp,设f[i,j]表示已经考虑了前i位(异或和0),石子和为j的方案数;转移时考虑下一位出现的被统计1的个数k,k为偶数,带系数C((m+1)/2,k),即m+1/2个奇数堆中出现k个1的方案。

    最后将没有分配的n-m-i个石子放入m/2个偶数堆中,可以一个都不放。

    实现

    #include <bits/stdc++.h>
    #define ll long long 
    using namespace std;
    
    const int N=2e5+10;
    const int mod=1e9+9;
    
    int n,m;
    ll fc[N],fv[N],f[20][N];
    
    inline ll c(int n,int m) {
    	return fc[n]*fv[m]%mod*fv[n-m]%mod;
    }
    
    int main() {
    	scanf("%d%d",&n,&m);
    	fc[0]=fc[1]=fv[0]=fv[1]=1;
    	for(int i=2; i<=n+m; ++i) fv[i]=fv[mod%i]*(mod-mod/i)%mod;
    	for(int i=2; i<=n+m; ++i) fv[i]=fv[i-1]*fv[i]%mod,fc[i]=fc[i-1]*i%mod;
    	f[19][0]=1;
    	for(int i=19; i; --i) 
    	for(int j=0; j<=n-m; ++j) if(f[i][j]) 
    		for(int k=0; k<=(m+1)/2&&j+(k<<(i-1))<=n-m; k+=2) 
    			(f[i-1][j+(k<<(i-1))]+=f[i][j]*c((m+1)/2,k)%mod)%=mod;		
    	int ans=c(n,m);
    	for(int i=0; i<=n-m; ++i) ans=(ans-f[0][i]*c(n-m-i+m/2,m/2)+mod)%mod;
    	printf("%d
    ",ans);
    	return 0;
    }
    

    暴力

    #include <bits/stdc++.h>
    #include <typeinfo>
    #define ll long long 
    using namespace std;
    
    const int N=255;
    const int M=55;
    const int mod=1e9+9;
    
    int n,m,T;
    int f[M][N][N<<2];
    inline void add(int&x,int y) {
        if((x+=y)>=mod) x-=mod;
    }
    
    int main() {
        scanf("%d%d",&n,&m);
        for(T=1; T<n; T<<=1);
        f[0][0][0]=1;
        for(int j=0; j<m; ++j) {
            int(&F)[N][N<<2]=f[j+1];
            int(&G)[N][N<<2]=f[j];
            if(j&1) {
                for(int i=j; i<=n; ++i) 
                for(int p=0; p<T; ++p) if(G[i][p])
                for(int k=i+1; k<=n; ++k) 
                    add(F[k][p],G[i][p]);
            }
            else {
                for(int i=j; i<=n; ++i) 
                for(int p=0; p<T; ++p) if(G[i][p])
                for(int k=i+1; k<=n; ++k) 
                    add(F[k][p^(k-i-1)],G[i][p]);
            }
        }
        int ans=0;
        for(int i=m; i<=n; ++i) 
        for(int p=1; p<T; ++p) 
            add(ans,f[m][i][p]);
        printf("%d
    ",ans);
        return 0;
    }
    
  • 相关阅读:
    【 数据结构(C语言)】栈的应用——行编辑程序
    【 数据结构 (C语言)】栈的应用(二)——括号匹配问题
    节点
    页面加载--延迟加载
    雅黑php 探针
    Swiper 触屏滑动切换
    tab 选择悬停展示
    翻牌抽奖功能讲解
    公告信息滚动功能
    织梦提交表单不进行跳转
  • 原文地址:https://www.cnblogs.com/nosta/p/10982021.html
Copyright © 2020-2023  润新知