結果

問題 No.753 最強王者決定戦
ユーザー tsutajtsutaj
提出日時 2018-11-10 11:30:30
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 649 ms / 1,000 ms
コード長 2,344 bytes
コンパイル時間 693 ms
コンパイル使用メモリ 78,616 KB
実行使用メモリ 4,376 KB
最終ジャッジ日時 2023-08-15 10:54:15
合計ジャッジ時間 4,274 ms
ジャッジサーバーID
(参考情報)
judge15 / judge13
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 611 ms
4,376 KB
testcase_01 AC 635 ms
4,376 KB
testcase_02 AC 523 ms
4,376 KB
testcase_03 AC 649 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <iostream>
#include <vector>
#include <algorithm>
#include <bitset>
using namespace std;
using ll = long long int;

const int N = 16;
ll match[20][20], dp[1 << 16], frac[20];

int main() {
    frac[0] = 1;
    for(int i=1; i<=N; i++) {
        frac[i] = frac[i-1] * i;
    }

    for(int i=0; i<N; i++) {
        for(int j=0; j<N; j++) {
            cin >> match[i][j];
            if(i > j) match[i][j] = -match[j][i];
        }
    }

    dp[(1<<N) - 1] = 1;
    for(int bit=(1<<N)-1; bit>0; bit--) {
        if(dp[bit] == 0) continue;
        int popBit = __builtin_popcount(bit);

        vector<int> alive;
        for(int i=0; i<N; i++) if(bit >> i & 1) alive.push_back(i);

        vector<int> perm(popBit);
        for(int i=0; i<popBit; i++) {
            if(i < popBit / 2) perm[i] = 0;
            else perm[i] = 1;
        }

        do {
            // 勝者と敗者
            vector<int> winner, loser;
            int par = 0;
            for(int i=0; i<popBit; i++) {
                if(perm[i] == 1) {
                    winner.push_back(alive[i]);
                    par |= (1 << alive[i]);
                }
                else {
                    loser.push_back(alive[i]);
                }
            }

            // 敗者の集合
            int M = popBit / 2;
            vector<ll> dp2(1 << M);
            dp2[0] = 1;
            for(int mask=0; mask<(1 << M); mask++) {
                int idx = __builtin_popcount(mask);
                // winner[idx] vs loser[jdx]
                for(int jdx=0; jdx<M; jdx++) {
                    if(mask >> jdx & 1) continue;
                    int nmask = mask | (1 << jdx);

                    int A = winner[idx], B = loser[jdx];
                    if(match[A][B] == 1) {
                        dp2[nmask] += dp2[mask];
                    }
                }
            }

            ll addVal = dp2[(1 << M) - 1] * dp[bit];
            addVal *= (1 << M);
            // addVal *= frac[M];
            dp[par] += addVal;
        }while(next_permutation(perm.begin(), perm.end()));
    }

    /*
    for(int i=0; i<(1<<N); i++) {
        if(dp[i] == 0) continue;
        cout << "dp[" << bitset<N>(i).to_string() << "] = " << dp[i] << endl;
    }
    */

    for(int i=0; i<N; i++) {
        cout << dp[1 << i] / 2 << endl;
    }
    return 0;
}
0