bzoj3209

参考链接：

http://blog.csdn.net/popoqqq/article/details/39936291

感觉数位dp什么的，只要把它分解成一个个的函数，就可以有比较清楚的思路。注意一下细节就可以了。

（省赛时，遇到过数位dp，那时候还不知道是这么个东西，但是凭着直觉，就硬写了下来，写的思路一团糟，但是幸亏ac了，后来知道这是数位dp！）

（意外？不，在加图索家的家训里没有意外这回事，意外只是懦夫为自己找的理由。   ——江南 《龙族》）

2015.10.27

今天的晚饭是小土豆和西北风。

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<iostream>
using namespace std;
#define N 60
#define MOD 10000007

long long int c[N][N];
long long int ans;

long long int km(long long int x,long long int y){
    long long int ans=1;

    while(y){
        if(y%2){
            ans=ans*x%MOD;
        }
        y=y/2;
        x=x*x%MOD;
    }

    return ans;
}

void getpro(int sum,int templength,int length){
    int start=0;

    if(length-1==templength){
        start=1;
    }
    for(int i=start;i<=templength;i++){
        ans=ans*km(sum+i,c[templength][i])%MOD;
    }

    return;
}

int main(){
    for(int i=0;i<N;i++){
        c[i][0]=1;
        c[i][i]=1;
        for(int j=1;j<i;j++){
            c[i][j]=c[i-1][j-1]+c[i-1][j];
        }
    }

    long long int n;
    long long int tempn;

    while(scanf("%lld",&n)!=EOF){
        int length=log2(n)+1;
        int sum=0;

        ans=1;
        for(int i=length-1;i>=0;i--){
            tempn=n>>i;
            if(tempn%2){
                getpro(sum,i,length);
                sum++;
            }
        }
        ans=ans*sum%MOD;//忘了这一步，得出的答案是不对的。n本身并没有通过上面的for循环计算进去，所以要在这里补充一下。

        printf("%lld\n",ans);
    }

    return 0;
}