LightOJ 1032 Fast Bit Calculations

题目链接：http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=25909

数位DP专题：http://blog.csdn.net/chy20142109/article/details/50930004

题意：给定范围内，所有数的二进制表示中含连续1的位的总数，比如 01110 第3,4位就是符合计数要求的位，如果一个位是1，且它的前一位也是1，那么它
就算一个。

思路：可以直接递推解决。
dp[i][0] = dp[i-1][0] + dp[i-1][1];
dp[i][1] = dp[i-1][0] + dp[i-1][1] + 2^(i-2);
因为从第i位开始是11那么后面有i-2位，可以随便取，所以"11"就会出现2^(i-2)次。
统计[0,x]范围内的数
对于第一类计数就是dp[len][0]
第二类计数要注意一点就是处理到第i位时，如果i+1位也是1，那么就要额外加上k次，k = x-∑（bit[j]*2^(j-1)）j>=i
也就是说前面的数位已经固定了，后面几位从全部0开始到上限情况还有多少种取值方式。    

#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <sstream>
#include <queue>
#include <utility>
using namespace std;

#define rep(i,j,k) for (int i=j;i<=k;i++)
#define Rrep(i,j,k) for (int i=j;i>=k;i--)

#define Clean(x,y) memset(x,y,sizeof(x))
#define LL long long
#define ULL unsigned long long
#define inf 0x7fffffff
#define mod %100000007

int T;
LL n;
LL dp[35][2];
LL Pow[40];

void init()
{
    Clean(dp,0);
    LL t = 0;
    Pow[0] = 1;
    rep(i,1,32) Pow[i] = Pow[i-1]<<1;

    // dp[i][0] = dp[i-1][0] + dp[i-1][1];
    // dp[i][1] = dp[i-1][0] + dp[i-1][1] + 2^(i-2);
    rep(i,2,31)
    {
        dp[i][1] = t + Pow[i-2];
        dp[i][0] = t;
        t = dp[i][1] + dp[i][0];
    }
}

LL cal(LL x)
{
    LL temp = x;
    int bit[35];
    Clean(bit,0);
    int len = 0;
    while( x )
    {
        bit[++len] = x % 2;
        x /= 2;
    }
    LL ans = 0;
    ans+=dp[len][0];

    temp-=Pow[len-1];
    Rrep(i,len-1,1)
    {
        temp-=bit[i]*Pow[i-1];

        if ( bit[i] )
        {
            ans+=dp[i][0];
            if ( bit[i+1] ) ans+=temp+1; //查漏补缺
        }
    }

    return ans;
}

int main()
{
    cin>>T;
    init();
    rep(kase,1,T)
    {
        scanf("%lld",&n);
        printf("Case %d: %lld\n",kase,cal(n));
    }
    return 0;
}