結果
問題 | No.2825 Sum of Scores of Sets of Specified Sections |
ユーザー |
![]() |
提出日時 | 2024-07-27 04:51:51 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 417 ms / 3,000 ms |
コード長 | 1,628 bytes |
コンパイル時間 | 534 ms |
コンパイル使用メモリ | 82,208 KB |
実行使用メモリ | 80,884 KB |
最終ジャッジ日時 | 2024-07-27 04:51:58 |
合計ジャッジ時間 | 6,175 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 40 ms
53,212 KB |
testcase_01 | AC | 417 ms
80,280 KB |
testcase_02 | AC | 388 ms
80,688 KB |
testcase_03 | AC | 390 ms
80,884 KB |
testcase_04 | AC | 110 ms
76,952 KB |
testcase_05 | AC | 54 ms
67,856 KB |
testcase_06 | AC | 45 ms
60,572 KB |
testcase_07 | AC | 79 ms
76,744 KB |
testcase_08 | AC | 68 ms
71,896 KB |
testcase_09 | AC | 55 ms
67,316 KB |
testcase_10 | AC | 55 ms
67,088 KB |
testcase_11 | AC | 50 ms
62,636 KB |
testcase_12 | AC | 46 ms
60,544 KB |
testcase_13 | AC | 56 ms
67,188 KB |
testcase_14 | AC | 115 ms
77,180 KB |
testcase_15 | AC | 116 ms
77,280 KB |
testcase_16 | AC | 115 ms
77,060 KB |
testcase_17 | AC | 119 ms
77,232 KB |
testcase_18 | AC | 116 ms
77,348 KB |
testcase_19 | AC | 115 ms
77,276 KB |
testcase_20 | AC | 115 ms
77,088 KB |
testcase_21 | AC | 112 ms
77,168 KB |
testcase_22 | AC | 115 ms
77,520 KB |
testcase_23 | AC | 117 ms
77,580 KB |
testcase_24 | AC | 86 ms
76,580 KB |
testcase_25 | AC | 83 ms
76,620 KB |
testcase_26 | AC | 82 ms
76,456 KB |
testcase_27 | AC | 87 ms
76,456 KB |
testcase_28 | AC | 40 ms
53,608 KB |
testcase_29 | AC | 38 ms
53,244 KB |
testcase_30 | AC | 39 ms
53,504 KB |
testcase_31 | AC | 39 ms
53,076 KB |
testcase_32 | AC | 209 ms
78,056 KB |
ソースコード
# I:{1,...,H}の部分集合,J:{1,...,W}の部分集合でI,Jのsizeが同じ組み合わせすべてに対して # det(A') * det(B')の総和が答え.X':Iに含まれる行,Jに含まれる列としたXの部分行列. # 解説 # コーシー・ビネの定理 https://ja.wikipedia.org/wiki/%E3%82%B3%E3%83%BC%E3%82%B7%E3%83%BC%E3%83%BB%E3%83%93%E3%83%8D%E3%81%AE%E5%85%AC%E5%BC%8F # X = [[I,B^t],[-A,I]]とするとIだけででてくる1を引けばいい. mod = 998244353 def mul(A,B): l = len(A) m = len(A[0]) n = len(B[0]) C = [[0] * n for i in range(l)] for i in range(l): for j in range(n): for k in range(m): C[i][j] += A[i][k] * B[k][j] C[i][j] %= mod return C def det(A): n = len(A) res = 1 for i in range(n): if A[i][i] == 0: for j in range(i+1,n): if A[i][j] != 0: for k in range(i,n): A[i][k],A[j][k] = A[j][k],A[i][k] res *= -1 res %= mod break if A[i][i] == 0: return 0 a = A[i][i] res = res * a % mod a_inv = pow(a,mod - 2,mod) for j in range(i,n): A[i][j] = A[i][j] * a_inv % mod for j in range(i+1,n): a = A[j][i] for k in range(i,n): A[j][k] -= A[i][k] * a % mod A[j][k] %= mod return res H,W = map(int,input().split()) A = [list(map(int,input().split())) for h in range(H)] B = [[0 for h in range(H)] for w in range(W)] for h in range(H): b = list(map(int,input().split())) for w in range(W): B[w][h] = b[w] C = mul(A,B) if H > W: C = mul(A,B) for h in range(H): C[h][h] += 1 print((det(C) - 1) % mod)