結果
問題 | No.241 出席番号(1) |
ユーザー |
|
提出日時 | 2020-09-25 23:26:35 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 36 ms / 2,000 ms |
コード長 | 2,480 bytes |
コンパイル時間 | 104 ms |
コンパイル使用メモリ | 13,056 KB |
実行使用メモリ | 11,648 KB |
最終ジャッジ日時 | 2024-06-28 08:08:59 |
合計ジャッジ時間 | 2,591 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 29 |
ソースコード
import sysimport mathfrom collections import defaultdictfrom collections import dequesys.setrecursionlimit(1000000)MOD = 10 ** 9 + 7input = lambda: sys.stdin.readline().strip()NI = lambda: int(input())NMI = lambda: map(int, input().split())NLI = lambda: list(NMI())SI = lambda: input()INF = 10 ** 20class Dinic:""" 最大流(Dinic) """def __init__(self, n):self.n = nself.links = [[] for _ in range(n)]self.depth = Noneself.progress = Nonedef add_link(self, _from, to, cap):self.links[_from].append([cap, to, len(self.links[to])])self.links[to].append([0, _from, len(self.links[_from]) - 1])def bfs(self, s):from collections import dequedepth = [-1] * self.ndepth[s] = 0q = deque([s])while q:v = q.popleft()for cap, to, rev in self.links[v]:if cap > 0 and depth[to] < 0:depth[to] = depth[v] + 1q.append(to)self.depth = depthdef dfs(self, v, t, flow):if v == t:return flowlinks_v = self.links[v]for i in range(self.progress[v], len(links_v)):self.progress[v] = icap, to, rev = link = links_v[i]if cap == 0 or self.depth[v] >= self.depth[to]:continued = self.dfs(to, t, min(flow, cap))if d == 0:continuelink[0] -= dself.links[to][rev][0] += dreturn dreturn 0def max_flow(self, s, t):flow = 0while True:self.bfs(s)if self.depth[t] < 0:return flowself.progress = [0] * self.ncurrent_flow = self.dfs(s, t, INF)while current_flow > 0:flow += current_flowcurrent_flow = self.dfs(s, t, INF)def main():N = NI()A = [NI() for _ in range(N)]s, t = 2*N, 2*N+1din = Dinic(2*N+2)for i in range(N):din.add_link(s, i, 1)din.add_link(i+N, t, 1)hate = A[i]for j in range(N):if j == hate: continuedin.add_link(i, j+N, 1)if din.max_flow(s, t) != N:print(-1)exit()for i in range(N):for cap, to, _ in din.links[i]:if cap == 0:print(to - N)if __name__ == "__main__":main()