hdu4389 X mod f(x) 数位DP

题意很简单，就不说了。

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

x%f(x)，看起来很复杂，但是仔细一想，f(x)的变换范围只有[1,81]

那么，x%f(x) == 0，不妨变成 x%k==0 && f(x)==k ，其中k枚举。

剩下的，就简单了。数位DP，写来写去都那样。

（ 那个预处理优化，偶手残，没法写成好看的记忆化搜索 o(╯□╰)o  ）

//#pragma comment(linker, "/STACK:102400000,102400000")
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<cmath>
#include<cctype>
#include<string>
#include<algorithm>
#include<iostream>
#include<ctime>
#include<map>
#include<set>
using namespace std;
#define MP(x,y) make_pair((x),(y))
#define PB(x) push_back(x)
typedef long long LL;
//typedef unsigned __int64 ULL;
/* ****************** */
const int INF=1000111222;
const double INFF=1e200;
const double eps=1e-8;
const LL mod=1000000007;
const int NN=2510;
const int MM=100010;
/* ****************** */

int dp[11][91][91][91];
int a[15];
LL pow_10[15];
const int limit=82;

//dp[i][j][k][k1];
//表示i位数，和值为j，模上k（这个枚举），余数为k1的数的个数

void init()
{
    memset(dp,0,sizeof(dp));
    int i,j,k,k1,jj,zt;

    pow_10[0]=1;
    for(i=1;i<=10;i++)
    {
        pow_10[i]=pow_10[i-1]*10;
    }

    for(k=1;k<limit;k++)
    {
        dp[0][0][k][0]=1;
        for(i=0;i<10;i++)
            for(j=0;j<limit;j++)
                for(k1=0;k1<limit-1;k1++)
                    if(dp[i][j][k][k1])
                    {
                        for(jj=0;jj<10;jj++)
                        {
                            zt=(k1*10+jj)%k;
                            dp[i+1][j+jj][k][zt]+=dp[i][j][k][k1];
                        }
                    }
    }
}

int dfs(int pos,int sum,int chu,int yu,bool fg)
{
    if(!fg)
    {
        int fuck=pow_10[pos]*yu%chu;
        if(fuck!=0)
            fuck=chu-fuck;
        return dp[pos][sum][chu][fuck];
    }
    if(pos==0)
    {
        if(sum==0 && yu==0)
            return 1;
        return 0;
    }
    int i,ans=0;
    for(i=0;i<=a[pos];i++)
    {
        if(sum>=i)
            ans+=dfs(pos-1,sum-i,chu,(yu*10+i)%chu,i==a[pos]);
    }
    return ans;
}

int solve(int x)
{
    int i,tol=0;
    int ans=0;
    do
    {
        a[++tol]=x%10;
        x/=10;
    }while(x);
    for(i=1;i<limit;i++)
    {
        ans+=dfs(tol,i,i,0,1);
    }
    return ans;
}

int main()
{
    init();
    int cas,ee=0;
    int a,b,ans;
    scanf("%d",&cas);
    while(cas--)
    {
        scanf("%d%d",&a,&b);
        ans=solve(b)-solve(a-1);
        printf("Case %d: ",++ee);
        printf("%d\n",ans);
    }
    return 0;
}