結果

問題 No.2405 Minimal Matrix Decomposition
ユーザー chineristAC
提出日時 2023-08-04 22:24:32
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,242 bytes
コンパイル時間 195 ms
コンパイル使用メモリ 82,304 KB
実行使用メモリ 81,676 KB
最終ジャッジ日時 2024-10-14 20:24:20
合計ジャッジ時間 9,001 ms
ジャッジサーバーID
(参考情報)
judge4 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 44 WA * 2
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

import sys
input = lambda :sys.stdin.readline().rstrip()
mod = int(input())
N,M = map(int,input().split())
A = [list(map(int,input().split())) for i in range(N)]
original = [[A[i][j] for j in range(M)] for i in range(N)]
rank = 0
B = [[] for i in range(N)]
C = []
for j in range(M):
pivot = -1
for i in range(j,N):
if A[i][j]!=0:
pivot = i
break
if pivot == -1:
continue
rank += 1
C.append([A[pivot][j] for j in range(M)])
inv = pow(A[pivot][j],mod-2,mod)
for i in range(N):
if i == pivot:
continue
t = A[i][j] * inv % mod
B[i].append(t)
for k in range(j,M):
A[i][k] -= A[pivot][k] * t % mod
A[i][k] %= mod
#print(A[i])
B[pivot].append(1)
for j in range(M):
A[pivot][j] = 0
#print()
if rank == 0:
print(2)
print(N,1)
for i in range(N):
print(0)
print(1,M)
print(*[0]*M)
exit()
if N*M <= rank * (N+M):
print(1)
print(N,M)
for i in range(N):
print(*original[i])
else:
print(2)
print(N,rank)
for i in range(N):
print(*B[i])
print(rank,M)
for i in range(rank):
print(*C[i])
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0