HDU 4389 X mod f(x) 數位DP

題意：求一個範圍內 能被自身每一位數字和整除的數的個數。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <stack>
using namespace std;
#define FOR(i,l,r) for(int i=(l);i<=(r);++i)
#define REP(i,n) for(int i=0;i<(n);++i)
#define DSC(i,r,l) for(int i=(r);i>=(l);--i)
#define N 83
int dp[10][N][N][N];
int f[15];
void init()
{   f[0]=1;
    f[1]=1;
    for(int i=2;i<=10;i++) f[i]=f[i-1]*10;
    memset(dp,0,sizeof(dp));
    FOR(i,0,81) dp[0][0][i][0]=1;
    FOR(i,0,8)
        FOR(j,0,i*9)
            FOR(k,1,81)
                REP(q,k)
                {
                    if(!dp[i][j][k][q]) continue;
                    for(int p=0;p<10 && j+p<N;++p)
                        dp[i+1][j+p][k][ (q*10+p)%k ]+=dp[i][j][k][q];//同餘公式..?還沒弄懂
                }
}

int num[15];
int solve(int x)
{
    if(x==0) return 0;
    int ans=0,sum=0,temp=0,len=0;
    if(x==1e9) ans++,x--;
    if(x==1e9-1) ans++,x--;
    x++;
    while(x)
    {
        num[++len]=x%10;
        x/=10;
    }

    FOR(mod,1,9*len) //先枚舉mod
    {
        int sum=mod,temp=0;
        for(int i=len;i>=1;i--)
        {
            for(int j=0;j<num[i] && j<=sum;j++)
            {
                for(int res=0;res<mod;res++)
                {
                    if((temp+j*f[i]+res)%mod==0) //之前確定的位數加上枚舉的這位的數字的和,找一個res對mod取餘得0
                    ans+=dp[i-1][sum-j][mod][res];
                }
            }
            sum-=num[i];
            temp+=num[i]*f[i];
        }
    }
    return ans;
}

int main()
{
    int cas,cas1=1;
    int x,y;
    init();
    cin>>cas;
    while(cas--)
    {
        scanf("%d%d",&x,&y);
        printf("Case %d: %d\n",cas1++,solve(y)-solve(x-1));
    }
    return 0;
}