結果
問題 |
No.2898 Update Max
|
ユーザー |
![]() |
提出日時 | 2024-09-25 15:56:28 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,172 bytes |
コンパイル時間 | 1,833 ms |
コンパイル使用メモリ | 194,420 KB |
最終ジャッジ日時 | 2025-02-24 12:12:03 |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 23 WA * 5 |
ソースコード
#include <bits/stdc++.h> using namespace std; using ll = long long; const int N = 2e5 + 10, mod = 998244353; int n, a[N], cnt[N], mmax, ct, c; bool f[N]; ll ans, fact[N], factinv[N], inv[N]; ll C(int n, int m) { return fact[n] * factinv[m] % mod * factinv[n - m] % mod; } int main() { ios::sync_with_stdio(0), cin.tie(0); cin >> n; for (int i = 1; i <= n; i++) cin >> a[i], f[a[i]] = 1; for (int i = 1; i <= n; i++) ct += a[i] == -1; inv[1] = 1; for (int i = 2; i <= n; i++) { inv[i] = (mod - mod / i) * inv[mod % i] % mod; } fact[0] = factinv[0] = 1; for (int i = 1; i <= n; i++) { fact[i] = fact[i - 1] * i % mod; factinv[i] = factinv[i - 1] * inv[i] % mod; } for (int i = 1; i <= n; i++) cnt[i] = cnt[i - 1] + !f[i]; for (int i = 1; i <= n; i++) { mmax = max(mmax, a[i]); if (a[i] == -1) c++, (ans += (C(ct, c) - C(cnt[mmax - 1], c) + mod) % mod * fact[c - 1] % mod * fact[ct - c] % mod) %= mod; else if (mmax == a[i] && cnt[mmax - 1] >= c) (ans += C(cnt[mmax - 1], c) * fact[c] % mod * fact[ct - c] % mod) %= mod; } cout << ans << '\n'; return 0; }