[HAOI2010]计数，组合数+数位Dp

最新推荐文章于 2019-08-07 12:02:48 发布

Deep_Kevin

最新推荐文章于 2019-08-07 12:02:48 发布

阅读量376

点赞数

分类专栏：组合数

本文链接：https://blog.csdn.net/Deep_Kevin/article/details/83003066

版权

组合数专栏收录该内容

5 篇文章 0 订阅

订阅专栏

正题

[HAOI2010]计数

“给您一个字符串，要你求这个字符串在他的全排列中排第几”，题面很明显可以转化为这一个问题。

那么对于每一位，我们就可以计算当前这位小于它的所有全排列。

具体做法就是先把枚举到的j先删除，然后拿剩下的数进行全排列。具体的全排列公式是

$ans=\frac{n!}{\prod_{i=0}^na[i]!}$

但是很明显会爆long long，~~所以不想打高精度的情况下~~就有了第二条公式:

$C(m,a[0])*C(m-a[0],a[1])*...*C(m-\sum_{i=0}^8a[i],a[9])$

a[i]表示i有多少个。

考虑先把0放进去，那么就相当于在m个位置中找a[0]个位置。

在考虑吧1放进去，那么就相当于在m-a[0]个韦志中找a[1]个位置

那么就每一次算一下就可以了。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
using namespace std;

int s[55],a[10];
long long ans=0;
int n=0,t;

long long C(int x,int y){
	long long tot=1;
	for(int i=x;i>=x-y+1;i--) tot*=i;
	for(int i=y;i>=2;i--) tot/=i;
	return tot;
}

long long calc(){
	int m=t;
	long long tot=1;
	for(int i=0;i<=9;i++) if(a[i]) tot*=C(m,a[i]),m-=a[i];
	return tot;
}

int main(){
	char ch;
	while(scanf("%c",&ch)!=EOF) if(ch>='0' && ch<='9')s[++n]=ch-'0',a[s[n]]++;
	t=n;
	for(int i=1;i<=n;i++){
		t--;
		for(int j=0;j<s[i];j++)
			if(a[j]!=0) {a[j]--;ans+=calc();a[j]++;}
		a[s[i]]--;
	}
	printf("%lld",ans);
}