求含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;
}