結果

問題 No.2445 奇行列式
コンテスト
ユーザー srjywrdnprkt
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
0