Description
背景 众所周知,花神多年来凭借无边的神力狂虐各大 OJ、OI、CF、TC …… 当然也包括 CH 啦。 描述
话说花神这天又来讲课了。课后照例有超级难的神题啦…… 我等蒟蒻又遭殃了。 花神的题目是这样的 设 sum(i) 表示 i 的二进制表示中 1
的个数。给出一个正整数 N ,花神要问你 派(Sum(i)),也就是 sum(1)—sum(N) 的乘积。
Input
一个正整数 N。
Output
一个数,答案模 10000007 的值。
Sample Input
样例输入一
3
Sample Output
样例输出一
2
HINT
对于样例一,1*1*2=2;
数据范围与约定
对于 100% 的数据,N≤10^15
题解
打了个表后发现有点小小的意思。。
枚举1的个数i,那么可以发现有很多数的二进制1的个数是相同的,这可以搞个dp
设f[i][j]为二进制位数i,1的个数为j,首位是0的个数
g[i][j]为二进制位数i,1的个数为j,首位是1的个数
转移瞎搞
加的时候可以枚举一下n的二进制位,这个位是1的时候就进去搞一搞,因为我们不想超过范围对吧
快速幂get完就OK了
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
const LL mod=10000007;
LL pow_mod(LL a,LL b)
{
LL ret=1;
while(b)
{
if(b%2==1)ret=ret*a%mod;
a=a*a%mod;b/=2;
}
return ret;
}
LL f[110][110],g[110][110];
/*
0开头 位数i 1的个数为j
1开头 位数i 1的个数为j
*/
LL n,mx;
LL sol(LL x)
{
LL t=mx,c=0,ans=1;
while(t)
{
if((1LL<<(t-1)&x))
{
for(int i=1;i<=t;i++)
ans=ans*pow_mod(i+c,f[t][i])%mod;
if(c)ans=ans*c%mod;
c++;
}
t--;
}
return (ans*c)%mod;
}
int main()
{
scanf("%lld",&n);
mx=0;LL tmp=n;
while(tmp){mx++;tmp/=2;}
g[1][1]=1;f[1][0]=1;
for(int i=2;i<=mx;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]=g[i-1][j-1]+f[i-1][j-1];
}
printf("%lld\n",sol(n));
return 0;
}