結果

問題 No.1193 Penguin Sequence
ユーザー ianCKianCK
提出日時 2020-09-25 20:21:08
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 89 ms / 2,000 ms
コード長 2,300 bytes
コンパイル時間 490 ms
コンパイル使用メモリ 53,032 KB
実行使用メモリ 10,504 KB
最終ジャッジ日時 2023-09-10 14:44:37
合計ジャッジ時間 4,882 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge15
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 55 ms
10,464 KB
testcase_01 AC 89 ms
10,428 KB
testcase_02 AC 88 ms
10,444 KB
testcase_03 AC 89 ms
10,456 KB
testcase_04 AC 89 ms
10,448 KB
testcase_05 AC 89 ms
10,460 KB
testcase_06 AC 89 ms
10,400 KB
testcase_07 AC 88 ms
10,484 KB
testcase_08 AC 89 ms
10,396 KB
testcase_09 AC 89 ms
10,412 KB
testcase_10 AC 88 ms
10,460 KB
testcase_11 AC 51 ms
9,880 KB
testcase_12 AC 51 ms
9,984 KB
testcase_13 AC 76 ms
10,388 KB
testcase_14 AC 72 ms
10,392 KB
testcase_15 AC 83 ms
10,388 KB
testcase_16 AC 46 ms
9,456 KB
testcase_17 AC 1 ms
4,376 KB
testcase_18 AC 12 ms
9,024 KB
testcase_19 AC 88 ms
10,388 KB
testcase_20 AC 69 ms
10,212 KB
testcase_21 AC 56 ms
9,952 KB
testcase_22 AC 13 ms
9,036 KB
testcase_23 AC 47 ms
9,696 KB
testcase_24 AC 42 ms
9,616 KB
testcase_25 AC 21 ms
9,208 KB
testcase_26 AC 9 ms
8,924 KB
testcase_27 AC 75 ms
10,448 KB
testcase_28 AC 52 ms
9,940 KB
testcase_29 AC 76 ms
10,468 KB
testcase_30 AC 28 ms
9,388 KB
testcase_31 AC 24 ms
9,268 KB
testcase_32 AC 60 ms
10,212 KB
testcase_33 AC 43 ms
9,728 KB
testcase_34 AC 35 ms
9,460 KB
testcase_35 AC 43 ms
9,640 KB
testcase_36 AC 34 ms
9,468 KB
testcase_37 AC 71 ms
10,504 KB
testcase_38 AC 5 ms
8,756 KB
testcase_39 AC 5 ms
8,784 KB
testcase_40 AC 5 ms
8,796 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>
#include <vector>
#include <string.h>
#include <algorithm>
using namespace std;
#define PB push_back

typedef long long int ll;
constexpr int kN = int(2E5 + 10), kMod = 998244353;

struct BIT {
	int val[kN];
	void init() {memset(val, 0, sizeof(val));}
	void add(int pos) {
		while (pos < kN) {
			val[pos]++;
			pos += pos & -pos;
		}
		return ;
	}
	int ask(int pos) {
		int ans = 0;
		while (pos) {
			ans += val[pos];
			pos ^= pos & -pos;
		}
		return ans;
	}
};

ll Pow(ll a, ll b) {
	ll ans = 1;
	while (b) {
		if (b & 1) ans = ans * a % kMod;
		a = a * a % kMod;
		b >>= 1;
	}
	return ans;
}

ll Rev(ll n) {return Pow(n, kMod - 2);}

ll f[kN], inf[kN];

void pre() {
	f[0] = f[1] = inf[0] = inf[1] = 1;
	for (int i = 2; i < kN; i++) f[i] = f[i - 1] * i % kMod;
	inf[kN - 1] = Rev(f[kN - 1]);
	for (int i = kN - 1; i > 2; i--) inf[i - 1] = inf[i] * i % kMod;
	return ;
}

ll C(int n, int m) {return f[n] * inf[m] % kMod * inf[n - m] % kMod;}

int a[kN], rnk[kN]; 
ll same[kN];
BIT bit;

int main() {
	int n;
	ll ans = 0, tot = 1, tmp = 0, cnt = 0, ccnt = 0;
	vector<int> v;
	scanf("%d", &n);
	if (n == 1) {
		printf("0\n");
		return 0;
	}
	for (int i = 1; i <= n; i++) scanf("%d", &a[i]);
	pre();
	for (int i = 1; i <= n; i++) tot = tot * C(n, i) % kMod;
	// inside
	bit.init();
	for (int i = 1; i <= n; i++) v.PB(a[i]);
	sort(v.begin(), v.end());
	v.resize(unique(v.begin(), v.end()) - v.begin());
	for (int i = 1; i <= n; i++) rnk[i] = lower_bound(v.begin(), v.end(), a[i]) - v.begin() + 1;
	for (int i = 1; i <= n; i++) {
		cnt += (i - 1) - bit.ask(rnk[i]);
		bit.add(rnk[i]);
	}
	for (int i = 2; i <= n; i++) tmp += ll(i) * (i - 1) % kMod;
	tmp %= kMod;
	cnt %= kMod;
	ans = tmp * tot % kMod * Rev(n) % kMod * Rev(n - 1) % kMod * cnt % kMod;
	// between
	for (int i = 1; i <= n; i++) same[rnk[i]]++;
	for (int i = 1; i <= n; i++) if (same[rnk[i]]) {
		ccnt += (n - same[rnk[i]]) * same[rnk[i]] % kMod;
		same[rnk[i]] = 0;
	}
	ccnt %= kMod;
	ccnt = ccnt * Rev(2) % kMod;
	tmp = ll(n) * (n + 1) % kMod * Rev(2) % kMod;
	tmp = tmp * tmp % kMod;
	for (int i = 1; i <= n; i++) tmp -= ll(i) * i % kMod;
	tmp %= kMod;
	if (tmp < 0) tmp += kMod;
	tmp = tmp * Rev(2) % kMod;
	ans += tmp * tot % kMod * ccnt % kMod * Rev(n) % kMod * Rev(n) % kMod;

	printf("%lld\n", ans % kMod);
}
0