結果
問題 |
No.2520 L1 Explosion
|
ユーザー |
![]() |
提出日時 | 2023-03-28 22:46:32 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 376 ms / 2,000 ms |
コード長 | 1,155 bytes |
コンパイル時間 | 418 ms |
コンパイル使用メモリ | 82,304 KB |
実行使用メモリ | 148,224 KB |
最終ジャッジ日時 | 2024-07-26 14:22:41 |
合計ジャッジ時間 | 6,622 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 |
ソースコード
import sys input = sys.stdin.readline MOD = 998244353 def solve(): n, m = map(int, input().split()) X1 = []; X2 = []; Xst = set([]) Y1 = []; Y2 = []; Yst = set([]) for x, y, h in (map(int, input().split()) for _ in [0] * n): d = m - h x1, y1 = (x - d) - y, (x - d) + y x2, y2 = (x + d) - y, (x + d) + y X1.append(x1); X2.append(x2); Xst |= set([x1, x2]) Y1.append(y1); Y2.append(y2); Yst |= set([y1, y2]) X = sorted(Xst); Y = sorted(Yst) hX = {}; hY = {} for i, (x, y) in enumerate(zip(X, Y)): hX[x] = hY[y] = i h, w = len(hX), len(hY) dp = [[0] * w for _ in [0] * h] for x1, y1, x2, y2 in zip(X1, Y1, X2, Y2): dp[hX[x1]][hY[y1]] += 1 dp[hX[x1]][hY[y2]] -= 1 dp[hX[x2]][hY[y1]] -= 1 dp[hX[x2]][hY[y2]] += 1 for i in range(h): for j in range(1, w): dp[i][j] += dp[i][j - 1] for i in range(1, h): for j in range(w): dp[i][j] += dp[i - 1][j] ans = [0] * (n + 1) for i in range(h - 1): for j in range(w - 1): ans[dp[i][j]] += ((X[i + 1] - X[i]) % MOD) * ((Y[j + 1] - Y[j]) % MOD) % MOD ans[dp[i][j]] %= MOD div2 = pow(2, MOD - 2, MOD) for i in range(1, n + 1): print(ans[i] * div2 % MOD) solve()