結果
問題 |
No.1690 Power Grid
|
ユーザー |
![]() |
提出日時 | 2021-09-25 11:00:20 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,368 ms / 3,000 ms |
コード長 | 3,106 bytes |
コンパイル時間 | 2,669 ms |
コンパイル使用メモリ | 200,180 KB |
最終ジャッジ日時 | 2025-01-24 18:05:56 |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 25 |
ソースコード
#include <bits/stdc++.h> //#include<boost/multiprecision/cpp_int.hpp> //#include<boost/multiprecision/cpp_dec_float.hpp> // #include <atcoder/all> using namespace std; #define rep(i, a) for (int i = (int)0; i < (int)a; ++i) #define repl(i, a) for (long long i = (long long)0; i < (long long)a; ++i) #define rrep(i, a) for (int i = (int)a; i > -1; --i) #define REP(i, a, b) for (int i = (int)a; i < (int)b; ++i) #define REPL(i, a, b) for (long long i = (long long)a; i < (long long)b; ++i) #define RREP(i, a, b) for (int i = (int)a; i > b; --i) #define pb push_back #define eb emplace_back #define all(x) x.begin(), x.end() #define rall(x) x.rbegin(), x.rend() #define popcount __builtin_popcount #define popcountll __builtin_popcountll #define fi first #define se second using ll = long long; constexpr ll mod = 1e9 + 7; constexpr ll mod_998244353 = 998244353; constexpr ll INF = 1LL << 60; // #pragma GCC target("avx2") // #pragma GCC optimize("O3") // #pragma GCC optimize("unroll-loops") //using lll=boost::multiprecision::cpp_int; // using Double=boost::multiprecision::number<boost::multiprecision::cpp_dec_float<128>>;//仮数部が1024桁 template <class T> inline bool chmin(T &a, T b) { if (a > b) { a = b; return true; } return false; } template <class T> inline bool chmax(T &a, T b) { if (a < b) { a = b; return true; } return false; } ll mypow(ll x, ll n, const ll &p = -1) { // x^nをmodで割った余り if (p != -1) { x = (x % p + p) % p; } ll ret = 1; while (n > 0) { if (n & 1) { if (p != -1) ret = (ret * x) % p; else ret *= x; } if (p != -1) x = (x * x) % p; else x *= x; n >>= 1; } return ret; } struct myrand{ random_device seed; mt19937 mt; myrand():mt(seed()){} int operator()(int a,int b){//[a,b) uniform_int_distribution<int>dist(a,b-1); return dist(mt); } }; // using namespace atcoder; //------------------------ //------------------------ //------------------------ //------------------------ //------------------------ void solve() { int n,m,k; cin>>n>>m>>k; vector<int>a(n); rep(i,n)cin>>a[i]; vector<vector<ll>>dist(n,vector<ll>(n,INF)); rep(i,n)dist[i][i]=0; rep(i,m){ int x,y,z; cin>>x>>y>>z; x--;y--; dist[x][y]=dist[y][x]=z; } rep(k,n){ rep(i,n){ rep(j,n){ chmin(dist[i][j],dist[i][k]+dist[k][j]); } } } vector<vector<ll>>dp(n+1,vector<ll>(1<<18,INF)); dp[0][0]=0; rep(i,n){ rep(bit,1<<n){ if(dp[i][bit]==INF)continue; chmin(dp[i+1][bit],dp[i][bit]); rep(j,n){ if(bit&1<<j)continue; ll res=!bit?0:INF; rep(k,n){ if(bit&1<<k){ chmin(res,dist[j][k]); } } chmin(dp[i+1][bit|1<<j],dp[i][bit]+res+a[j]); } } } ll ans=INF; rep(bit,1<<n){ if(popcount(bit)==k){ chmin(ans,dp[n][bit]); } } cout<<ans<<"\n"; } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cout << fixed << setprecision(15); solve(); return 0; }