bzoj4918: 回文数对

Description

给定区间[L,R]，请统计有多少对整数A,B(L<=A,B<=R)满足A xor B的值在二进制表示下，去掉所有前导0后是回文串

Input

第一行包含一个正整数T(1<=T<=100)，表示测试数据的组数。

每组数据包含一行两个整数L,R(0<=L<=R<=10^12)，含义如题面所述。

Output

对于每组数据输出一行一个整数，即满足条件的整数对个数。

枚举位数，从两侧向中间数位dp

#include<cstdio>
#include<cstring>
#define F(i,n) for(int i=0;i<n;++i)
typedef long long i64;
const int M=;
int T;
i64 L,R;
i64 f[M+][][][][];
i64 min(i64 a,i64 b){return a<b?a:b;}
inline int tr(int a,int b){return a==?b:a;}
i64 cal(i64 m1,i64 m2){
    if(m1<||m2<)return ;
    i64 mn=min(m1,m2);
    i64 ans=;
    for(int p1=M;p1>=-;--p1){
        memset(f,,sizeof(f));
        f[p1+][(m1>>p1+)==(mn>>p1+)][][(m2>>p1+)==(mn>>p1+)][]=;
        f[p1+][][][][]+=mn>>(p1+);
        int l,r;
        for(l=p1,r=;l>=r;--l,++r){
            int xb1=(m1>>l&)-,xb2=(m1>>r&)-;
            int yb1=(m2>>l&)-,yb2=(m2>>r&)-;
            F(x1,)F(x2,)F(y1,)F(y2,)
            for(int z=(l==p1);z<;++z)
            F(t1,)F(t2,)if(l>r||t1==t2){
                f[l][tr(x1,t1-xb1)][tr(t2-xb2,x2)][tr(y1,(t1^z)-yb1)][tr((t2^z)-yb2,y2)]+=f[l+][x1][x2][y1][y2];
            }
        }
        F(x1,)F(x2,)F(y1,)F(y2,)if(tr(x1,x2)!=&&tr(y1,y2)!=)ans+=f[l+][x1][x2][y1][y2];
    }
    return ans;
}
int main(){
    for(scanf("%d",&T);T;--T){
        scanf("%lld%lld",&L,&R);
        printf("%lld\n",cal(R,R)+cal(L-,L-)-cal(L-,R)-cal(R,L-));
    }
    return ;
}