poj 3252 Round Numbers(数位dp+组合计数)

题意：询问区间[L,R]内中的二进制数0>=1的数的个数。

思路：网上很多都直接推导一下就算出来了。。。。我只会用数位dp搞一下，再计数Orz……觉得自己数学思维实在太差了、、、、dp[i][j][k]表示当前在第i位，前面的有j个0，有k个1的合法的数有多少个。当前面有1的时候，可以算一下至少需要多少个0才能满足要求，比如现在有n位，至少需要m个0，那么结果就为 Σ{i=m to n}（C[n][i]）。

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<set>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
const int maxn=35;
ll C[maxn][maxn];
void Init()
{
    memset(C,0,sizeof(C));
    C[0][0]=1;
    for(int i=1;i<maxn;++i)
    {
        C[i][0]=1;
        for(int j=1;j<=i;++j)
            C[i][j]=C[i-1][j-1]+C[i-1][j];
    }
}
int num[maxn];
ll dp[maxn][maxn][maxn];
ll f(int pos,int zc,int oc,bool limit)
{
    if(pos==-1) return zc>=oc;
    if(!limit&&dp[pos][zc][oc]!=-1) return dp[pos][zc][oc];
    ll ans=0;
    if(oc&&!limit)
    {
        int tmp=zc-oc;
        tmp=(pos+1)-tmp;
        int x=tmp/2+((tmp&1)?1:0);
        if(tmp<0) x=0;
        for(int i=x;i<=pos+1;++i)
            ans+=C[pos+1][i];
        return dp[pos][zc][oc]=ans;
    }
    int last=limit?num[pos]:1;
    for(int i=0;i<=last;++i)
    {
        ans+=f(pos-1,zc+(oc?(i^1):0),oc+i,limit&&i==last);
    }
    if(!limit) dp[pos][zc][oc]=ans;
    return ans;
}
ll cal(ll x)
{
    memset(dp,0xff,sizeof(dp));
    int len=0;
    while(x)
    {
        num[len++]=(x&1);
        x>>=1;
    }
    return f(len-1,0,0,true);
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    Init();
    ll S,T;
    while(cin>>S>>T)
    {
        cout<<cal(T)-cal(S-1)<<endl;
    }
    return 0;
}