[BZOJ2425][HAOI2010]计数（组合数学）

首先，统计出数字串中$0,1,...,9$出现的次数$cnt[0],cnt[1],...,cnt[9]$。
先考虑用其中的所有非零数字能构成多少个允许含前导$0$的$len$位数。
可以想到，现在$len$位中选出$cnt[1]$位设为$1$，然后在剩下的$len-cnt[1]$位中选出$cnt[2]$位设为$2$，再在剩下的$len-cnt[1]-cnt[2]$位中选出$cnt[3]$位设为$3$…，即结果为：
$\prod_{i=1}^9C_{len-\sum_{j=1}^{i-1}cnt[j]}^{cnt[i]}$
而对于不含前导$0$的，就只要枚举最高位，就可以用上面的式子了（有些小细节，具体见代码实现）。
代码实现如下（$lowe(len,0/1)$表示构成$len$位数，第二个参数为$0$表示可以含前导$0$，为$1$表示不能含前导$0$）：

ll lowe(int len, bool first) {
    int i, j; ll ans = 1; if (!first) {
        for (i = 1; i <= 9; len -= cnt[i++])
            ans *= C[len][cnt[i]];
        return ans;
    }
    ans = 0; for (i = 1; i <= 9; i++) {
        cnt[i]--; ans += lowe(len - 1, 0);
        cnt[i]++;
    }
    return ans;
}

继续考虑求位数相等且比原数字串小的数字串个数。
考虑最高位的取值，分两部分计算：
1、最高位的取值比原数字串的最高位小。这样后面位的值不管取多少，都是不会大于等于原数字串的。这时候可以用上面的方法求得。
2、最高位的取值等于原数字串的最高位。这时候就要求去掉最高位之后比原数字串小。但可以发现，这是和原问题相同的一个子问题，可以递归处理（在标程中使用的是非递归）。
注意判断前导$0$。
组合数可以使用公式$C_n^m=C_{n-1}^m+C_{n-1}^{m-1}$预处理。
代码：

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
const int N = 1005, R = 12;
int n, num[N], tot, cnt[R];
char s[N]; ll C[N][N], Ans;
void init() {
    int i, j; for (i = 0; i <= 1000; i++) C[i][0] = 1;
    for (i = 1; i <= 1000; i++) for (j = 1; j <= i; j++)
        C[i][j] = C[i - 1][j] + C[i - 1][j - 1];
}
ll lowe(int len, bool first) {
    int i, j; ll ans = 1; if (!first) {
        for (i = 1; i <= 9; len -= cnt[i++])
            ans *= C[len][cnt[i]];
        return ans;
    }
    ans = 0; for (i = 1; i <= 9; i++) {
        cnt[i]--; ans += lowe(len - 1, 0);
        cnt[i]++;
    }
    return ans;
}
void solve() {
    int i, j; for (i = 1; i <= n; i++) {
        for (j = (i == 1); j < num[i]; j++) {
            if (!cnt[j]) continue;
            cnt[j]--; Ans += lowe(n - i, 0);
            cnt[j]++;
        }
        cnt[num[i]]--;
    }
}
int main() {
    int i; scanf("%s", s + 1); n = strlen(s + 1);
    for (i = 1; i <= n; i++) {
        cnt[num[i] = s[i] - '0']++;
        if (num[i]) tot++;
    }
    init(); for (i = tot; i < n; i++) Ans += lowe(i, 1);
    cout << (solve(), Ans) << endl;
    return 0;
}