P4317 花神的数论题

求含1个数的乘积，也可以遍历从1个到二进制有效最长的长度，分别用快速幂求出每部分的值
mst成-1能过 但是mst成0就过不了
十分神奇

//cyc
#pragma GCC optimize("Ofast")
#pragma GCC target("avx,avx2,fma")
#pragma GCC optimization ("unroll-loops")
#include<bits/stdc++.h>
#define rep(i,a,n) for(int i=a;i<=n;i++)
#define per(i,a,n) for(int i=n;i>=a;i--)
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define mst(a) memset(a,-1,sizeof a)
#define int long long
using namespace std;
typedef pair<int,int> pii;

/*
 ____ ____   _____ .____________   _____
/_   /_   | /  |  ||   ____/_   | /  |  |
 |   ||   |/   |  ||____  \ |   |/   |  |_
 |   ||   /    ^   /       \|   /    ^   /
 |___||___\____   /______  /|___\____   |
               |__|      \/          |__|

*/
const int mod = 10000007;
int a[63];
int dp[64][64];//pos cnt

int q_pow(int x,int n)
{
    int res=1;
    while(n)
    {
        if(n&1)
        {
            res=res*x%mod;
        }
        x=x*x%mod;
        n>>=1;
    }
    return res;
}
int dfs(int pos,int now,int cnt,bool tp)//now 目前是1的位数 cnt目标1的个数
{
    if(pos==-1){
        return now==cnt;
    }
    if(!tp&&dp[pos][now]!=-1) return dp[pos][now];
    int ret=0;
    int bound=tp?a[pos]:1;
    for(int i=0;i<=bound;i++){
        ret+=dfs(pos-1,now+(i==1),cnt,tp&&(i==bound));
    }
    if(!tp)dp[pos][now]=ret;
    return ret;
}
signed main()
{
    ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    int n;
    cin>>n;
    int c=0;
    while(n){
        a[c]=n&1;
        n>>=1;
        c++;
    }
   // cout<<c<<endl;
    int ans=1;
    for(int i=1;i<=c;i++){
        mst(dp);
        ans=ans*q_pow(i,dfs(c-1,0,i,1))%mod;
    }
    cout<<ans<<endl;

}