結果

問題 No.1690 Power Grid
ユーザー tabae326
提出日時 2021-09-24 22:44:51
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 254 ms / 3,000 ms
コード長 3,425 bytes
コンパイル時間 2,752 ms
コンパイル使用メモリ 205,448 KB
最終ジャッジ日時 2025-01-24 17:37:18
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define rep(i, srt, end) for (long long i = (srt); i < (long long)(end); i++)
// Ref: https://qiita.com/ysuzuki19/items/d89057d65284ba1a16ac
#define dump(var)  do{std::cerr << #var << " : ";view(var);}while(0)
template<typename T> void view(T e){std::cerr << e << "\n";}
template<typename T> void view(const std::vector<T>& v){for(const auto& e : v){ std::cerr << e << " "; } std::cerr << "\n";}
template<typename T> void view(const std::vector<std::vector<T> >& vv){ std::cerr << "\n"; for(const auto& v : vv){ view(v); } }
template<typename T> void dump_cout(const T& v) { for(long long i = 0; i < v.size(); i++) std::cout << v[i] << (i == v.size()-1 ? "\n" : " "); }

namespace wf {
    // There is no edge with negative cost
    template<class T>
    void warshallFloyed(std::vector<std::vector<T>> &d) {
        int n = d.size();
        for(int i = 0; i < n; i++) d[i][i] = 0;
        for(int i = 0; i < n; i++) {
            for(int j = 0; j < n; j++) {
                for(int k = 0; k < n; k++) {
                    d[j][k] = min(d[j][k], d[j][i] + d[i][k]);
                }
            }
        }
    }
    // There are edges with negative cost.
    // If graph has negative cycle, return false
    template<class T>
    bool warshallFloyedNeg(std::vector<std::vector<T>> &d, T inf) {
        bool res = true;
        int n = d.size();
        auto chmin = [&](T &a, T b) {
            if(a == inf) a = b;
            else if(b == inf) ;
            else a = std::min(a, b);
        };
        auto addinf = [&](T a, T b) -> T {
            if(a == inf || b == inf) return inf;
            else return a + b;
        };
        for(int i = 0; i < n; i++) chmin(d[i][i], 0);
        for(int i = 0; i < n; i++) {
            for(int j = 0; j < n; j++) {
                for(int k = 0; k < n; k++) {
                    chmin(d[j][k], addinf(d[j][i], d[i][k]));
                }
            }
        }
        for(int i = 0; i < n; i++) if(d[i][i] < 0) res = false;
        return res;
    }
};

void solve() {
    ll n, m, k;
    cin >> n >> m >> k;
    vector<ll> a(n);
    rep(i, 0, n) cin >> a[i];
    const ll inf = 1e18;
    vector<vector<pair<ll,ll>>> G(n);
    vector dist(n, vector<ll>(n, inf));
    rep(i, 0, m) {
        ll x, y, z;
        cin >> x >> y >> z;
        x--; y--;
        G[x].push_back({y, z});
        G[y].push_back({x, z});
        dist[x][y] = z;
        dist[y][x] = z;
    }
    wf::warshallFloyed(dist);
    vector dp(n+1, vector<ll>(1<<n, -1));
    auto tsp = [&](auto tsp, ll cur, ll mask) -> ll {
        if(dp[cur][mask] != -1) return dp[cur][mask];
        if(cur == k) return 0;
        ll res = inf;
        rep(i, 0, n) {
            if(mask & (1<<i)) continue;
            if(cur != 0) {
                ll mn = inf;
                rep(j, 0, n) {
                    if((mask & (1<<j)) && i != j) {
                        mn = min(mn, dist[i][j]);
                    }
                }
                if(mn == inf) continue;
                res = min(res, tsp(tsp, cur+1, mask | (1<<i)) + mn + a[i]);
            } else {
                res = min(res, tsp(tsp, cur+1, mask | (1<<i)) + a[i]);
            }
        }
        return dp[cur][mask] = res;
    };
    cout << tsp(tsp, 0, 0) << endl;
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    solve();
    return 0;
}
0