bzoj 3209 花神的数论题 —— 数位DP

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=3209

算是挺简单的数位DP吧，但还是花了好久才弄明白...

又参考了博客：https://blog.csdn.net/sunshinezff/article/details/51049132

先预处理开头是 0/1 的 i 位数中有 j 个1的数的数量，然后按位一边限制大小，一边分成0开头或1开头统计答案。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
ll mx,f[][],g[][],ans=,mod=;
ll pw(ll a,ll b)
{
    ll ret=;
    for(;b;b>>=1ll,a=(a*a)%mod)
        if(b&)ret=(ret*a)%mod;
    return ret;
}
ll cal(ll n)
{
    int t=,c=;
    for(;(1ll<<t)<=n;t++);//t是最高位
    for(;t;t--)//当前是t位数
        if((1ll<<(t-))&n)//限制大小
        {
            for(int i=;i<=t;i++)
                (ans*=pw(i+c,f[t][i]))%=mod;//本位取0
            if(c) (ans*=c)%=mod; c++;//c表示前面已经有几个1 //加上本位取1的一种情况
        }
    return (ans*c)%mod;
}
int main()
{
    scanf("%lld",&mx);
    f[][]=; g[][]=;
    for(int i=;i<=;i++)
        for(int j=;j<=i;j++)
        {
            f[i][j]=f[i-][j]+g[i-][j];
            if(j)g[i][j]=g[i-][j-]+f[i-][j-];
        }
    printf("%lld\n",cal(mx));
    return ;
}