[파이썬 | BOJ | 17779] 게리맨더링2

https://www.acmicpc.net/problem/17779

17779번: 게리맨더링 2

문제

재현시의 시장 구재현은 지난 몇 년간 게리맨더링을 통해서 자신의 당에게 유리하게 선거구를 획정했다. 견제할 권력이 없어진 구재현은 권력을 매우 부당하게 행사했고, 심지어는 시의 이름도 재현시로 변경했다. 이번 선거에서는 최대한 공평하게 선거구를 획정하려고 한다.

재현시는 크기가 N×N인 격자로 나타낼 수 있다. 격자의 각 칸은 구역을 의미하고, r행 c열에 있는 구역은 (r, c)로 나타낼 수 있다. 구역을 다섯 개의 선거구로 나눠야 하고, 각 구역은 다섯 선거구 중 하나에 포함되어야 한다. 선거구는 구역을 적어도 하나 포함해야 하고, 한 선거구에 포함되어 있는 구역은 모두 연결되어 있어야 한다. 구역 A에서 인접한 구역을 통해서 구역 B로 갈 수 있을 때, 두 구역은 연결되어 있다고 한다. 중간에 통하는 인접한 구역은 0개 이상이어야 하고, 모두 같은 선거구에 포함된 구역이어야 한다.

선거구를 나누는 방법은 다음과 같다.

1. 기준점 (x, y)와 경계의 길이 d1, d2를 정한다. (d1, d2 ≥ 1, 1 ≤ x < x+d1+d2 ≤ N, 1 ≤ y-d1 < y < y+d2 ≤ N)
2. 다음 칸은 경계선이다.
   1. (x, y), (x+1, y-1), ..., (x+d1, y-d1)
   2. (x, y), (x+1, y+1), ..., (x+d2, y+d2)
   3. (x+d1, y-d1), (x+d1+1, y-d1+1), ... (x+d1+d2, y-d1+d2)
   4. (x+d2, y+d2), (x+d2+1, y+d2-1), ..., (x+d2+d1, y+d2-d1)
3. 경계선과 경계선의 안에 포함되어있는 곳은 5번 선거구이다.
4. 5번 선거구에 포함되지 않은 구역 (r, c)의 선거구 번호는 다음 기준을 따른다.
   • 1번 선거구: 1 ≤ r < x+d1, 1 ≤ c ≤ y
   • 2번 선거구: 1 ≤ r ≤ x+d2, y < c ≤ N
   • 3번 선거구: x+d1 ≤ r ≤ N, 1 ≤ c < y-d1+d2
   • 4번 선거구: x+d2 < r ≤ N, y-d1+d2 ≤ c ≤ N

아래는 크기가 7×7인 재현시를 다섯 개의 선거구로 나눈 방법의 예시이다.

x = 2, y = 4, d1 = 2, d2 = 2

x = 2, y = 5, d1 = 3, d2 = 2

x = 4, y = 3, d1 = 1, d2 = 1

구역 (r, c)의 인구는 A[r][c]이고, 선거구의 인구는 선거구에 포함된 구역의 인구를 모두 합한 값이다. 선거구를 나누는 방법 중에서, 인구가 가장 많은 선거구와 가장 적은 선거구의 인구 차이의 최솟값을 구해보자.

입력

첫째 줄에 재현시의 크기 N이 주어진다.

둘째 줄부터 N개의 줄에 N개의 정수가 주어진다. r행 c열의 정수는 A[r][c]를 의미한다.

출력

첫째 줄에 인구가 가장 많은 선거구와 가장 적은 선거구의 인구 차이의 최솟값을 출력한다.

제한

• 5 ≤ N ≤ 20
• 1 ≤ A[r][c] ≤ 100

풀이

N 제한이 20까지인데, 탐색해야 할 변수가 x, y, d1, d2 총 4개이므로 4중 for문을 돌려도 20 ^ 4 = 16만으로 시간복잡도가 충분하다

 

따라서 모든 x y d1 d2를 탐색해 나가면서 모든 구역을 나누는 수를 탐색하고, 그 중에서 가장 인구 차이가 적은 것을 고르면 된다.

 

x y d1 d2로 인해서 마름모의 꼭지점이 결정되고, 모든 꼭지점이 범위 내에 존재하면, 주어진 선거구 범위에 맞게 배열에 1, 2, 3, 4를 부여해서 지역구를 결정하고, 미리 주어진 인구수 지도에 맞게 각 선거구 별 인구수 계산을 하면 된다.

 

 

# [파이썬 | BOJ | 17779] 게리맨더링2
 
import sys
read = sys.stdin.readline
 
def myPrint(arr):
    for a in arr:
        print(a)
    print()
 
def solve():
    minAns = sys.maxsize
 
    for x in range(N):
        for y in range(N):
            for d1 in range(1, N+1):
                for d2 in range(1, N+1):
                    # 나머지 세 모서리가 범위 내에 있어야 한다.
                    if (0 <= x+d1 < N and 0 <= y-d1 < N) and (0 <= x+d2 < N and 0 <= y+d2 < N) and (0 <= x+d1+d2 < N and 0 <= y-d1+d2 < N):
                        section = [[5 for _ in range(N)] for _ in range(N)]
                        
                        #1번 구역 설정
                        for i in range(d1+1):
                            r1, c1 = x+i, y-i
                            for tr in range(r1):
                                section[tr][c1] = 1
 
                        r2, c2 = x+d1, y-d1
                        for tc in range(c2):
                            for tr in range(r2):
                                section[tr][tc] = 1
 
                        #2번 구역 설정
                        for i in range(d2+1):
                            r1, c1 = x+i, y+i
                            for tc in range(c1+1, N):
                                section[r1][tc] = 2
 
                        r2, c2 = x, y
                        for tr in range(r2):
                            for tc in range(c2+1, N):
                                section[tr][tc] = 2
 
                        #3번 구역 설정
                        for i in range(d2+1):
                            r1, c1 = x+d1+i, y-d1+i
                            for tc in range(c1):
                                section[r1][tc] = 3
 
                        r2, c2 = x+d1+d2, y-d1+d2
                        for tr in range(r2+1, N):
                            for tc in range(c2):
                                section[tr][tc] = 3
 
                        #4번 구역 설정
                        for i in range(d1+1):
                            r1, c1 = x+d2+i, y+d2-i
                            for tr in range(r1+1, N):
                                section[tr][c1] = 4
 
                        r2, c2 = x+d2, y+d2
                        for tc in range(c2+1, N):
                            for tr in range(r2+1, N):
                                section[tr][tc] = 4
 
                        cal = [0 for _ in range(6)]
 
                        for r in range(N):
                            for c in range(N):
                                sec = section[r][c]
                                cal[sec] += population[r][c]
 
                        ans = max(cal[1:]) - min(cal[1:])
                        if minAns > ans:
                            minAns = ans
                            #print(section, population, ans)
    
    print(minAns)
 
N = int(read())
 
population = []
for _ in range(N):
    population.append(list(map(int, read().split())))
solve()