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 10000007long 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;}