薇小薇

”喵喵喵?“

[BZOJ2425][HAOI2010]计数(数位+组合数学)

题目:

我是超链接

题解:

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

我们先把所有有的数字记录下来,删掉零相当于把0移到最前面,所以和其他数字没什么区别
当各个数字个数分别为a0,a1...a9时,自由排列能产生的数的个数为(a0+a1+...+a9)!a0!a1!...a9!
然后我们就可以开始数位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);
}
阅读更多
版权声明:口胡文章,请谨慎转载╮(╯▽╰)╭ https://blog.csdn.net/Blue_CuSO4/article/details/79971061
文章标签: 数位 组合数学
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

不良信息举报

[BZOJ2425][HAOI2010]计数(数位+组合数学)

最多只允许输入30个字

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭