洛谷P4317 花神的数论题

版权声明:欢迎转载(标记出处),写得不好还请多指教 https://blog.csdn.net/quan_tum/article/details/82154428

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

a ϕ p 1 ( mod 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);
}

猜你喜欢

转载自blog.csdn.net/quan_tum/article/details/82154428