bzoj4918: 回文数对

#include<cstdio>
#include<cstring>
#define F(i,n) for(int i=0;i<n;++i)
typedef long long i64;
const int M=40;
int T;
i64 L,R;
i64 f[M+5][3][3][3][3];
i64 min(i64 a,i64 b){return a<b?a:b;}
inline int tr(int a,int b){return a==1?b:a;}
i64 cal(i64 m1,i64 m2){
    if(m1<0||m2<0)return 0;
    i64 mn=min(m1,m2);
    i64 ans=0;
    for(int p1=M;p1>=-1;--p1){
        memset(f,0,sizeof(f));
        f[p1+1][(m1>>p1+1)==(mn>>p1+1)][1][(m2>>p1+1)==(mn>>p1+1)][1]=1;
        f[p1+1][0][1][0][1]+=mn>>(p1+1);
        int l,r;
        for(l=p1,r=0;l>=r;--l,++r){
            int xb1=(m1>>l&1)-1,xb2=(m1>>r&1)-1;
            int yb1=(m2>>l&1)-1,yb2=(m2>>r&1)-1;
            F(x1,2)F(x2,3)F(y1,2)F(y2,3)
            for(int z=(l==p1);z<2;++z)
            F(t1,2)F(t2,2)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+1][x1][x2][y1][y2];
            }
        }
        F(x1,2)F(x2,3)F(y1,2)F(y2,3)if(tr(x1,x2)!=2&&tr(y1,y2)!=2)ans+=f[l+1][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-1,L-1)-cal(L-1,R)-cal(R,L-1));
    }
    return 0;
}