bzoj1799[Ahoi2009]self 同类分布

[Ahoi2009]self 同类分布

Time Limit: 50 Sec  Memory Limit: 64 MB
Submit: 1558  Solved: 687
[Submit][Status][Discuss]

Description

给出a,b，求出[a,b]中各位数字之和能整除原数的数的个数。

Input

 

Output

 

Sample Input

10 19

Sample Output

3

HINT

【约束条件】1 ≤ a ≤ b ≤ 10^18

Source

Day1

 

题解：

　　数位dp，f[i][j][k][0/1]表示到第i位数位和为j，在模sum意义下的余数为k，是否卡上界的数的个数。 

#include<cstring>
#include<cmath>
#include<algorithm>
#include<iostream>
#include<cstdio>

#define N 166
#define ll long long
using namespace std;

int mod,tot,cnt;
int a[25],mark[2][25][N][N]; 
ll f[2][25][N][N];
ll dp(int p,int d,int s,int v)
{
    if (!d) return !s && !v;
    if (mark[p][d][s][v]==tot) return f[p][d][s][v];
    mark[p][d][s][v]=tot; ll t=0;
    int i,l=max(0,s-(d-1)*9),r=min((p)?9:a[d],s);
    for (i=l; i<=r; i++) t+=dp(p|(i<a[d]),d-1,s-i,(v*10+i)%mod);
    return f[p][d][s][v]=t;
}
ll solve(ll x)
{
    ll t=0;
    for (cnt=0; x; x/=10) a[++cnt]=x%10;
    for (mod=1; mod<=cnt*9; mod++)
    {
        tot++; 
        t+=dp(0,cnt,mod,0);
    }
    return t;
}
int main()
{
    ll x,y; scanf("%lld%lld",&x,&y);
    printf("%lld
",solve(y)-solve(x-1));
}