結果
問題 | No.2445 奇行列式 |
ユーザー |
![]() |
提出日時 | 2023-08-25 22:53:25 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 774 ms / 3,000 ms |
コード長 | 990 bytes |
コンパイル時間 | 2,000 ms |
コンパイル使用メモリ | 197,312 KB |
最終ジャッジ日時 | 2025-02-16 14:16:04 |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 20 |
ソースコード
#include <bits/stdc++.h> using namespace std; using ll = long long; int main(){ ll N, B, M, x, c, d; cin >> N >> B; vector<vector<ll>> A(N, vector<ll>(N)); for (int i=0; i<N; i++){ for (int j=0; j<N; j++){ cin >> A[i][j]; A[i][j] %= B; } } M = 1<<N; vector<vector<ll>> dp(M, vector<ll>(2)); dp[0][0] = 1; for (int i=0; i<N; i++){ x = M-1 - ((1<<i)-1); } for (int i=0; i<M; i++){ for (int j=0; j<N; j++){ if (i & 1<<j) continue; x = M-1 - ((1<<j)-1); c = __builtin_popcount(x & i); //(x以上のものでiに含まれるものの数) d = __builtin_popcount(i); //(iに含まれるものの要素数) for (int k=0; k<2; k++){ dp[i|(1<<j)][(k+c)%2] += (dp[i][k] * A[d][j]) % B; dp[i|(1<<j)][(k+c)%2] %= B; } } } cout << dp[M-1][1] << endl; return 0; }