BZOJ 3209 数位DP

思路：
先预处理出来组合数
按位做 枚举sum[x]是多少

注意Mod不是一个质数

//By SiriusRen
#include <cstdio>
using namespace std;
#define int long long
int n,top,w[66],C[66][66],ans=1,cnt,Mod=10000007;
int pow(int x,int y){
    int res=x,tmp=1;
    while(y){
        if(y&1)tmp=(tmp*res)%Mod;
        res=(res*res)%Mod;
        y>>=1;
    }return tmp;
}
signed main(){
    scanf("%lld",&n);
    for(;n;n>>=1)w[++top]=n&1;
    for(int i=0;i<=65;i++)C[i][0]=1;
    for(int i=1;i<=65;i++)
        for(int j=1;j<=65;j++)
            C[i][j]=C[i-1][j-1]+C[i-1][j];
    for(;top;top--)if(w[top]){
        for(int i=(cnt==0);i<top;i++)
            ans=ans*pow(cnt+i,C[top-1][i])%Mod;
        cnt++;
    }printf("%lld\n",(ans*cnt)%Mod);
}