結果

問題 No.2898 Update Max
ユーザー MitarushiMitarushi
提出日時 2024-09-19 21:42:24
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 473 ms / 2,000 ms
コード長 1,604 bytes
コンパイル時間 144 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 221,500 KB
最終ジャッジ日時 2024-09-19 21:50:18
合計ジャッジ時間 9,891 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 103 ms
135,552 KB
testcase_01 AC 102 ms
135,424 KB
testcase_02 AC 87 ms
135,296 KB
testcase_03 AC 130 ms
152,704 KB
testcase_04 AC 126 ms
152,576 KB
testcase_05 AC 124 ms
152,924 KB
testcase_06 AC 124 ms
152,704 KB
testcase_07 AC 126 ms
152,448 KB
testcase_08 AC 222 ms
220,960 KB
testcase_09 AC 240 ms
221,500 KB
testcase_10 AC 241 ms
221,496 KB
testcase_11 AC 223 ms
221,112 KB
testcase_12 AC 218 ms
221,388 KB
testcase_13 AC 264 ms
215,420 KB
testcase_14 AC 264 ms
215,924 KB
testcase_15 AC 271 ms
215,276 KB
testcase_16 AC 305 ms
215,672 KB
testcase_17 AC 290 ms
215,420 KB
testcase_18 AC 389 ms
195,052 KB
testcase_19 AC 385 ms
193,016 KB
testcase_20 AC 386 ms
194,036 KB
testcase_21 AC 401 ms
194,308 KB
testcase_22 AC 367 ms
195,168 KB
testcase_23 AC 469 ms
169,624 KB
testcase_24 AC 457 ms
169,880 KB
testcase_25 AC 449 ms
169,892 KB
testcase_26 AC 458 ms
169,492 KB
testcase_27 AC 473 ms
169,884 KB
testcase_28 AC 319 ms
169,360 KB
testcase_29 AC 87 ms
135,680 KB
testcase_30 AC 89 ms
135,296 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import bisect

n = int(input())
a = list(map(int, input().split()))
a_set = set(a)

assert 1 <= n <= 2 * 10**5
assert len(a) == n
assert all(x == -1 or 1 <= x <= n for x in a)
assert len(a_set - {-1}) == sum(x != -1 for x in a)

mod = 998244353
max_frac = 10 ** 6
frac = [1]
for i in range(1, max_frac + 1):
    frac.append(frac[-1] * i % mod)

def comb(a, b):
    if a < b or a < 0 or b < 0:
        return 0
    return frac[a] * pow(frac[b] * frac[a - b], -1, mod) % mod

a_set = set(a)
not_a = sorted([i for i in range(1, n + 1) if i not in a_set])
not_a_len = len(not_a)

prefix_max = 0
prefix_neg = 0
ans = 0
for i, x in enumerate(a):
    if x == -1:
        not_a_pos = max(bisect.bisect_left(not_a, prefix_max), prefix_neg)
        ans += (comb(not_a_len, prefix_neg + 1) - comb(not_a_pos, prefix_neg + 1)) * frac[not_a_len - prefix_neg - 1] * frac[prefix_neg]
        # print(i, x, (comb(not_a_len + 2, not_a_pos - prefix_neg + 2) - comb(not_a_pos + 1, not_a_pos - prefix_neg + 2)) * frac[not_a_len - prefix_neg] * frac[prefix_neg])
        # print(f"{i=} {x=} {not_a_pos=} {prefix_neg=} {ans=}")
        # print(comb(not_a_len, prefix_neg + 1) - comb(not_a_pos, prefix_neg), frac[not_a_len - prefix_neg - 1], frac[prefix_neg])
    else:
        if prefix_max < x:
            not_a_pos = bisect.bisect_left(not_a, x)
            ans += comb(not_a_pos, prefix_neg) * frac[not_a_len - prefix_neg] * frac[prefix_neg]
        # print(i, x, comb(not_a_pos, prefix_neg) * frac[not_a_len - prefix_neg] * frac[prefix_neg])
    prefix_max = max(prefix_max, x)
    prefix_neg += x == -1

print(ans % mod)
0