JZOJ 排列统计

排列统计

\(Description\)

对于给定的一个长度为n的序列{B[n]},问有多少个序列{A[n]}对于所有的i满足:A[1]~A[i]这i个数字中有恰好B[i]个数字小等于i。其中{A[n]}为1~n的一个排列,即1~n这n个数字在序列A[I]中恰好出现一次。 数据保证了至少有一个排列满足B序列。

\(Input\)

输入的第1行为一个正整数N,表示了序列的长度。
第2行包含N个非负整数,描述了序列{B[i]}。

\(Output\)

输出仅包括一个非负整数,即满足的{A[i]}序列个数。

\(Sample Input\)

3
0 1 3

\(Sample Output\)

3

\(Hint\)

【样例说明】
  对于A序列为1~3的全排列分别对应的B序列如下(冒号左边为A序列,冒号右边为对应B的序列)
  1 2 3:1 2 3
  1 3 2:1 1 3
  2 1 3:0 2 3
  2 3 1:0 1 3
  3 1 2:0 1 3
  3 2 1:0 1 3
  所以有3个满足的A序列。
【数据说明】
  对于20%的数据,有N≤8;
  对于30%的数据,有N≤11且答案不大于20000;
  对于50%的数据,有N≤100;
  对于100%的数据,有N≤2000。

解题思路

其实很容易发现,当且仅当 \(B_i - B_{i-1} = {0 , 1 , 2}\) 时才有解
因为考虑到 \(B_i\)\(B_{i-1}\) 的关系,前者是后者产生的数列再加入一个新的数
为了方便表达,记产生的数列为 \(a_{1..n}\)
那么新加入一个数,记为x , 它至多可以给 \(B_i\) 贡献 1(\(x <= i\)),此时可贡献的数从上限 \(i-1\) 提升到 \(i\) ,如果之前有一个数恰好等于 i,那么它又可以给 \(B_i\) 贡献 1
然后就找不到可贡献的了,即 \(B_i + 2 \geq B_{i-1}\) 时才可能构造数列

现在考虑如何求得答案

\(f_i\) 表示前 \(1..i\) 满足条件的数列数目,\(k = B_{i-1} - B_i\)

\[f_i = \left \{ \begin{aligned} f_{i-1} & & (k=0) \\ f_{i-1}[(i-B_{i-1})+(i-1-B_{i-1})] & & (k=1) \\ f_{i-1}(i-B_{i-1}-1)^2 & & (k=2) \end{aligned} \right. \]

要用高精度!!!

代码参考

#include<cstdio>
using namespace std;
typedef long long LL;

int b[2005] , n , ans[10005];

inline void mul(int x)
{
	int g = 0;
	for(register int i = 1; i <= ans[0]; i++)
	{
		ans[i] = ans[i] * x + g;
		g = ans[i] / 10;
		ans[i] = ans[i] % 10;
	}
	while(g)
	{
		ans[++ans[0]] = g;
		g = ans[ans[0]] / 10;
		ans[ans[0]] = ans[ans[0]] % 10;
	}
}

int main()
{
	scanf("%d" , &n);
 	for(register int i = 1; i <= n; i++) scanf("%d" , &b[i]);
 	ans[0] = ans[1] = 1;
 	for(register int i = 1; i <= n; i++)
 	{
 		if (b[i] - b[i-1] == 1) mul(i - b[i-1] + i - 1 - b[i - 1]);
 		else if (b[i] - b[i-1] == 2) mul((i - b[i - 1] - 1) * (i - b[i - 1] - 1));
	}
	for(register int i = ans[0]; i; i--) printf("%d" , ans[i]);
}
原文地址:https://www.cnblogs.com/leiyuanze/p/13306024.html