bzoj 4918: 回文数对

传送门

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>
#include<algorithm>
#define ll long long
using namespace std;

ll n,m,T,f[][];
ll min(ll a,ll b){return a<b?a:b;}
ll work(ll n,ll a,ll b){
    ll MMH=;
    memset(f,,sizeof(f));
    f[][]=;
    for (ll i=;i<=(n>>);i++)
    for (ll qa=;qa<;qa++)
    for (ll qb=;qb<;qb++)
    for (ll ha=;ha<;ha++)
    for (ll hb=;hb<;hb++)
    if ((qa^qb)==(ha^hb)&&(i!=||qa!=qb)){
        ll QA=(a>>n-i)&,QB=(b>>n-i)&;
        ll HA=(a>>i-)&,HB=(b>>i-)&;

        for (ll k=;k<;k++){
            ll K=k;
            if (QA<qa&&!(k&)) continue;
            if (QB<qb&&!(k&)) continue;

            if (QA>qa) K|=;
            if (QB>qb) K|=;

            if (HA>ha) K|=;else if (HA<ha) K&=;
            if (HB>hb) K|=;else if (HB<hb) K&=;
            f[i][K]+=f[i-][k];
        }
    }

    if (n&){
        ll i=n+>>;
        for (ll _a=;_a<;_a++)
        for (ll _b=;_b<;_b++)
        if (n!=||_a!=_b){
            ll A=(a>>i-)&,B=(b>>i-)&;

            for (ll k=;k<;k++){
                ll K=k;
                if (A<_a&&!(k&)) continue;
                if (B<_b&&!(k&)) continue;

                if (A>_a) K|=;
                if (B>_b) K|=;
                f[i][K]+=f[i-][k];
            }
        }
    }
    ll i=n+>>;
    for (ll k=;k<;k++)
    if (((k&)||(k&))&&((k&)||(k&))) MMH+=f[i][k];
    return MMH;
}
ll Mavis(ll a,ll b){
    if (b<) return ;
    ll i,j,k,mmh=;
    for (i=,j=,k=;i;i--){
        j<<=;k<<=;j|=(a>>i)&;k|=(b>>i)&;
        if (j==k) mmh+=j*work(i,(1ll<<i)-,(1ll<<i)-)+work(i,a&((1ll<<i)-),b&((1ll<<i)-));else mmh+=k*work(i,(1ll<<i)-,(1ll<<i)-)+work(i,(1ll<<i)-,b&((1ll<<i)-));
    }
    return mmh+b+;
}
int main(){
    scanf("%lld",&T);
    while(T--){
        scanf("%lld%lld",&n,&m);
        printf("%lld\n",Mavis(m,m)-*Mavis(m,n-)+Mavis(n-,n-));
    }
}