結果
問題 |
No.1621 Sequence Inversions
|
ユーザー |
![]() |
提出日時 | 2021-07-11 17:39:43 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,344 ms / 3,000 ms |
コード長 | 1,079 bytes |
コンパイル時間 | 297 ms |
コンパイル使用メモリ | 82,304 KB |
実行使用メモリ | 281,088 KB |
最終ジャッジ日時 | 2024-07-17 15:54:44 |
合計ジャッジ時間 | 32,998 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
mod = 998244353 f = [[[0] * 2501 for _ in range(101)] for _ in range(101)] for i in range(101): f[0][i][0] = 1 for n in range(1, 101): for k in range(101): for m in range(2501): if n * k < m: break f[n][k][m] += f[n - 1][k][m] if m - n >= 0 and k > 0: f[n][k][m] += f[n][k - 1][m - n] if f[n][k][m] >= mod: f[n][k][m] -= mod t, k = map(int, input().split()) a = list(map(int, input().split())) if k > t * (t - 1) // 2: print(0) exit() a.sort() cnt = [0] * 100 n = 0 for i in range(t - 1): cnt[n] += 1 if a[i] != a[i + 1]: n += 1 cnt[n] += 1 n += 1 dp = [[0] * 5000 for _ in range(n + 1)] dp[0][0] = 1 cursum = 0 for i in range(1, n + 1): x = cnt[i - 1] for j in range(5000): for l in range(cursum * x + 1): if j - l < 0: break dp[i][j] += f[x][cursum][l] * dp[i - 1][j - l] if dp[i][j] >= mod: dp[i][j] %= mod cursum += x print(dp[n][k])