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

题目：

我是超链接

题解：

数位dp的思想应该是把下界的都处理好了，然后做紧贴上界的计数

我们先把所有有的数字记录下来，删掉零相当于把0移到最前面，所以和其他数字没什么区别
当各个数字个数分别为$a0,a1...a9$时，自由排列能产生的数的个数为${(a_0+a_1+...+a_9)!\over a_0!a_1!...a_9!}$
然后我们就可以开始数位DP了
假设当前位是i，那么0~i-1这些数放在第一位时，后面就可以自由排列了，方案数就可以直接算出，所以只需要接着令当前位也为i，然后继续搞下一位就可以了
这里的阶乘最大的约是50!，肯定存不下，我们可以分解质因数来存储。

代码：

#include <cstdio>
#include <cstring>
#define LL long long
using namespace std;
LL ans;int have[15],A[15],a[55];char st[55];
int pri[15]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47};
void fj(int a,int v)
{
    int lj=0;
    while (a!=1)
    {
        while (a%pri[lj]==0)
        {
            A[lj]+=v;
            a/=pri[lj];
        }
        lj++;
    }
}
void jc(LL a,int v){for (int i=2;i<=a;i++) fj(i,v);}
LL qurry()
{
    LL r=1;
    for (int i=0;i<15;i++) 
      for (int j=1;j<=A[i];j++) r*=(LL)pri[i];
    return r;
}
int main()
{
    scanf("%s",st+1);
    int n=strlen(st+1);
    for (int i=1;i<=n;i++) a[i]=st[i]-'0',have[a[i]]++;
    for (int i=1;i<=n;i++)
    {
        for (int j=0;j<a[i];j++)
          if (have[j])
          {
            have[j]--;
            memset(A,0,sizeof(A));
            jc(n-i,1);
            for (int k=0;k<=9;k++) jc(have[k],-1);
            ans+=qurry();
            have[j]++;
          }
        have[a[i]]--;
    }
    printf("%lld",ans);
}