結果

問題 No.165 四角で囲え！
ユーザー nebukuro09nebukuro09
提出日時 2017-01-05 17:06:21
言語 D
(dmd 2.106.1)
結果
WA  
実行時間 -
コード長 2,521 bytes
コンパイル時間 857 ms
コンパイル使用メモリ 106,536 KB
実行使用メモリ 6,220 KB
最終ジャッジ日時 2023-09-03 00:24:13
合計ジャッジ時間 2,140 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 AC 1 ms
4,500 KB
testcase_02 AC 1 ms
4,376 KB
testcase_03 AC 1 ms
4,376 KB
testcase_04 AC 1 ms
4,380 KB
testcase_05 WA -
testcase_06 AC 4 ms
4,380 KB
testcase_07 AC 1 ms
4,376 KB
testcase_08 AC 1 ms
4,376 KB
testcase_09 AC 1 ms
4,376 KB
testcase_10 AC 1 ms
4,376 KB
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 AC 40 ms
4,416 KB
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

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 h, int w, int i, int j) {
    return acm[i+h][j+w] - acm[i+h][j] - acm[i][j+w] + acm[i][j];
}

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 hi = N*N+1;
    int lo = 0;
    while (hi - lo > 1) {
        int mid = (hi + lo) / 2;
        bool ok;
        a:foreach (h; 1..mid+1) {
            if (mid % h != 0) continue;
            int w = mid / h;
            foreach (i; 0..N-h+1) {
                foreach (j; 0..N-w+1) {
                    if (acm2d(acm_score, h, w, i, j) <= B &&
                        acm2d(acm_exist, h, w, i, j) > 0) {
                        ok = true;
                        break a;
                    }
                }
            }
        }
        if (ok) lo = mid;
        else hi = mid;
    }

    int ans = 0;
    foreach (h; 1..lo+1) {
        if (lo % h != 0) continue;
        int w = lo / h;
        foreach (i; 0..N-h+1) {
            foreach (j; 0..N-w+1) {
                if (acm2d(acm_score, h, w, i, j) <= B)
                    ans = max(ans, acm2d(acm_exist, h, w, i, j));
            }
        }
    }

    ans.writeln;
}
0