Palindrome Function HDU - 6156 数位DP

#include<cstdio> 
#include<iostream>
int dp[39][39][37]={0};
int an[39][37]={0};
using namespace std;
int work(int x,int k)
{
	if(x==0) return 0;
	int num[40]={0};
	int cnt=0;
	while(x)
	{
		++cnt;
		num[cnt] =x%k;
		x /=k;
	}
	int ans=0;
	ans +=an[cnt-1][k];
	for(int i=1;i<num[cnt];i++) ans +=dp[cnt][i][k]; 
	//累加 x.....x x<a 
	int len=cnt,l=1,flag=1;
	while((cnt -=2)>0)
	{
		for(int i=0;i<num[len-l];i++) ans +=dp[cnt][i][k];
		l++;
    }
    /* abcdefg 
	   累加 axx....axx  x<b
	         abx....bax x<c
			   abcxcba  x<e 答案 显然成立
			   x=e 可能成立 
	*/
    int i=len/2-1;
    while(num[len-i]==num[1+i]&&i>=0) i--;
    flag = i==-1 ?1: num[len-i]>num[1+i] ? 0:1;
     /*
	  验证 x=e 的情况
	*/ 
    return ans+flag;
}
int main()
{
	
	for(int k=2;k<=36;k++)
	 {
	   for(int i=1;i<=2;i++)
	    for(int j=0;j<k;j++)
        {
	     dp[i][j][k]=1;
	     if(j>0) an[i][k]+=dp[i][j][k]; 
        }
	    an[2][k] +=an[1][k];
     }
    for(int q=2;q<=36;q++)
   {
   	int m=2e9+7,cc=0;
   	while(m) cc++,m /=q;
	for(int i=3;i<=cc;i++)
	{
		an[i][q] +=an[i-1][q];
		 for(int j=0;j<q;j++)
		 {
		  for(int k=0;k<q;k++)
		   {
		   	 dp[i][j][q] +=dp[i-2][k][q];
		   }
		   an[i][q] +=j==0?0:dp[i][j][q];
	     }
	}
   }
	int l,r,l1,r1,t;
	scanf("%d",&t);
	for(int j=1;j<=t;j++)
	{
	scanf("%d%d%d%d",&l,&r,&l1,&r1);
		long long ans=0;
		for(int i=l1;i<=r1;i++)
		{
			long long d= (work(r,i) -work(l-1,i));
			d *=(i-1);
			ans +=d;
			ans +=r-l+1;
		}
		printf("Case #%d: %lld\n",j,ans);
    }
	return 0;
}