BZOJ2425: [HAOI2010]计数

题目大意：给你一个最多50位的数，求在所有通过由他去掉0和重新排列数字顺序能产生的数中，比他小的有多少个

首先位数少的肯定有去掉0的操作，我们可以默认是把0挪到了前面
问题就变成了把这些数字重新排列顺序能产生多少个比他小的数
首先需要知道当各个数字个数分别为$a_{0},a_{1}...a_{9}$时，自由排列能产生的数的个数为$\frac{(a_{0}+a_{1}+...+a_{9})!}{a_{0}!a_{1}!...a_{9}!}$
然后我们就可以开始数位DP了
假设当前位是i，那么0~i-1这些数放在第一位时，后面就可以自由排列了，方案数就可以直接算出，所以只需要接着令当前位也为i，然后继续搞下一位就可以了

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
char s[101];
long long n,a[51];
long long re[10];
long long A[16],P[16]={0,2,3,5,7,11,13,17,19,23,29,31,37,41,43,47};
long long ANS;
void change(long long x,long long k)
{
    long long j=0;
    while(x!=1)
    {
        j++;
        while(x%P[j]==0)
        {
            x/=P[j];
            A[j]+=k;
        }
    }
}
void jc(long long x,long long k)
{
    long long i;
    for(i=2;i<=x;i++)
    change(i,k);
}
long long cal()
{
    long long i,j;
    long long ret=1;
    for(i=1;i<=15;i++)
    for(j=1;j<=A[i];j++)
    ret=ret*P[i];
    return ret;
}
int main()
{
    scanf("%s",s+1);
    n=strlen(s+1);
    long long i,j,k;
    for(i=1;i<=n;i++)
    {
        a[i]=s[i]-48;
        re[a[i]]++;
    }
    for(i=1;i<=n;i++)
    {
        for(j=0;j<a[i];j++)
        if(re[j])
        {
            re[j]--;
            memset(A,0,sizeof(A));
            jc(n-i,1);
            for(k=0;k<=9;k++)
            jc(re[k],-1);
            ANS+=cal();
            re[j]++;
        }
        re[a[i]]--;
    }
    printf("%lld",ANS);
}