오늘은 스타트링크에 다니는 사람들이 모여서 축구를 해보려고 한다. 축구는 평일 오후에 하고 의무 참석도 아니다. 축구를 하기 위해 모인 사람은 총 N명이고 신기하게도 N은 짝수이다. 이제 N/2명으로 이루어진 스타트 팀과 링크 팀으로 사람들을 나눠야 한다.
BOJ를 운영하는 회사 답게 사람에게 번호를 1부터 N까지로 배정했고, 아래와 같은 능력치를 조사했다. 능력치 Sij는 i번 사람과 j번 사람이 같은 팀에 속했을 때, 팀에 더해지는 능력치이다. 팀의 능력치는 팀에 속한 모든 쌍의 능력치 Sij의 합이다. Sij는 Sji와 다를 수도 있으며, i번 사람과 j번 사람이 같은 팀에 속했을 때, 팀에 더해지는 능력치는 Sij와 Sji이다.
N=4이고, S가 아래와 같은 경우를 살펴보자.
i\j 1 2 3 4
1 1 2 3
2 4 5 6
3 7 1 2
4 3 4 5
예를 들어, 1, 2번이 스타트 팀, 3, 4번이 링크 팀에 속한 경우에 두 팀의 능력치는 아래와 같다.
스타트 팀: S12 + S21 = 1 + 4 = 5
링크 팀: S34 + S43 = 2 + 5 = 7
1, 3번이 스타트 팀, 2, 4번이 링크 팀에 속하면, 두 팀의 능력치는 아래와 같다.
스타트 팀: S13 + S31 = 2 + 7 = 9
링크 팀: S24 + S42 = 6 + 4 = 10
축구를 재미있게 하기 위해서 스타트 팀의 능력치와 링크 팀의 능력치의 차이를 최소로 하려고 한다. 위의 예제와 같은 경우에는 1, 4번이 스타트 팀, 2, 3번 팀이 링크 팀에 속하면 스타트 팀의 능력치는 6, 링크 팀의 능력치는 6이 되어서 차이가 0이 되고 이 값이 최소이다.
첫째 줄에 N(4 ≤ N ≤ 20, N은 짝수)이 주어진다. 둘째 줄부터 N개의 줄에 S가 주어진다. 각 줄은 N개의 수로 이루어져 있고, i번 줄의 j번째 수는 Sij 이다. Sii는 항상 0이고, 나머지 Sij는 1보다 크거나 같고, 100보다 작거나 같은 정수이다.
첫째 줄에 스타트 팀과 링크 팀의 능력치의 차이의 최솟값을 출력한다.
왠지 수학적인 편법이 있을거 같지만,
백트래킹을 공부하기 위한 것이니 정공법으로 가자.
일단 재귀로 조합을 구한 다음, 차이값을 비교해 최소값을 도출하자.
from itertools import combinations
import sys
input = sys.stdin.readline
n = int(input())
grid = tuple(tuple(map(int, input().split())) for _ in range(n))
start, link = [], []
ans = sys.maxsize
def solve(depth, idx):
global ans
if depth == n//2:
link = [i for i in range(n) if i not in start]
spower, lpower = 0,0
for x,y in combinations(start,2):
spower += grid[x][y] + grid[y][x]
for x,y in combinations(link,2):
lpower += grid[x][y] + grid[y][x]
ans = min(ans, abs(spower-lpower))
return
for i in range(idx,n):
start.append(i)
solve(depth+1, i+1)
start.pop()
solve(0,0)
print(ans)
재귀함수 대신 itertools
의 combinations을 주로 이용한 풀이다.
from itertools import combinations
import sys
input = sys.stdin.readline
n = int(input())
grid = tuple(tuple(map(int, input().split())) for _ in range(n))
start, link = [], []
ans = sys.maxsize
for cb in combinations(range(n),n//2):
start = cb
link = tuple(set(range(n)) - set(cb))
spower, lpower = 0,0
for x,y in combinations(start,2):
spower += grid[x][y] + grid[y][x]
for x,y in combinations(link,2):
lpower += grid[x][y] + grid[y][x]
ans = min(ans, abs(spower-lpower))
print(ans)