ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

K - Binary Indexed Tree HDU - 5921 (数位dp)

2021-02-21 20:01:24  阅读:186  来源: 互联网

标签:Binary HDU 150 int ll Tree pos long define


#include<stdio.h>
#include<string.h>
#define ll long long
ll w[150],q[150],dp[150][2][2][2];
#define mod 1000000007
ll getans(int pos,int flaga,int flagb){
	if(pos<0)return 1;
	ll ans=1;
	if(flaga==1)ans=(ans*q[pos])%mod;
	else ans=(ans*((1ll<<(pos+1))%mod))%mod;
	if(flagb==1)ans=(ans*q[pos])%mod;
	else ans=(ans*((1ll<<(pos+1))%mod))%mod;
	return ans;
}
ll dfs(int pos,int flaga,int flagb,int flag){
	if(pos<0)return 0;
	if(dp[pos][flaga][flagb][flag]!=-1)return dp[pos][flaga][flagb][flag];
	int maxa=1,maxb=1;
	if(flaga==1)maxa=w[pos];
	if(flagb==1)maxb=w[pos];
	ll ans=0;
	for(int i=0;i<=maxa;i++){
		for(int j=0;j<=maxb;j++){
			ans=(ans+dfs(pos-1,flaga&(i==maxa),flagb&(j==maxb),flag|(i!=j)))%mod;
			ans=(ans+(flag|(i!=j))*(i+j)*getans(pos-1,flaga&(i==maxa),flagb&(j==maxb)))%mod;
		}
	}
	dp[pos][flaga][flagb][flag]=ans;
	return ans;
}
ll solve(ll n){
	int cnt=0;
	while(n!=0){
		w[cnt++]=n%2;
		n=n/2;
	}
	q[0]=w[0];
	for(int i=1;i<cnt;i++){
		q[i]=(q[i-1]+(w[i]<<i))%mod;
	}
	for(int i=0;i<cnt;i++){
		q[i]++;
	}
	return dfs(cnt-1,1,1,0);
}
int main(){
	int t;
	scanf("%d",&t);
	int cas=0;
	while(t--){
		memset(dp,-1,sizeof(dp));
		ll n;
		scanf("%lld",&n);
		ll temp=(mod+1)/2;
		ll ans=solve(n)*temp%mod;
		printf("Case #%d: %lld\n",++cas,ans);
	}
}

标签:Binary,HDU,150,int,ll,Tree,pos,long,define
来源: https://www.cnblogs.com/League-of-cryer/p/14426915.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有