結果
問題 | No.2263 Perms |
ユーザー | titia |
提出日時 | 2023-04-08 01:45:23 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
AC
|
実行時間 | 225 ms / 2,000 ms |
コード長 | 2,718 bytes |
コンパイル時間 | 77 ms |
コンパイル使用メモリ | 13,184 KB |
実行使用メモリ | 11,392 KB |
最終ジャッジ日時 | 2024-10-02 22:14:15 |
合計ジャッジ時間 | 4,825 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 29 ms
11,008 KB |
testcase_01 | AC | 26 ms
11,136 KB |
testcase_02 | AC | 28 ms
11,008 KB |
testcase_03 | AC | 28 ms
11,008 KB |
testcase_04 | AC | 32 ms
10,880 KB |
testcase_05 | AC | 32 ms
11,008 KB |
testcase_06 | AC | 155 ms
11,136 KB |
testcase_07 | AC | 116 ms
11,264 KB |
testcase_08 | AC | 90 ms
11,136 KB |
testcase_09 | AC | 187 ms
11,264 KB |
testcase_10 | AC | 225 ms
11,136 KB |
testcase_11 | AC | 122 ms
11,264 KB |
testcase_12 | AC | 128 ms
11,264 KB |
testcase_13 | AC | 165 ms
11,136 KB |
testcase_14 | AC | 160 ms
11,264 KB |
testcase_15 | AC | 143 ms
11,264 KB |
testcase_16 | AC | 37 ms
11,008 KB |
testcase_17 | AC | 40 ms
11,136 KB |
testcase_18 | AC | 36 ms
11,008 KB |
testcase_19 | AC | 42 ms
11,136 KB |
testcase_20 | AC | 117 ms
11,264 KB |
testcase_21 | AC | 126 ms
11,008 KB |
testcase_22 | AC | 109 ms
11,008 KB |
testcase_23 | AC | 61 ms
11,136 KB |
testcase_24 | AC | 39 ms
11,008 KB |
testcase_25 | AC | 126 ms
11,392 KB |
testcase_26 | AC | 141 ms
11,136 KB |
testcase_27 | AC | 57 ms
11,008 KB |
testcase_28 | AC | 40 ms
11,264 KB |
testcase_29 | AC | 48 ms
11,264 KB |
testcase_30 | AC | 33 ms
11,008 KB |
testcase_31 | AC | 31 ms
11,136 KB |
testcase_32 | AC | 39 ms
11,008 KB |
testcase_33 | AC | 37 ms
11,008 KB |
testcase_34 | AC | 29 ms
10,880 KB |
testcase_35 | AC | 30 ms
11,008 KB |
testcase_36 | AC | 32 ms
11,008 KB |
testcase_37 | AC | 31 ms
11,008 KB |
testcase_38 | AC | 107 ms
11,264 KB |
testcase_39 | AC | 28 ms
11,008 KB |
testcase_40 | AC | 29 ms
11,136 KB |
ソースコード
import sys input = sys.stdin.readline from collections import defaultdict,deque N,M=map(int,input().split()) A=[list(map(int,input().split())) for i in range(N)] ALIST=[] for tests in range(M): # dinic法 V=N*2+2 EDGE=[defaultdict(int) for i in range(V)] for i in range(N): for j in range(N): if A[i][j]>0: EDGE[i][N+j]=1 start=N*2 goal=N*2+1 for i in range(N): EDGE[start][i]=1 for i in range(N): EDGE[N+i][goal]=1 #print(EDGE) ANS=0 while True: # まずBFSする DIS=[-1]*V Q=deque([start]) DIS[start]=0 EDGE2=[[] for i in range(V)] while Q: x=Q.popleft() for to in EDGE[x]: if EDGE[x][to]==0: continue if DIS[to]==-1: DIS[to]=DIS[x]+1 Q.append(to) EDGE2[x].append(to) elif DIS[to]==DIS[x]+1: EDGE2[x].append(to) if DIS[goal]==-1: break # BFSしたときのEDGEを使ってDFSする MINCOST=[float("inf")]*V NOW=start ROUTE=[-1]*V while NOW!=-1: # DFS cost=MINCOST[NOW] if NOW==goal: ANS+=cost i=goal while i!=start: # goalからたどり,Routeを使ってEDGEの更新 j=ROUTE[i] if EDGE[j][i]==cost: NOW=j EDGE2[j].pop() EDGE[j][i]-=cost # 使ったルートをいけなく更新 MINCOST[j]-=cost EDGE[i][j]+=cost # 逆向きに進めるようにする. i=j continue if EDGE2[NOW]: to=EDGE2[NOW][-1] ROUTE[to]=NOW MINCOST[to]=min(cost,EDGE[NOW][to]) NOW=to else: if NOW==start: break EDGE2[ROUTE[NOW]].pop() NOW=ROUTE[NOW] #print(ANS) ANS=[-1]*N for i in range(N): for j in range(N): if EDGE[i+N][j]>0: ANS[j]=i+1 A[j][i]-=1 break if -1 in ANS: print(-1) exit() ALIST.append(ANS) for i in range(N): for j in range(N): if A[i][j]>0: print(-1) exit() for ans in ALIST: print(*ans)