結果

問題 No.753 最強王者決定戦
コンテスト
ユーザー vjudge1
提出日時 2026-06-19 10:42:35
言語 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  
実行時間 460 ms / 1,000 ms
コード長 1,803 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,181 ms
コンパイル使用メモリ 181,468 KB
実行使用メモリ 13,696 KB
最終ジャッジ日時 2026-06-19 10:42:43
合計ジャッジ時間 4,237 ms
ジャッジサーバーID
（参考情報） 		judge1_1 / judge3_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 4
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
long long dp[5][1<<16][16];
int win[16][16];
int main(){
    int mat[16][16];
    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){
                                    long long lw=dp[k-1][L][x];
                                    long long 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{
                                    long long rw=dp[k-1][R][x];
                                    long long 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';
}
0