結果
問題 |
No.2445 奇行列式
|
ユーザー |
👑 ![]() |
提出日時 | 2023-08-25 22:49:14 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,076 ms / 3,000 ms |
コード長 | 959 bytes |
コンパイル時間 | 405 ms |
コンパイル使用メモリ | 82,432 KB |
実行使用メモリ | 100,608 KB |
最終ジャッジ日時 | 2024-12-24 10:15:10 |
合計ジャッジ時間 | 8,419 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
def solve(): N, B = map(int, input().split()) A = [None] * N for i in range(N): A[i] = list(map(int, input().split())) A[i] = [a % B for a in A[i]] popcount = [0] * (1 << N) for S in range(1, 1 << N): x = S & (-S) popcount[S] = popcount[S ^ x] + 1 bit = lambda S,k : (S>>k) & 1 DP0 = [0] * (1 << N); DP0[0] = 1 DP1 = [0] * (1 << N) for S in range(1 << N): t = 0 i = popcount[S] for k in range(N): if not bit(S, k): T = S | (1 << k) if t % 2 == 0: DP0[T] += DP0[S] * A[i][k]; DP0[T] %= B DP1[T] += DP1[S] * A[i][k]; DP1[T] %= B else: DP0[T] += DP1[S] * A[i][k]; DP0[T] %= B DP1[T] += DP0[S] * A[i][k]; DP1[T] %= B t += 1 return DP1[-1] #================================================== print(solve())