[BZOJ3209]花神的数论题

最新推荐文章于 2020-06-10 21:19:08 发布

原创最新推荐文章于 2020-06-10 21:19:08 发布 · 539 阅读

1 ·

CC 4.0 BY-SA版权

BZOJ 同时被 3 个专栏收录

159 篇文章

订阅专栏

31 篇文章

订阅专栏

数位DP

4 篇文章

订阅专栏

本文深入探讨了数位DP算法的应用实例，通过具体代码详细解释了如何利用数位DP解决特定问题，包括状态转移方程的设计及实现细节。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

原题地址

数位DP…
感觉已经被数位DP虐习惯了…

#include <cstdio>
typedef long long ll;
const ll N=101;
const ll TOMOD=10000007;
ll a,len,Ans=1;
ll b[N],ans[N];
ll f[N][N];

void cal(ll x,ll tag){
    if(x==1){
        if(b[x]) ans[tag+1]++;
        ans[tag]++;
        return ;
    }
    if(!b[x]){
        cal(x-1,tag);
        return ;
    }
    for(ll i=0;i<=x-1;i++) ans[i+tag]+=f[x-1][i];
    cal(x-1,tag+1);
}

ll qpow(ll x,ll p){
    if(!p) return 1;
    ll y=qpow(x,p>>1);
    y=(y*y)%TOMOD;
    if(p&1) y=(y*x)%TOMOD;
    return y;
}

int main(){
    scanf("%lld",&a);
    while(a){
        b[++len]=a&1;
        a>>=1;
    }
    f[1][0]=f[1][1]=1;
    for(ll i=2;i<=len;i++){
        for(ll j=0;j<=len;j++){
            f[i][j]=f[i-1][j];
            if(j) f[i][j]+=f[i-1][j-1];
        }
    }
    cal(len,0);
    for(ll i=1;i<=len;i++) Ans=(Ans*qpow(i,ans[i]))%TOMOD;
    printf("%lld\n",Ans);

    return 0;
}