結果

問題 No.2405 Minimal Matrix Decomposition
ユーザー chineristACchineristAC
提出日時 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])
0