結果

問題 No.1690 Power Grid
コンテスト
ユーザー hir355
提出日時 2021-09-24 22:32:10
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 151 ms / 3,000 ms
コード長 3,405 bytes
コンパイル時間 2,414 ms
コンパイル使用メモリ 199,828 KB
最終ジャッジ日時 2025-01-24 17:32:54
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

// #include <atcoder/modint>
#include <bits/stdc++.h>

using namespace std;
// using namespace atcoder;

constexpr long long INF_LL = 2000000000000000000LL;
constexpr int INF = 1000000000;
constexpr long long MOD = 998244353;

#define all(x) x.begin(), x.end()
#define REP(i, a, b) for(int i = a; i < b; i++)
#define rep(i, n) REP(i, 0, n)

typedef long long ll;
typedef pair<int, int> P;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<P> vp;
typedef vector<ll> vl;

int dx[4] = {0, -1, 0, 1};
int dy[4] = {1, 0, -1, 0};
int sign[2] = {1, -1};

template <class T> bool chmax(T &a, T b) {
    if(a < b) {
        a = b;
        return 1;
    }
    return 0;
}
template <class T> bool chmin(T &a, T b) {
    if(a > b) {
        a = b;
        return 1;
    }
    return 0;
}

ll modpow(ll a, ll b, ll m) {
    if(b == 0)
        return 1;
    ll t = modpow(a, b / 2, m);
    if(b & 1) {
        return (t * t % m) * a % m;
    } else {
        return t * t % m;
    }
}

template <class T>
T gcd(T m, T n) {
    if(n == 0) return m;
    return gcd(n, m % n);
}

struct edge {
    int to;
    ll cost;
    edge(int t, ll c) { to = t, cost = c; }
};

typedef vector<vector<edge>> graph;

// using mint = modint998244353;
// constexpr int MAX_COM = 100001;
// mint fac[MAX_COM], ifac[MAX_COM];
// void initfac() {
//     fac[0] = ifac[0] = 1;
//     REP(i, 1, MAX_COM) fac[i] = i * fac[i - 1];
//     REP(i, 1, MAX_COM) ifac[i] = 1 / fac[i];
// }
// mint nCr(int n, int r){
//     if(r < 0 || n < r) return 0;
//     return fac[n] * ifac[n - r] * ifac[r];
// }

// typedef ll S;
// S op(S x, S y){ return min(x, y); }
// S e(){ return INF_LL; }
// typedef ll F;
// S mapping(F f, S x){ 
//     if(f == -1) return x;
//     return f;
// }
// F composition(F f, F g){ 
//     if(f == -1) return g;
//     return f; 
// }
// F id(){ return -1; }

template< typename T >
void warshall_floyd(vector<vector< T >> &g, T INF) {
  for(int k = 0; k < g.size(); k++) {
    for(int i = 0; i < g.size(); i++) {
      for(int j = 0; j < g.size(); j++) {
        if(g[i][k] == INF || g[k][j] == INF) continue;
        g[i][j] = min(g[i][j], g[i][k] + g[k][j]);
      }
    }
  }
}

ll dp[262144];

void solve(){
    int n, m, k;
    cin >> n >> m >> k;
    rep(i, 1 << n) dp[i] = INF_LL;
    vi a(n);
    rep(i, n) cin >> a[i];
    vector<vl> mat(n, vl(n, INF_LL));
    rep(i, n) mat[i][i] = 0;
    rep(i, m){
        int x, y, z;
        cin >> x >> y >> z;
        mat[x - 1][y - 1] = z;
        mat[y - 1][x - 1] = z;
    }
    warshall_floyd(mat, INF_LL);
    dp[0] = 0;
    rep(i, 1 << n){
        rep(j, n){
            if(i & (1 << j)){
                continue;
            }
            ll d = INF_LL;
            rep(k, n){
                if(i & (1 << k)){
                    chmin(d, mat[j][k]);
                }
            }
            if(d == INF_LL) d = 0;
            chmin(dp[i | (1 << j)], dp[i] + d + a[j]);
        }
    }
    ll ans = INF_LL;
    rep(i, 1 << n){
        int cnt = 0;
        rep(j, n){
            if(i & (1 << j)){
                cnt++;
            }
        }
        if(cnt == k) chmin(ans, dp[i]);
    }
    cout << ans << endl;
}

int main(){
    cin.tie(0);
    std::ios_base::sync_with_stdio(false);
    std::cout << std::fixed << std::setprecision(16);

    // initfac();
    int t;
    t = 1;
    // cin >> t;
    while (t--)
    {
        solve();
    }
}
0