关于数位DP

#include<iostream>
#define int long long
using namespace std;
const int mod=10000007;
int n,digit[65],top,ans=1,c[65][65];
int quickpow(int base,int exp){
    int tmp=base,ans=1;
    for(int i=0;((long long)1<<i)<=exp;i++){
        if(((long long)1<<i)&exp) (ans*=tmp)%=mod;
        (tmp*=tmp)%=mod;
    }
    return ans;
}
signed main(){
    cin>>n;
    ans=1;
    c[0][0]=1;
    for(int i=1;i<=62;i++){
        for(int j=0;j<=i;j++){
            c[i][j]=c[i-1][j]+c[i-1][j-1];
        }
    }
    while(n){
        digit[++top]=n%2;
        n/=2;
    }
    for(int i=2;i<=top;i++){
        int tmp=i;
        for(int j=top;j>=1;j--){
            if(digit[j]){
                (ans*=quickpow(i,c[j-1][tmp]))%=mod;
                tmp--;
            }
            if(!tmp){
                (ans*=i)%=mod;
                break;
            }
        }
    }
    cout<<ans;
}