码迷,mamicode.com
首页 > 其他好文 > 详细

[HAOI2010] 计数 - 数位dp,组合数

时间:2020-02-16 16:28:24      阅读:69      评论:0      收藏:0      [点我收藏+]

标签:高精度   全排列   cal   define   康托展开   int   sig   cpp   turn   

你有一组非零数字(不一定唯一),你可以在其中插入任意个0,这样就可以产生无限个数。比如说给定{1,2},那么可以生成数字12,21,102,120,201,210,1002,1020,等等。

现在给定一个数,问在这个数之前有多少个数。(注意这个数不会有前导0).

Solution

可重复康托展开

  • 常用数位 dp 套路,枚举哪一位开始比原数小,前方唯一而后方算全排列

  • 回避高精度的全排列数计算,考虑到阶乘最多大约算到 50, 开个因子计数桶,然后在桶上操作就可以了

#include <bits/stdc++.h>
using namespace std;

#define int long long
const int N = 1005;
char str[N];
int n;

int buc[N];

void push(int x) {
    for(int i=2;i<=x;i++) {
        while(x%i==0) buc[i]++, x/=i;
    }
}

void pop(int x) {
    for(int i=2;i<=x;i++) {
        while(x%i==0) buc[i]--, x/=i;
    }
}

int calc(vector <int> v) {
    int sum=0;
    for(int i=0;i<10;i++) sum+=v[i];
    for(int i=2;i<=sum;i++) push(i);
    for(int i=0;i<10;i++) {
        for(int j=2;j<=v[i];j++) pop(j);
    }
    int ans=1;
    for(int i=2;i<=50;i++) while(buc[i]) buc[i]--, ans*=i;
    return ans;
}

signed main() {
    cin>>str+1;
    n=strlen(str+1);
    vector <int> v;
    for(int i=0;i<10;i++) v.push_back(0);
    for(int i=1;i<=n;i++) v[str[i]-'0']++;
    int ans=0;
    for(int i=1;i<=n;i++) {
        for(int j=0;j<str[i]-'0';j++) {
            if(v[j]==0) continue;
            v[j]--;
            ans+=calc(v);
            v[j]++;
        }
        v[str[i]-'0']--;
    }
    cout<<ans<<endl;
}

[HAOI2010] 计数 - 数位dp,组合数

标签:高精度   全排列   cal   define   康托展开   int   sig   cpp   turn   

原文地址:https://www.cnblogs.com/mollnn/p/12317129.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!