Happy April Fools' Day! (Click Here!)

백준8131 Ploughing

작성일 | In

문제 링크

  • http://icpc.me/8131

문제 출처

  • 2005/2006 POI Stage2 6번

사용 알고리즘

  • 슬라이딩 윈도우

시간복잡도

  • $O(NM)$

풀이

열심히 직사각형을 깎아내다보면, 마지막에는 가로 혹은 세로 길이가 1인 직사각형을 없애게 됩니다. 마지막에 없애는 직사각형의 크기를 $1 \times K$ 또는 $K \times 1$이라고 하면, 직사각형을 총 $(N-1) + (M-K)$ 또는 $(N-K) + (N-1)$번 없애게 됩니다. $(N+M)$은 일정하므로 $K$를 최대화시키는 것이 이 문제의 목적입니다.

마지막에 가로로 길쭉한 직사각형을 하나 남겨놓는다고 합시다. 세로로 길쭉한 직사각형을 남기는 것이 최적인 경우에는 배열을 회전시켜서 한 번 더 해주면 됩니다.

마지막까지 남겨놓을 가로로 길쭉한 직사각형의 양쪽 끝 지점을 고정해주면, 그런 직사각형을 남겨놓을 수 있는지 $O(N+M)$에 판별할 수 있습니다. 총 $O(M^2)$가지의 경우가 있으므로 $O(M^2(N+M))$정도에 문제를 풀 수 있지만 N과 M이 너무 큽니다.

하지만 투포인터 느낌으로 끝 지점을 관리해주면 총 $O(M)$개의 경우만 봐주면 되고, $O(M(N+M))$에 해결할 수 있습니다.

배열을 회전시키고 이 작업을 한 번 더 해주면 $O(N^2 + M^2)$ 정도 시간에 문제를 풀 수 있습니다.

전체 코드

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
#include <bits/stdc++.h>
#define x first
#define y second
#define all(v) v.begin(), v.end()
#define compress(v) sort(all(v)), v.erase(unique(all(v)), v.end())
using namespace std;

typedef long long ll;

int n, m; ll k;
ll inp[2020][2020];
ll a[2020][2020];
ll tmp[2020][2020];

inline ll query(int x, int xx, int y, int yy){ return a[xx][yy] - a[x-1][yy] - a[xx][y-1] + a[x-1][y-1]; }

int chk(int l, int r){
    int t = r - l + 1;
    int a = 1, b = n, c = 1, d = m;
    for(int i=0; i<n+m-t-1; i++){
        if(b-a+1 > 1 && query(a, a, c, d) <= k){ a++; continue; }
        if(b-a+1 > 1 && query(b, b, c, d) <= k){ b--; continue; }
        if(c < l && query(a, b, c, c) <= k){ c++; continue; }
        if(d > r && query(a, b, d, d) <= k){ d--; continue; }
        return 0;
    }
    if(query(a, b, c, d) <= k) return 1;
    return 0;
}

int solve(){
    for(int i=1; i<=n; i++) for(int j=1; j<=m; j++) a[i][j] = inp[i][j] + a[i-1][j] + a[i][j-1] - a[i-1][j-1];
    int l = 1, r = 1;
    int ans = 1e9+7;
    while(1){
        int res = chk(l, r);
        if(res) ans = min(ans, n+m-(r-l+1));
        if(l == m && r == m) break;
        if(r == m) l++;
        else if(l == r) r++;
        else if(res) r++;
        else l++;
    }
    return ans;
}

int main(){
    ios_base::sync_with_stdio(false); cin.tie(nullptr);
    cin >> k >> m >> n;
    for(int i=1; i<=n; i++) for(int j=1; j<=m; j++) cin >> inp[i][j], a[i][j] = inp[i][j];
    int t1 = solve();
    memcpy(tmp, inp, sizeof tmp);
    for(int i=1; i<=m; i++) for(int j=1; j<=n; j++) inp[i][j] = tmp[n-j+1][i];
    swap(n, m);
    int t2 = solve();
    cout << min(t1, t2);
}