洛谷P4317 花神的数论题

最初想法是打表找规律，当成数学题来做，看了这里之后发现跟自己的想法很像（以前好像做过类似题目），就是通过组合数求出$n$以内的$sum(i)$出现个数。因为模数不是质数，所以要用到欧拉定理

$$a^{\phi p}\equiv1\pmod p$$

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int mo=10000007;
ll n,ans,c[61][61];
inline ll F(ll x){
    return x>=mo?x-mo:x;
}
inline ll qp(ll x,ll t){
    ll res=1;
    while(t){
        if(t&1) res=res*x%mo;
        x=x*x%mo;t>>=1;
    }return res;
}
void init(){
    for(ll i=0;i<=60;++i)c[i][0]=1;
    for(ll i=1;i<=60;++i)
    for(ll j=1;j<=i;++j)
    c[i][j]=(c[i-1][j-1]+c[i-1][j])%9988440;
}
int main(){
    scanf("%lld",&n);
    init();ans=1;ll t=0;
    for(ll i=60;i>=0;--i){
        if(!((n>>i)&1))continue;
        for(ll j=0;j<=i;++j)
        if(j+t)ans=ans*qp(j+t,c[i][j])%mo;
        ++t;
    }
    return !printf("%lld",ans*t%mo);
}