結果

問題 No.165 四角で囲え！
コンテスト
ユーザー nebukuro09
提出日時 2017-01-05 19:56:29
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 171 ms / 5,000 ms
コード長 2,025 bytes
コンパイル時間 1,159 ms
コンパイル使用メモリ 125,388 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-12 06:16:19
合計ジャッジ時間 3,800 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import std.stdio, std.array, std.string, std.conv, std.algorithm;
import std.typecons, std.range, std.random, std.math, std.container;
import std.numeric, std.bigint, core.bitop, core.stdc.stdio, std.bitmanip;


int acm2d (ref int[][] acm, int r1, int c1, int r2, int c2) {
    return acm[r2+1][c2+1] - acm[r2+1][c1] - acm[r1][c2+1] + acm[r1][c1];
}

void main() {
    alias Tuple!(int, "x", int, "y", int, "p") Point;
    int N, B;
    Point[] p;

    scanf("%d %d", &N, &B);
    p = new Point[](N);
    foreach (i; 0..N) {
        int x, y, po;
        scanf("%d %d %d", &x, &y, &po);
        p[i] = Point(x, y, po);
    }

    p.sort!"a.x < b.x"();
    int prev = p[0].x;
    p[0].x = 0;
    foreach (i; 1..N) {
        int cur = p[i].x;
        p[i].x = (p[i].x == prev) ? p[i-1].x : p[i-1].x + 1;
        prev = cur;
    }

    p.sort!"a.y < b.y"();
    prev = p[0].y;
    p[0].y = 0;
    foreach (i; 1..N) {
        int cur = p[i].y;
        p[i].y = (p[i].y == prev) ? p[i-1].y : p[i-1].y + 1;
        prev = cur;
    }


    auto acm_score = new int[][](N+1, N+1);
    auto acm_exist = new int[][](N+1, N+1);

    foreach (i; 0..N) {
        acm_score[p[i].x+1][p[i].y+1] = p[i].p;
        acm_exist[p[i].x+1][p[i].y+1] = 1;
    }


    foreach (i; 0..N+1) {
        foreach (j; 0..N) {
            acm_score[i][j+1] += acm_score[i][j];
            acm_exist[i][j+1] += acm_exist[i][j];
        }
    }

    foreach (j; 0..N+1) {
        foreach (i; 0..N) {
            acm_score[i+1][j] += acm_score[i][j];
            acm_exist[i+1][j] += acm_exist[i][j];
        }
    }

    int ans = 0;
    foreach (r1; 0..N) {
        foreach (r2; r1..N) {
            int c1, c2;
            while (c2 < N) {
                if (acm2d(acm_score, r1, c1, r2, c2) <= B) {
                    ans = max(ans, acm2d(acm_exist, r1, c1, r2, c2));
                    c2++;
                } else{
                    c1++;
                    if (c2 < c1) c2 = c1;
                }
            }
        }
    }

    ans.writeln;
}
0