bzoj3209: 花神的数论题（数位Dp）

题目传送门
神。

解法：
又是肉dalao推荐的神题。
想了半天发现就算我知道每个的sum我要求派还是会超时啊。
于是我根本没有往同一个sum可能有多个的地方想。
结果正解就是这样的。。

数位Dp的思想：
f[i][j]表示二进制有i位，最高位为0，一共有j个1
g[i][j]表示二进制有i位，最高位为1，一共有j个1。
那么转移就很容易啦。

然后由于有n的限制。
每到一位的时候我们默认前面的位都跟n的一样。
那么我们从最高位开始。如果当前位为1的话。
我们就有两种选择。填0或1。具体看代码吧。

代码实现：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
ll g[110][110],f[110][110];
ll len;ll n;const ll mod=10000007;
ll pow_mod(ll a,ll b) {
    ll ans=1;
    while(b!=0) {
        if(b%2==1)ans=(ans*a)%mod;
        b/=2;a=(a*a)%mod;
    }return ans;
}
int main() {
    scanf("%lld",&n);
    len=0;ll x=n;while(x!=0) {x/=2;len++;}
    f[1][0]=1;g[1][1]=1;
    for(int i=2;i<=len;i++) 
        for(int j=0;j<=i;j++) {
            f[i][j]=f[i-1][j]+g[i-1][j];
            if(j!=0)g[i][j]=f[i-1][j-1]+g[i-1][j-1];
        }

    ll ans=1;ll c=0;  //c表示前面的位有多少个1
    for(ll t=len;t>=1;t--) {
        if((1ll<<(t-1)&n)) {
            for(int i=1;i<=t;i++) ans=(ans*pow_mod(i+c,f[t][i]))%mod;  //后面的位填0保证不超过了。
            if(c) ans=(ans*c)%mod;c++;  //因为前面记录的是后面的位一定有1的，那么还有一种情况就是后面全部都是0.那么整个数只有前面有1，也就是c个1。
        }
    }printf("%lld\n",(ans*c)%mod);
    return 0;
}