[bzoj3209]花神的数论题

[bzoj3209]花神的数论题

---

DP[i][1/0][j]表示考虑到二进制第i位，是否顶上界，1有j个的方案数

然后快速幂算答案

• 代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int m=60;
ll n;

ll DP[80][2][80];
int num[80];
ll mod=10000007;

ll qpow(ll x,ll b){
	ll ret=1;
	while(b){
		if(b%2)	ret=(ret*x)%mod;
		b/=2;
		x=(x*x)%mod;
	}
	return ret;
}

inline void dp(){
	for(int i=1;i<=60;i++,n/=2) num[i]=n%2;
	DP[60][1][0]=1ll;
	for(int i=59;i;i--){
		if(num[i]){
			for(int j=0;j<=m;j++){//选0
				DP[i][0][j]+=DP[i+1][0][j]+DP[i+1][1][j];
			}
			for(int j=1;j<=m;j++){//选1
				DP[i][1][j]=DP[i+1][1][j-1];
				DP[i][0][j]+=DP[i+1][0][j-1];
			}
		}else{
			for(int j=0;j<=m;j++){//选0
				DP[i][0][j]+=DP[i+1][0][j];
				DP[i][1][j]+=DP[i+1][1][j];
			}
			for(int j=1;j<=m;j++){//选1
				DP[i][0][j]+=DP[i+1][0][j-1];
			}
		}
	}
	ll ans=1;
	ll p=1;
	for(int i=1;i<=m;++i,++p){
		ll num=DP[1][0][i]+DP[1][1][i];
		ans=(ans*qpow(p,num))%mod;
	}
	printf("%lld\n",ans);
}

int main()
{
	scanf("%lld",&n);
	dp();
}