結果
問題 | No.165 四角で囲え! |
ユーザー |
|
提出日時 | 2015-05-30 01:35:49 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 983 ms / 5,000 ms |
コード長 | 2,105 bytes |
コンパイル時間 | 608 ms |
コンパイル使用メモリ | 82,304 KB |
実行使用メモリ | 78,720 KB |
最終ジャッジ日時 | 2024-12-26 03:33:07 |
合計ジャッジ時間 | 10,494 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 19 |
ソースコード
def read_data():N, B = map(int, input().split())Xs = []Ys = []Ps = []for i in range(N):x, y, p = map(int, input().split())Xs.append(x)Ys.append(y)Ps.append(p)return N, B, Xs, Ys, Psdef solve(N, B, Xs, Ys, Ps):if len(Xs) == 0:return 0Xs = compress(Xs)Ys = compress(Ys)if max(Xs) > max(Ys):Xs, Ys = Ys, Xscum, cumP = preprocess(Xs, Ys, Ps)return find_max(cum, cumP, B)def compress(As):Avals = list(set(As))Avals.sort()val2idx = {a:i for i, a in enumerate(Avals)}return [val2idx[a] for a in As]def preprocess(Xs, Ys, Ps):max_x = max(Xs)max_y = max(Ys)mat = [[0] * (max_y + 2) for _ in range(max_x + 2)]for x, y, p in zip(Xs, Ys, Ps):mat[x + 1][y + 1] = pcum = [[0] * (max_y + 2) for _ in range(max_x + 2)]cumP = [[0] * (max_y + 2) for _ in range(max_x + 2)]for x in range(1, max_x + 2):cumx = cum[x]cumx_prev = cum[x - 1]cumPx = cumP[x]cumPx_prev = cumP[x - 1]matx = mat[x]tmp = 0tmpP = 0for y in range(1, max_y + 2):if matx[y]:tmp += 1tmpP += matx[y]cumx[y] = cumx_prev[y] + tmpcumPx[y] = cumPx_prev[y] + tmpPreturn cum, cumPdef find_max(cum, cumP, B):nrow = len(cum)ncol = len(cum[0])max_n = 0for t in range(nrow - 1):cum_t = cum[t]cumP_t = cumP[t]for b in range(t + 1, nrow):n = find_max_n(cum_t, cumP_t, cum[b], cumP[b], B, ncol)if n > max_n:max_n = nreturn max_ndef find_max_n(t, Pt, b, Pb, B, ncol):max_n = 0left = 0for right in range(1, ncol):p = Pb[right] - Pt[right]while p - Pb[left] + Pt[left] > B:left += 1n = b[right] - t[right] - b[left] + t[left]if n > max_n:max_n = nreturn max_nif __name__ == '__main__':N, B, Xs, Ys, Ps = read_data()print(solve(N, B, Xs, Ys, Ps))