喵哈哈村的秘境探险(四)
发布时间: 2017年4月25日 20:24 最后更新: 2017年4月25日 20:25 时间限制: 1000ms 内存限制: 128M
喵哈哈村的一堆人在前往北京的路上,发现了一个洞穴。由于好奇心大作,于是准备前往洞穴进行探险。
走啊走,喵哈哈村一群人发现了一个先知,这个先知想知道:
设 sum(i) 表示 i 的二进制表示中 1 的个数。给出一个正整数 N ,先知要问你,sum(1)到sum(n)的乘积是多少。
一个正整数 N。
N<=1e15
一个数,答案模 10000007 的值。
复制
8
24
当前位是1的话,假如把它变成0,呢低位便可以随意取,呢么组合数一发就好了。。
官方题解:http://www.cnblogs.com/qscqesze/p/6760700.html
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
typedef long long ll;
#define Mod 10000007
ll a[65],C[65][65],len;
void init()
{
int i,j;
for(i=0;i<60;i++)
C[i][0]=1;
for(i=1;i<60;i++)
for(j=1;j<=i;j++)
C[i][j]=C[i-1][j]+C[i-1][j-1];
}
ll Pow(ll x,ll y)
{
ll res=1;
while(y)
{
if(y%2)
res=res*x%Mod;
x=x*x%Mod;
y/=2;
}
return res;
}
ll work(ll x)
{
int i;
ll sum=0;
for(i=len;i>0;i--)
{
if(a[i])
{
sum+=C[i-1][x];
x--;
}
if(x<0)
break;
}
return sum;
}
int main(void)
{
init();
ll n,i,tmp,ans=1;
scanf("%lld",&n);
tmp=n+1;
while(tmp)
{
a[++len]=tmp%2;
tmp>>=1;
}
for(i=1;i<=len;i++)
ans=ans*Pow(i,work(i))%Mod;
printf("%lld\n",ans);
return 0;
}