[bzoj3209]花神的数论题

[bzoj3209]花神的数论题

---

DP[i][1/0][j]表示考虑到二进制第i位，是否顶上界，1有j个的方案数

然后快速幂算答案

• 代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int m=60;
ll n;

ll DP[80][2][80];
int num[80];
ll mod=10000007;

ll qpow(ll x,ll b){
    ll ret=1;
    while(b){
        if(b%2) ret=(ret*x)%mod;
        b/=2;
        x=(x*x)%mod;
    }
    return ret;
}

inline void dp(){
    for(int i=1;i<=60;i++,n/=2) num[i]=n%2;
    DP[60][1][0]=1ll;
    for(int i=59;i;i--){
        if(num[i]){
            for(int j=0;j<=m;j++){//选0
                DP[i][0][j]+=DP[i+1][0][j]+DP[i+1][1][j];
            }
            for(int j=1;j<=m;j++){//选1
                DP[i][1][j]=DP[i+1][1][j-1];
                DP[i][0][j]+=DP[i+1][0][j-1];
            }
        }else{
            for(int j=0;j<=m;j++){//选0
                DP[i][0][j]+=DP[i+1][0][j];
                DP[i][1][j]+=DP[i+1][1][j];
            }
            for(int j=1;j<=m;j++){//选1
                DP[i][0][j]+=DP[i+1][0][j-1];
            }
        }
    }
    ll ans=1;
    ll p=1;
    for(int i=1;i<=m;++i,++p){
        ll num=DP[1][0][i]+DP[1][1][i];
        ans=(ans*qpow(p,num))%mod;
    }
    printf("%lld\n",ans);
}

int main()
{
    scanf("%lld",&n);
    dp();
}