Entirely Unsorted Sequences (BAPC2018 I)

Entirely Unsorted Sequences (BAPC2018 I)

题意：问含有$n(n<5\ast 10^3)$个元素的序列的全错（含重复元素）排列方案数量。

定义：全错排列为对于每个元素，比他大的都不在他的左边，比他小的都不在他的右边。

题解：

先考虑不含有重复元素的情况，令$dp[n]$为含有$n$个元素的序列的全错排列方案数,
那么有递推关系：$dp[n]=n!-\sum _{i=1}^{n}dp[i-1]\ast (n-i)!$
反过来想，考虑($n$个元素的序列的全错排列方案数）为（$n$个元素的全排列方案数)减去(至少有一个元素被摆放正确的方案数)，问题就变成了要求至少有一个元素被摆放正确的方案数，我们枚举排序后的元素，令第$i$个元素是排序后第一个正确被摆放的元素，那么第$i$个元素之前则必须全为错排，即$dp[i-1]$，且第$i$个元素之后的元素可以被任意摆放，即$(n-i)!$，最后相乘即可。

推广：

现在考虑含有重复元素的情况，我没只要把上式中的（$n!$）换成（$n$个元素的全排列的方案数）即可。
公式为 $cnt(n)=\frac{n!}{c_1\ast c_2\ast \cdots \ast c_k}$，其中$k$为不同元素个数，$c_i$为每个元素的数量,
例如序列含有10个元素的序列$[1223444445]$，则$cnt(n)=\frac{10!}{1!\ast 2!\ast 1!\ast 5!\ast 1!}$（简单理解为全排列除以每个元素内部重复排列的方案数）
同时还需要把$(n-i)!$换成$i+1$到$n$这个区间内元素的全排列方案数。
这个推广的整个过程可以使用动态规划完成，即在$O(n^{2}log(n))$的时间内，求出$cnt[i][j]$，即区间$i$到$j$内元素的全排列方案数。

代码：

typedef long long ll;
const int MAXN = 5e3 + 10;
ll fac[MAXN] = { 1,1 };
ll inv[MAXN] = { 1,1 };
ll a[MAXN], dp[MAXN];

ll num[MAXN][MAXN];
int main() {
	int n;
	cin >> n;
	for (int i = 2; i <= n; i++) {
		fac[i] = fac[i - 1] * i % MOD;
		inv[i] = inv[MOD % i] * (MOD - MOD / i) % MOD;
	}
	
	for (int i = 1; i <= n; i++) cin >> a[i];
	sort(a, a + n + 1);
	for (int i = 1; i <= n; i++) {
		map<int, int> ma;
		num[i][i - 1] = 1;
		for (int j = i; j <= n; j++) {
			num[i][j] = num[i][j - 1] * (j - i + 1) % MOD * inv[++ma[a[j]]] % MOD;
		}
	}

	dp[0] = 1, num[n + 1][n] = 1;
	for (int i = 1; i <= n; i++) {
		ll sum = 0;
		for (int j = 1; j <= i; j++) {
			sum += dp[j - 1] * num[j+1][i] % MOD;
			sum %= MOD;
		}
		dp[i] = (num[1][i] - sum + MOD ) % MOD;
	}
	
	cout << dp[n] << endl;
}