【数位DP】HDU 6156 Palindrome Function

http://acm.hdu.edu.cn/showproblem.php?pid=6156

【AC】

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod=1e9+7;
int num[50];
int L,R,l,r;
ll query(int x,int k)
{
    if(x==0) return 0;
    int cnt=0;
    int cp=x;
    while(cp)
    {
        num[cnt++]=cp%k;
        cp/=k;
    }
    int base=k,sum=1,tot=0;
    for(int i=1;i<cnt;i++)
    {
        tot+=base-sum;
        if(i%2==0)
        {
            sum=base;
            base*=k;
        }
    }
    int p=0;
    for(int i=cnt-1;i>=cnt/2;i--)
    {
        p=p*k+num[i];
    }
    tot+=p-sum;
    int o=(cnt%2==0?cnt/2:cnt/2+1);
    for(int i=o;i<cnt;i++)
    {
        p=p*k+num[i];
    }
    if(p<=x) tot+=1;
    return 1ll*tot*(k-1)+(ll)x;
} 
int main()
{
    int T;
    scanf("%d",&T);
    int cas=0;
    while(T--)
    {
        scanf("%d%d%d%d",&L,&R,&l,&r);
        ll ans=0;
        for(int i=l;i<=r;i++)
        {
            ans+=query(R,i)-query(L-1,i);
        }
        printf("Case #%d: %I64d
",++cas,ans);
    }
    return 0;
}