結果

問題 No.753 最強王者決定戦
コンテスト
ユーザー Pachicobue
提出日時 2018-11-11 02:24:22
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 472 ms / 1,000 ms
コード長 1,524 bytes
コンパイル時間 1,982 ms
コンパイル使用メモリ 193,416 KB
最終ジャッジ日時 2025-01-06 16:24:18
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 4
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using ll = long long;
constexpr int N = 16;
constexpr int MASK = 1 << N;
int win[N][N];
ll dp[N][MASK];
int bit[MASK][N];
int pc[MASK];
int main()
{
    for (int i = 0; i < MASK; i++) {
        for (int j = 0; j < N; j++) {
            if (i & (1 << j)) { bit[i][pc[i]++] = j; }
        }
    }
    for (int i = 0; i < N; i++) {
        for (int j = 0, c; j < N; j++) {
            std::cin >> c;
            if (i < j) { win[i][j] = c == 1, win[j][i] = c != 1; }
        }
    }
    auto next_permutation = [](const int mask) {
        const int x = mask & -mask, y = mask + x;
        return (((mask & ~y) / x) >> 1) | y;
    };
    for (int i = 0; i < N; i++) { dp[i][1 << i] = 1; }
    for (int i = 1; i < MASK; i++) {
        const int b = pc[i], h = b / 2;
        if (b % 2 == 1) { continue; }
        int ind = (1 << h) - 1;
        for (; ind < (1 << (b - 1)); ind = next_permutation(ind)) {
            int sub = 0;
            for (int j = 0; j < b; j++) {
                if (ind & (1 << j)) { sub += (1 << bit[i][j]); }
            }
            const int other = i - sub;
            for (int ji = 0; ji < b / 2; ji++) {
                const int j = bit[sub][ji];
                for (int ki = 0; ki < b / 2; ki++) {
                    const int k = bit[other][ki];
                    dp[win[j][k] ? j : k][i] += dp[j][sub] * dp[k][other] * 2;
                }
            }
        }
    }
    for (int i = 0; i < N; i++) { std::cout << dp[i][MASK - 1] << "\n"; }
    return 0;
}
0