这道题细节并不算太多,但是求方案数的时候一定不要乱取模!
如果非要取模的话也要遵循欧拉定理.
code:
#include <bits/stdc++.h> #define N 57 #define ll long long #define mod 10000007 #define setIO(s) freopen(s".in","r",stdin) using namespace std; ll f[N][2][N]; int num[N],cnt; int qpow(int x,ll y) { int tmp=1; for(;y;y>>=1,x=(ll)x*x%mod) if(y&1) tmp=(ll)tmp*x%mod; return tmp; } void init() { f[1][0][0]=f[1][1][1]=1; for(int i=2;i<=56;++i) { for(int j=0;j<=i;++j) { f[i][0][j]+=f[i-1][0][j]+f[i-1][1][j]; if(j) f[i][1][j]+=f[i-1][0][j-1]+f[i-1][1][j-1]; } } } ll calc(ll x,int t) { memset(num,0,sizeof(num)),cnt=0; while(x) { num[++cnt]=x&1; x>>=1; } // 先计算不满足 cnt 位的所有数字. ll ans=0; int cn=0; for(int i=1;i<cnt;++i) ans+=f[i][1][t]; for(int i=cnt;i>=1;--i) { if(num[i]&&i!=cnt&&cn<=t) ans+=f[i][0][t-cn]; cn+=num[i]; } return ans; } int main() { // setIO("input"); ll n,ans=1; init(); scanf("%lld",&n); for(int i=2;i<=56;++i) { ans=ans*1ll*qpow(i,calc(n+1,i))%mod; } printf("%lld\n",ans); return 0; }