The 2021 ICPC Asia Jinan Regional Contest - C Optimal Strategy - 2021ICPC济南站C题 组合数学

题目大意

Ena和Mizuki正在玩一个游戏。

他们面前有n个项目，编号从1到n。第i个项目的价值是

 .Ena和Mizuki轮流操作，而Ena先操作。在一次移动中，玩家选择一个尚未被拿走的物品并将其带走。当所有物品都被拿走时，游戏结束。任何一方的目标都是使他们拿走的物品的价值之和最大。

鉴于双方都以最佳方式行动，有多少个可能的游戏过程？由于这个数字可能太大，你应该输出998244353的模数。

如果存在某个整数i(1≤i≤n)，使第i次移动中拿走的物品的指数不同，则认为两个过程不同。

输入描述

The first line contains an integer n (1≤n≤106).

The second line contains n integers a1​,a2​,…,an​ (1≤ai​≤n).

输出描述

Output the answer.

Sample Input 1

3
1 2 2

Sample Output 1

4

Sample Input 2

6
1 3 2 2 3 1

Sample Output 2

120

Sample Input 3

12
1 1 4 5 1 4 1 9 1 9 8 10

Sample Output 3

28800

思路

在比赛过程中，我们队用了不同于北京大学题解的方法解开这道题。由题意可知，我们在操作的时候可以取走最大的奇数个元素，也可以取走比这个数大的任意一个数。我们可以采取一个贪心策略，我们先取走大的元素，判断接下来所有元素的放置位置，然后我们再取小的元素。很显然，可以放置的位置一个有cnt[i]/2+1个（cnt[i]为我们当前枚举的元素的个数），使用挡板法，可知接下来所有元素的插空放置方案一共有C(cnt[i]/2+1+len,cnt[n]/2+1)(len是剩余元素的个数)。这样我们就可以确定了已经放置的元素和剩余元素的相对位置。放置完当前值最大的所有元素后，我们再取出在剩余元素中值最大的所有元素，用同样的办法插空，从而计算出最后结果。

可能我上面描述不是特别直观，下面我们直接上图来表示样例2是如何计算出来的。

 

 如上图所示，在样例2中，我们先放置2个值为3的元素，接下来放置剩余元素，很显然，剩余元素只能排在方框的位置，如果把剩余元素放置在两个3中间是不可能的，因为如果先手取了1个3，然后3就变成奇数个了，如果后手取了一个小于3的数，就导致了后手少得到了一个3，等到的分数就不可能为最大值。

如上图所示，剩余元素有4个，放在两个空位中，因此我们可以用挡板法计算如何放置这4个元素。接下来，我们只要确定剩下没有讨论的4个元素的相对位置即可。 

最后，在计算结果时，排列当前枚举的元素时，排列的方案有cnt[i]!种方案，不要忘记乘回来了，否则会导致最后计算结果偏小。

下面直接放上AC代码

#include<bits/stdc++.h>

#define ll long long
#define int long long 

using namespace std;

const int N = 5e6 + 5;
const int M = 2e3 + 5;
const ll mod = 998244353;

void io() {
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
}

int sum = 0;
int n;
int b[N];
int len = 0;
int a[N];
int cnt[N];
int fac[N];
int inv[N];

int qpow(int a, int b) {
	if (b == 0)return 1;
	if (b == 1)return a;
	if (b & 1)return (a * qpow(a * a % mod, b >> 1)) % mod;
	return qpow(a * a % mod, b >> 1) % mod;
}

void init() {
	fac[0] = 1;
	for (int i = 1; i < N; i++) {
		fac[i] = fac[i - 1] * (i) % mod;
	}
	inv[N - 1] = qpow(fac[N - 1], mod - 2);
	int i = N - 2;
	while (i) {
		inv[i--] = inv[i + 1] * (i + 1) % mod;
	}
	inv[0]=1;
}

int C(int a, int b)
{
	return (fac[a] % mod) * (inv[b] * inv[a - b] % mod) % mod;
}
signed main() {
    io();
	init();
	cin >> n;
	for (int i = 1; i <= n; i++) {
		cin >> b[i];
	}
	sort(b + 1, b + n + 1, [](int x, int y) {return x > y; });
	len=0;
	for (int i = 1; i <= n; i++) {
		if (b[i] != b[i - 1]) {
			cnt[++len] = 1;
		}
		else {
			cnt[len]++;
		}
	}
	int ans = 1;int Len = n;
	for (int i = 1; i < len; i++) {
		int group = cnt[i] >> 1;
		Len -= cnt[i];
		ans *= C( Len + group , group);
		ans %= mod;
	}
	for (int i = 1; i <= len; i++) {
		ans *= fac[cnt[i]];
		ans %= mod;
	}
	cout << ans << endl;
	return 0;
}