結果
問題 | No.1193 Penguin Sequence |
ユーザー | SSRS |
提出日時 | 2020-08-22 16:47:18 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,427 bytes |
コンパイル時間 | 2,233 ms |
コンパイル使用メモリ | 179,256 KB |
実行使用メモリ | 18,276 KB |
最終ジャッジ日時 | 2024-04-23 10:54:05 |
合計ジャッジ時間 | 12,456 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | WA | - |
testcase_01 | WA | - |
testcase_02 | WA | - |
testcase_03 | WA | - |
testcase_04 | WA | - |
testcase_05 | WA | - |
testcase_06 | WA | - |
testcase_07 | WA | - |
testcase_08 | WA | - |
testcase_09 | WA | - |
testcase_10 | WA | - |
testcase_11 | AC | 187 ms
11,664 KB |
testcase_12 | AC | 195 ms
11,796 KB |
testcase_13 | AC | 305 ms
15,996 KB |
testcase_14 | WA | - |
testcase_15 | AC | 342 ms
17,200 KB |
testcase_16 | AC | 154 ms
8,804 KB |
testcase_17 | AC | 2 ms
5,376 KB |
testcase_18 | AC | 23 ms
5,376 KB |
testcase_19 | WA | - |
testcase_20 | AC | 268 ms
14,772 KB |
testcase_21 | AC | 212 ms
12,608 KB |
testcase_22 | AC | 29 ms
5,376 KB |
testcase_23 | AC | 168 ms
11,124 KB |
testcase_24 | AC | 143 ms
9,960 KB |
testcase_25 | AC | 64 ms
6,268 KB |
testcase_26 | AC | 15 ms
5,376 KB |
testcase_27 | WA | - |
testcase_28 | AC | 197 ms
11,936 KB |
testcase_29 | AC | 308 ms
15,984 KB |
testcase_30 | AC | 91 ms
7,616 KB |
testcase_31 | AC | 73 ms
6,944 KB |
testcase_32 | AC | 249 ms
13,236 KB |
testcase_33 | AC | 154 ms
10,440 KB |
testcase_34 | AC | 119 ms
9,096 KB |
testcase_35 | AC | 160 ms
10,448 KB |
testcase_36 | AC | 116 ms
8,696 KB |
testcase_37 | AC | 279 ms
15,452 KB |
testcase_38 | AC | 2 ms
5,376 KB |
testcase_39 | AC | 2 ms
5,376 KB |
testcase_40 | AC | 1 ms
5,376 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; const long long MOD = 998244353; long long modsub(long long a, long long b){ return (a % MOD - b % MOD + MOD) % MOD; } long long modpow(long long a, long long b){ long long ans = 1; while (b > 0){ if (b % 2 == 1){ ans *= a; ans %= MOD; } a *= a; a %= MOD; b /= 2; } return ans; } long long modinv(long long a){ return modpow(a, MOD - 2); } vector<long long> mf = {1}; vector<long long> mfi = {1}; long long modfact(int n){ if (mf.size() > n){ return mf[n]; } else { for (int i = mf.size(); i <= n; i++){ long long next = mf.back() * i % MOD; mf.push_back(next); mfi.push_back(modinv(next)); } return mf[n]; } } long long modfactinv(int n){ if (mfi.size() > n){ return mfi[n]; } else { return modinv(modfact(n)); } } long long modbinom(int n, int k){ if (n < 0 || k < 0 || k > n){ return 0; } else { return modfact(n) * modfactinv(k) % MOD * modfactinv(n - k) % MOD; } } long long inversion_number(vector<int> &A){ int N = A.size(); vector<int> B = A; sort(B.begin(), B.end()); map<int, int> mp; for (int i = 0; i < N; i++){ mp[B[i]] = i + 1; } for (int i = 0; i < N; i++){ A[i] = mp[A[i]]; } long long ans = 0; vector<int> BIT(N + 1, 0); for (int i = 0; i < N; i++){ ans += i; int j = A[i]; while (j > 0){ ans -= BIT[j]; j -= j & -j; } j = A[i]; while (j <= N){ BIT[j]++; j += j & -j; } } return ans; } int main(){ int N; cin >> N; vector<int> A(N); for (int i = 0; i < N; i++){ cin >> A[i]; } long long all = 1; for (int i = 1; i <= N; i++){ all *= modbinom(N, i); all %= MOD; } long long E1 = 0; long long X = inversion_number(A); X %= MOD; for (int i = 2; i <= N; i++){ long long p = modbinom(N - 2, i - 2) * modinv(modbinom(N, i)) % MOD; E1 += X * p; E1 %= MOD; } long long sum = 0; long long sqsum = 0; for (int i = 1; i <= N; i++){ sum += i; sum %= MOD; sqsum += (long long) i * i; sqsum %= MOD; } long long half = modsub(sum * sum, sqsum) * modinv(2) % MOD; long long prob = half * modinv(N) % MOD * modinv(N) % MOD; map<int, int> mp; for (int i = 0; i < N; i++){ mp[A[i]] += 1; } long long Y = (long long) N * (N - 1) / 2; for (auto P : mp){ Y -= (long long) P.second * (P.second - 1) / 2; } long long E2 = prob * Y % MOD; cout << ((E1 + E2) % MOD) * all % MOD << endl; }