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

正解：數位$dp$。

數位$dp$水題。設$f[i][j][k]$表示前$i$位，二進制和為$j$，當前這位有沒有限制的答案。

直接用記憶化搜索轉移即可。

 1 #include <bits/stdc++.h>
 2 #define il inline
 3 #define RG register
 4 #define ll long long
 5 #define
 
 rhl (10000007)
 6 
 7 using namespace std;
 8 
 9 int f[70][70][2],st[70],top;
10 ll n;
11 
12 il int gi(){
13   RG int x=0,q=1; RG char ch=getchar();
14   while ((ch<‘0‘ || ch>‘9‘) && ch!=‘-‘) ch=getchar();
15   if (ch==‘-‘) q=-1,ch=getchar();
16   while (ch>=‘0‘ && ch<=‘
 
9‘) x=x*10+ch-48,ch=getchar();
17   return q*x;
18 }
19 
20 il int dfs(RG int pos,RG int tot,RG int lim){
21   if (!pos) return tot?tot:1;
22   if (f[pos][tot][lim]!=-1) return f[pos][tot][lim];
23   RG int end=lim?st[pos]:1,&res=f[pos][tot][lim]; res=1;
24   for (RG int i=0;i<=end;++i)
25     res=1LL*res*dfs(pos-1,tot+i,lim && i==end)%rhl;
26   return res;
27 }
28 
29 int main(){
30 #ifndef ONLINE_JUDGE
31   freopen("god.in","r",stdin);
32   freopen("god.out","w",stdout);
33 #endif
34   cin>>n; memset(f,-1,sizeof(f));
35   while (n) st[++top]=n&1,n>>=1;
36   cout<<dfs(top,0,1); return 0;
37 }

bzoj3209 花神的數論題