喵哈哈村的秘境探险（四）-(组合数)

喵哈哈村的秘境探险（四）

发布时间: 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 的值。

样例输入1  复制

8

样例输出1

24

查看隐藏信息

选择语言

 C (GCC 4.8)  C++ (G++ 4.3)  Java (Oracle JDK 1.7)

当前位是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;
}