花神的数论题

分析：就是对于1~n中的每个数，把他二进制中含有的1的个位乘起来即可。n是小于1e15的，所以n的二进制位也不会超过70，所以说1~n中的所有数含有1的个数是1~70，所以对于k属于区间[1,70]，我们可以用数位DP求出二进制数中含有1的个数为k的数的数目sum，这样k的贡献就是sum个k连乘，也就是k的sum次方，这样我们枚举每一个二进制表示中可能含有的1的个数就可以得到最终的贡献，而求每一个二进制表示中可能含有的1的个数就是一个基本的数位DP，我们只需要在数位DP过程中记录当前已经出现多少个1以及最后需要出现多少个1即可。

但是需要注意两个问题：

（1）对于满足同一个条件的多组询问，则我们上次记忆化的f数组依旧能用于下一次，所以f数组只需要初始化一次，而本题每次数位DP满足的条件都不一样（二进制表示中含有的1的个数不同）,所以我们每次都需要重新初始化f数组

（2）因为答案涉及到取模，但是需要注意的是我们仅仅能在加法减法或者乘法过程中随意取模，但是我们不能在求指数的过程中取模，否则会出现错误
——————————————

 上代码：

             #include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
using namespace std;
typedef long long ll;
const int N=70,mod=10000007;
ll f[N][N];//表示后i位可以任选且前i位含有j个1的个数 
ll a[N];
ll qpow(ll a,ll b)
{
    ll ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%mod;
        a=a%mod*a%mod;
        b>>=1;
    }
    return ans;
}
ll dp(int pos,int cnt,int limit,int t)//因为求的是指数，所以在求DP过程中不能模mod 
{
    if(!pos) return cnt==t;
    if(cnt>t) return 0;//剪枝 
    if(!limit&&f[pos][cnt]!=-1) return f[pos][cnt];
    ll ans=0;
    int up=limit?a[pos]:1;
    for(int i=0;i<=up;i++)
        ans+=dp(pos-1,cnt+(i==1),limit&&(i==up),t);
    if(!limit) f[pos][cnt]=ans;
    return ans;
}
ll solve(ll x)
{
    ll ans=1;
    int pos=0;
    while(x)
    {
        a[++pos]=x%2;
        x/=2;
    }
    for(int i=1;i<=pos;i++)
    {
        memset(f,-1,sizeof f);//因为我们每次求的满足条件不一样，所以每次都需要初始化 
        ans=ans*(qpow(i,dp(pos,0,1,i)))%mod;
    }
    return ans;
}
int main()
{
    ll n;
    cin>>n;
    cout<<solve(n);
    return 0;
}