結果

問題 No.114 遠い未来
ユーザー 沙耶花
提出日時 2021-10-26 21:16:16
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3,707 ms / 5,000 ms
コード長 1,779 bytes
コンパイル時間 4,434 ms
コンパイル使用メモリ 265,336 KB
最終ジャッジ日時 2025-01-25 07:23:08
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <stdio.h>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf 1000000000
int main(){
    
    int N,M,T;
    cin>>N>>M>>T;
    
    vector dis(N,vector<int>(N,Inf));
    rep(i,M){
        int a,b,c;
        cin>>a>>b>>c;
        a--;b--;
        dis[a][b] = c;
        dis[b][a] = c;
    }
    rep(i,N)dis[i][i] = 0;
    
    rep(i,N){
        rep(j,N){
            rep(k,N){
                dis[j][k] = min(dis[j][k], dis[j][i]+dis[i][k]);
            }
        }
    }
    
    vector<int> v(T);
    rep(i,T){
        cin>>v[i];
        v[i]--;
    }
    sort(v.begin(),v.end());
    
    vector<int> cov;
    rep(i,N){
        if(!binary_search(v.begin(),v.end(),i))cov.push_back(i);
    }
    int ans = Inf;
    if(T<=15){
        vector dp(N,vector<int>(1<<T,Inf));
        rep(i,N){
            dp[i][0] = 0;
            
            rep(j,T){
                dp[i][1<<j] = dis[i][v[j]];
            }
        }
        rep(i,1<<T){
            if(__builtin_popcount(i)<=1)continue;
            for(int S = i&(i-1);S!=0;S = (S-1)&i){
                rep(j,N){
                    dp[j][i] = min(dp[j][i],dp[j][S] + dp[j][S^i]);
                }               
            }
            rep(j,N){
                rep(k,N){
                    dp[j][i] = min(dp[j][i],dp[k][i] + dis[k][j]);
                }
            }
        }
        rep(i,N)ans = min(ans,dp[i].back());
    }
    else{
        vector<pair<int,pair<int,int>>> E;
        rep(j,N){
            for(int k=j+1;k<N;k++){
                E.emplace_back(dis[j][k],make_pair(j,k));
            }
        }
        sort(E.begin(),E.end());
        vector<bool> f(N,false);
        rep(i,v.size())f[v[i]] = true;
        rep(i,1<<cov.size()){
            rep(j,cov.size()){
                if((i>>j)&1)f[cov[j]] = true;
            }
            
            dsu D(N);
            
            int temp = 0;
            rep(j,E.size()){
                int x = E[j].second.first,y = E[j].second.second;
                if(f[x]&&f[y]){
                    if(!D.same(x,y)){
                        D.merge(x,y);
                        temp += E[j].first;
                    }
                }
            }
            ans = min(ans,temp);
            rep(j,cov.size()){
                if((i>>j)&1)f[cov[j]] = false;
            }
        }
    }
    
    cout<<ans<<endl;
    
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0