結果
問題 | No.2255 Determinant Sum |
ユーザー |
![]() |
提出日時 | 2023-03-20 05:39:04 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 306 ms / 2,000 ms |
コード長 | 1,134 bytes |
コンパイル時間 | 263 ms |
コンパイル使用メモリ | 82,148 KB |
実行使用メモリ | 77,736 KB |
最終ジャッジ日時 | 2024-09-18 16:36:39 |
合計ジャッジ時間 | 3,530 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 23 |
ソースコード
def det(M, mod):M = [row[:] for row in M]N, sign, prev = len(M), 1, 1for i in range(N-1):#print(*M, sep="\n")#print()if M[i][i] == 0:swapto = next((j for j in range(i+1, N) if M[j][i] != 0), None)if swapto is None:return 0M[i], M[swapto], sign = M[swapto], M[i], -signprevinv = pow(prev, mod-2, mod)for j in range(i+1, N):for k in range(i+1, N):M[j][k] = (M[j][k]*M[i][i]%mod-M[j][i]*M[i][k]%mod)*previnv%modprev = M[i][i]return sign*M[-1][-1]%modT = int(input())for _ in range(T):n, p = map(int,input().split())a = [list(map(int,input().split())) for i in range(n)]if p > 2:print(0)continue# p = 2x = [0] * ny = [0] * nfor i in range(n):for j in range(n):if a[i][j] == -1:x[i] += 1y[j] += 1mode = 1for i in range(n):if x[i] > 1 or y[i] > 1:mode = 0breakif mode == 0:print(0)continuebx = []by = []for i in range(n):if x[i] == 0:bx.append(i)if y[i] == 0:by.append(i)m = len(bx)b = [[a[bx[i]][by[j]] for j in range(m)]for i in range(m)]if len(b) > 0:print(det(b, 2))else:print(1)