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


posted @ 2015-10-09 20:27  buzhidaohahaha  阅读(220)  评论(0编辑  收藏  举报