結果
問題 |
No.2898 Update Max
|
ユーザー |
![]() |
提出日時 | 2024-09-20 23:19:38 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,187 bytes |
コンパイル時間 | 307 ms |
コンパイル使用メモリ | 82,208 KB |
実行使用メモリ | 135,748 KB |
最終ジャッジ日時 | 2024-09-20 23:19:50 |
合計ジャッジ時間 | 8,937 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | WA * 3 |
other | WA * 28 |
ソースコード
mod = 998244353 n = 10**6 inv = [1 for j in range(n+1)] for a in range(2,n+1): # ax + py = 1 <=> rx + p(-x-qy) = -q => x = -(inv[r]) * (p//a) (r = p % a) res = (mod - inv[mod%a]) * (mod // a) inv[a] = res % mod fact = [1 for i in range(n+1)] for i in range(1,n+1): fact[i] = fact[i-1]*i % mod fact_inv = [1 for i in range(n+1)] fact_inv[-1] = pow(fact[-1],mod-2,mod) for i in range(n,0,-1): fact_inv[i-1] = fact_inv[i]*i % mod def binom(n,r): if n < r or n < 0 or r < 0: return 0 res = fact_inv[n-r] * fact_inv[r] % mod res *= fact[n] res %= mod return res N = int(input()) A = list(map(int,input().split())) X = [0 for i in range(N)] Q = [-1 for i in range(N + 1)] for i in range(N): if A[i] == -1: X[i] += 1 if i > 0: X[i] += X[i - 1] else: Q[A[i]] = i S = set(A) R = [] for i in range(1,N + 1): if i not in S: R.append(i) n = len(R) M = X[-1] m = M j = M ans = 0 for i in range(N,-1,-1): if Q[i] == -1: res = fact[j - 1] * fact[M - j] % mod ans = (ans + res * binom(M + 1,j) % mod) % mod j -= 1 else: n = X[Q[i]] ans = (ans + fact[m] * binom(n,m) % mod + fact[X[-1] - n]) % mod M = min(M,n) print(ans)