結果
問題 | No.2446 完全列 |
ユーザー |
![]() |
提出日時 | 2023-08-25 23:24:52 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 42 ms / 2,000 ms |
コード長 | 1,233 bytes |
コンパイル時間 | 275 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 11,392 KB |
最終ジャッジ日時 | 2024-12-24 14:23:31 |
合計ジャッジ時間 | 2,472 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 30 |
ソースコード
import sysfrom fractions import Fraction# print(sys.stdin.buffer.read())inp = map(int, sys.stdin.buffer.read().split())# l, m, n = map(int, input().split())l = next(inp)m = next(inp)n = next(inp)# N->M->L# A: M->L ... M * L# B: N->M ... N * Ma = [[0] * l for _ in range(m)]b = [[0] * m for _ in range(n)]for j in range(l):for i in range(m):a[i][j] = Fraction(next(inp), 1)for j in range(m):for i in range(n):b[i][j] = Fraction(next(inp), 1)def get_base(b):b = [bi[:] for bi in b]if not b:return bn = len(b)m = len(b[0])i = 0for j in range(m):for i0 in range(i, n):if b[i0][j]:breakelse:continueb[i], b[i0] = b[i0], b[i]c = 1 / b[i][j]for j0 in range(j, m):b[i][j0] *= cfor i0 in range(i + 1, n):c = b[i0][j]for j0 in range(j, m):b[i0][j0] -= b[i][j0] * ci += 1b = b[:i]return bb = get_base(b)for u in b:if any(sum(u[i] * a[i][j] for i in range(m)) for j in range(l)):print('No')exit()ok = m == len(b) + len(get_base(a))print('Yes' if ok else 'No')