結果

問題 No.753 最強王者決定戦
コンテスト
ユーザー vjudge1
提出日時 2026-06-19 10:47:16
言語 C++14
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=c++14 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 475 ms / 1,000 ms
コード長 1,983 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,147 ms
コンパイル使用メモリ 181,440 KB
実行使用メモリ 13,696 KB
最終ジャッジ日時 2026-06-19 10:47:27
合計ジャッジ時間 4,102 ms
ジャッジサーバーID
（参考情報） 		judge1_0 / judge2_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 4
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll dp[5][1<<16][16];
int win[16][16];
int mat[16][16];
int main(){
    for(int i = 0;i < 16;i++){
        for(int j = 0;j < 16;j++){
            cin >> mat[i][j];
        }
    }
    for(int i = 0;i < 16;i++){
        for(int j = 0;j < 16;j++){
            if(i == j)win[i][j] = i;
            else{
                int a = min(i,j),b = max(i,j);
                win[i][j]= mat[a][b] == 1?a:b;
            }
        }
    }
    for(int i = 0;i < 16;i++){
        dp[0][1<<i][i] = 1;
    }
    for(int k = 1;k < 5;k++){
        int s=1<<k,h=s>>1;
        for(int m = 0;m < (1<<16);m++){
            if(__builtin_popcount(m)==s){
                for(int L = m;;L = (L-1)&m){
                    if(__builtin_popcount(L) == h){
                        int R = m^L;
                        for(int x = 0;x < 16;x++){
                            if(m>>x&1)
                                if(L>>x&1){
                                    ll lw = dp[k-1][L][x];
                                    ll rs = 0;
                                    for(int y=0;y<16;y++)
                                        if(R>>y&1&&win[x][y]==x)
                                            rs += dp[k-1][R][y];
                                    dp[k][m][x] += lw*rs;
                                }else{
                                    ll rw = dp[k-1][R][x];
                                    ll ls = 0;
                                    for(int y=0;y<16;y++)
                                        if(L>>y&1&&win[x][y]==x)
                                            ls += dp[k-1][L][y];
                                    dp[k][m][x] += rw*ls;
                                }
                        }
                    }
                    if(L == 0)break;
                }
            }
        }
    }
    int f=(1<<16)-1;
    for(int i = 0;i < 16;i++){
        cout << dp[4][f][i] <<'\n';
    }
    return 0;
}
0