CF908G New Year and Original Order（DP，数位 DP）

又一次降智……

（数位 DP 原来可以写这么短，学到了）

问题可以转化为求数位中 $\ge k$ 的有恰好 $j$ 位的数的个数。设为 $c_{j,k}$。

那么答案就是：（考虑把 $k$ 的贡献拆开，比如 $9$ 的贡献拆成 $1$ 的贡献的 $9$ 倍，然后分配到 $1$ 到 $9$）

$$\sum_{1\le j\le n,1\le k\le 9}c_{j,k}\underbrace{111\dots111}_{j}$$

求 $c_{j,k}$ 可以数位 DP。（此处不是记忆化搜索的形式，所以方程会有点不同）

$f[i][j][k][l]$ 表示前 $i$ 位，$\ge k$ 的有恰好 $j$ 位，$l$ 是有没有顶到上界（就是后面要不要顶着 $n$ 枚举，类似记忆化搜索中的 $limit$）。

转移较为显然，不再赘述。

时间复杂度 $O(len^2)$。（有一个 $100$ 的常数）

#include<bits/stdc++.h>
using namespace std;
const int maxn=,mod=1e9+;
#define FOR(i,a,b) for(int i=(a);i<=(b);i++)
#define ROF(i,a,b) for(int i=(a);i>=(b);i--)
#define MEM(x,v) memset(x,v,sizeof(x))
inline int read(){
    int x=,f=;char ch=getchar();
    while(ch<'' || ch>'') f|=ch=='-',ch=getchar();
    while(ch>='' && ch<='') x=x*+ch-'',ch=getchar();
    return f?-x:x;
}
int n,f[maxn][maxn][][],ans;
char s[maxn];
int calc(int x){
    int pro=,s=;
    FOR(i,,x){
        s=(s+pro)%mod;
        pro=10ll*pro%mod;
    }
    return s;
}
int main(){
    scanf("%s",s+);
    n=strlen(s+);
    FOR(i,,) f[][][i][]=;
    FOR(i,,n) FOR(j,,i) FOR(k,,) FOR(l,,) if(l<k || j){
        if(l==s[i]-'') f[i][j][k][]=(f[i][j][k][]+f[i-][j-(l>=k)][k][])%mod;
        f[i][j][k][]=(f[i][j][k][]+f[i-][j-(l>=k)][k][])%mod;
        if(l<s[i]-'') f[i][j][k][]=(f[i][j][k][]+f[i-][j-(l>=k)][k][])%mod;
    }
    FOR(j,,n) FOR(k,,) ans=(ans+1ll*calc(j)*(f[n][j][k][]+f[n][j][k][])%mod)%mod;
    printf("%d\n",ans);
}