[Bzoj3209]花神的数论题（数位dp）

3209: 花神的数论题

---

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 2633  Solved: 1182
[Submit][Status][Discuss]

Description

---

背景

众所周知，花神多年来凭借无边的神力狂虐各大 OJ、OI、CF、TC …… 当然也包括 CH 啦。
描述
话说花神这天又来讲课了。课后照例有超级难的神题啦…… 我等蒟蒻又遭殃了。
花神的题目是这样的
设 sum(i) 表示 i 的二进制表示中 1 的个数。给出一个正整数 N ，花神要问你
派(Sum(i)),也就是 sum(1)—sum(N) 的乘积。

Input

---

一个正整数 N。

Output

---

一个数，答案模 10000007 的值。

Sample Input

---

样例输入一

Sample Output

---

样例输出一

HINT

---

对于样例一，1*1*2=2；

数据范围与约定

对于 100% 的数据，N≤10^15

Source

---

原创 Memphis

题解：

---

求一个数二进制有多少个1，再加上数据范围， 很显然让我们做数位dp。

然后我们发现不能一起转移，但我们可以求出<=N中 1的个数为k的所有数

 

最终答案就是。

关于dfs那个转移，就是很裸的模板题，看代码即可。

AC代码：

---

# include <iostream>
# include <cstdio>
# include <cstring>
using namespace std;
typedef long long LL;
const LL mod = ;
LL f[][],n;int data[],len;
LL dfs(int now,int num,bool lim,bool first)
{
    if(num < )return ;
    if(!now)return !num;
    if(!lim && !first && f[now][num] != -)return f[now][num];
    LL ret = ;int p = lim ? data[now] : ;
    for(int i = ;i <= p;i++)
    ret += dfs(now - ,num - i,lim && i == p,first && !i);
    if(!lim && !first)f[now][num] = ret;
    return ret;
}
LL cmd(LL k,LL x)
{
    LL c = ;
    while(k)
    {
        if(k & 1LL)c = c * x % mod;
        x = x * x % mod;
        k >>= 1LL;
    }
    return c;
}
LL calc(LL k)
{
    memset(f,-,sizeof f);
    len = ;
    while(k)data[++len] = k % ,k /= ;
    LL ret = ;
    for(int i = ;i <= len;i++)
    ret = ret * cmd(dfs(len,i,true,true),i) % mod;
    return ret;
}
int main()
{
    scanf("%lld",&n);
    printf("%lld\n",calc(n));
}