参考链接:
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;
}