首页 > 代码库 > ACdream 1154 Lowbit Sum (数位DP)

ACdream 1154 Lowbit Sum (数位DP)

Lowbit Sum

Time Limit: 2000/1000MS (Java/Others)Memory Limit: 128000/64000KB (Java/Others)
Submit   Status

Problem Description

long long ans = 0;
for(int i = 1; i <= n; i ++)
    ans += lowbit(i)
lowbit(i)的意思是将i转化成二进制数之后,只保留最低位的1及其后面的0,截断前面的内容,然后再转成10进制数
比如lowbit(7),7的二进制位是111,lowbit(7) = 1
6 = 110(2),lowbit(6) = 2,同理lowbit(4) = 4,lowbit(12) = 4,lowbit(2) = 2,lowbit(8) = 8

每输入一个n,求ans

Input

多组数据,每组数据一个n(1 <= n <= 10^9)

Output

每组数据输出一行,对应的ans

Sample Input

1
2
3

Sample Output

1
3
4



大致题意:中文的,都能看懂吧。。。


解题思路:这里利用了数论的一些小技巧,关于lowbit的规律。我发现这类题数论题上来直接暴力严重超时的,一般来说都有规律!!!

开始做的时候直接暴力打表,结果打个表都跑了十几秒,还是算了。。。

然后就潜心于找规律了,先随便输出了从1开始的几个连续数的lowbit值,还没啥感觉,后来又多输出了几组,才渐渐发现了规律——奇数的lowbit都是1;偶数的lowbit是先增后减的而且还是对称的,并且从两边向中间看的话,都是公比为2的等比数列。这样就可以计算了。如果n为偶数,偶数的可以转化为2*dp[n/2],然后再加上奇数的n/2个1,就可以了; n为奇数时,偶数的还是转化成2*dp[n/2],但是奇数的现在不是n/2个了,而是n/2 + 1个了。要想方便的总结一下,就可以写成dp[n] = 2*dp[n/2] + n/2 + (n%2);但是最近又发现了一种新的写法,那就是位运算的写法,位运算也可以实现乘除,而且比乘除运算要快,当然也能判别一个数的奇偶,可能是因为计算机本来就只能识别0和1的缘故吧,这些位运算就是直接对二进制数操作,所以更快。于是状态转移方程就可以写成dp[n] = 2*dp[n>>1] + (n>>1) + (n&1).




AC代码:

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

long long dp(int x){
	if(x == 1)  return 1;
	return  2*dp(x>>1) + (x>>1) + (x&1);
}

int main(){
	#ifndef ONLINE_JUDGE
	freopen("in.txt","r",stdin);
	#endif
	int n;
	while(~scanf("%d",&n)){
		printf("%lld\n", dp(n));
	}
	return 0;
}


ACdream 1154 Lowbit Sum (数位DP)