【bzoj 3209】花神的数论题（组合数学）

传送门biu~
令$an_i$代表二进制中$1$的个数为$i$的数有多少个。从高位向低位枚举，假设输入的数$n$的二进制第$i$位为$1$，那么设一个数$x$的第$i$位为$0$，前$i-1$位与$n$相同，那么无论比$i$位低的数位中哪一位为$1$，$x$都不会大于$n$，所以可以用组合数$C_{i-1}^k$更新答案。
Ps：$\varphi(10000007)=9988440$
说得不清楚看代码：

#include<bits/stdc++.h>
using namespace std;
const long long mod=10000007,phi=9988440;
long long n,ans=1;
long long c[51][51],an[51];
int a[51],now,cnt;
inline void getC(){
    for(int i=0;i<=50;++i){
        c[i][0]=1;
        for(int j=1;j<=i;++j)  c[i][j]=c[i-1][j]+c[i-1][j-1];
    }
}
inline long long ksm(long long p,long long n){
    long long re=1;
    while(n){
        if(n&1) re=re*p%mod;
        p=p*p%mod;
        n>>=1;
    }
    return re;
}
int main(){
    scanf("%lld",&n);
    getC();
    while(n){a[++cnt]=n&1;n>>=1;}
    for(int i=cnt;i>=1;--i){
        if(!a[i])   continue;
        for(int j=1;j<i;++j)    an[now+j]+=c[i-1][j];
        ++an[++now];
    }
    for(int i=1;i<=cnt;++i)     (ans*=ksm(i,an[i]%phi))%=mod;
    printf("%lld",ans);
    return 0;
}